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.  ,  ABSTRACT 

V. 

The  design  of  heat  pipes  involves  knowledge  of  phenomena 
such  as  surface  tension  forces,  wick  permeability*,  and  fluia 
vaporization  and  condensation.  Considerable  variability  in 
these  phenomena  has  been  observed  in  heat  pipe  experiments  v 
Thus,  a  probabilistic  design  model  for  predicting  heat  ;piv e 
heat  transfer  rate  .has-been  developed  taking  into;  considera¬ 
tion  uncertainty  in  the  prediction  of  the  above  phenomena. 

The  probabilistic  model  yields  a  mean,  a  standard  deviation, 
and  the  distribution  of  heat  transfer  rate  based  on  the 
means,  standard  deviations,  and  distributions  of  the  deoign 
parameters.  The  probabilistic  method  is  compared  to  experi¬ 
mental  data  from  heat  pipes  with  wire  mesh  wicks.  Mean  val¬ 
ues,  standard  deviations,  and  distributions  are  presented 
for  wick  permeability,  critical  radius,  area,  porosity,  tor¬ 
tuosity,  and  heat  transfer  rate.  A  technique  is  described 
for  making  wire  mesh  wicks.  The  probabilistic  design  model 
indicates  the  range  of  design  without  the  use  of  safety  -fac¬ 
tors.  ( 
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CHAPIER  1 


INTRODUCTION 


1.1  ;  Design  Methodology 


Optimization  and  reliability  of  design  have  been  the 


topics  of  growing  interest  to  engineers  in  recent  years.  A 
critical  review  of  present  design  methodology.-  that  based  on 
arbitrary  "ignorance  factors"  or  safety  factors,  is  under  way.’1' 
The  realization  that  design  parameters  are  usually  character¬ 


ized  by  some  statistical  distribution  of  values  rather  than 


by  a  single  value  indicates  that  probabilistic  methodology  is 


a  logical  alternative.  The  present  deterministic  (single  val¬ 
ued)  methods  are  special  cases  of  the  probabilistic  methodology 


when  the  parameter  variabilities  are  set  equal  to  zero.  Fig¬ 
ure  1.1  shows  a  comparison  between  the  conventional  determin¬ 


istic  design  and  probabilistic  design  methodologies.  In  actual 
physical  systems,  the  absence  of  design  parameter  variability 
is  indeed  a  rare  case. 


Figure  1.1.  Design  .Methodology. 


It  is  not  difficult  to  find  examples  of  design  parameter 
variation.  For  instance/  in  heat  transfer  we  describe  the 
rate  of  energy  radiation  from  a  body  as 


Q 

* 

A 


eobT- 


(1.1) 


where  Q  is  the  rate  of  energy  radiated,  T  is  the  absolute  tem¬ 
perature  of  the  body,  cr^is  a  universal  physical  constant-,  e 
is  a  property  of  an  emitting  surface,  and  A  is  the  surface 

2 

area  of  tne  emitting  body.  Values  of  £  are  quite  variable. 
Values  of  Ujj  and  T  may  depend  on  the  accuracy  of  the  instru¬ 
mentation  used.  The  functional  variation  of  Q  will  grow  as 
the  fourth  power  of  the  temperature  variation.  It  is  in¬ 
deed  optimistic  to  believe  that  the  calculated  Q  and  actual 
Q  are  identical.  Yet  that  is  what  the  deterministic  model 
implies.  The  objective  of  this  discussion  is  not  to  under¬ 
mine  the  deterministic  model  but  to  improve  upon  it.  The 
deterministic  model  has  been  an  engineering  tool  for  many 
years  and  has  worked  well.  The  probabilistic  model  uses 
the  deterministic  equation  to  find  the  mean  value  of  the 
design  result.  But  the  added  feature  of  the  probabilistic 
methodology  is  that  it  yields  a  statistical  distribution 
of  the  occurrence  of  the  design  result.  This  determines 
the  expected  range  of  the  design  result  or,  in  the  parti¬ 
cular  example  mentioned  (Equation  1.1) ,  the  probabilistic 
design  yields  a  bound  on  the  variation  of  Q.  We  are  now 
able  to  determine  a  range  in  which  the  actual  energy  ra¬ 
diated,  Q,  may  lie  instead  of  calculating  a  single  value 
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of  Q  and  hoping  that  it  is  similar  in  magnitude  to  the  actual 
Qi  In  many  cases,  an  agreement  of  plus  or  minus  30%  between 
the  calculated  result  and  the  actual  result  is  considered  to 
be  normal.  This  is  in  itself  an  admission  of  the  uncertainty 
in  modeling  natural  phenomenon.  The  probabilistic  approach  to 
design  gives  us  a  logical  measure  of  this  uncertainty. 

1.2  Design  of  Heat  Pipes 

The  application  of  deterministic  design  techniques  to  the 
design  of  heat  pipes  has  indicated  a  need  for  an  improved  tech¬ 
nique.  Uncertainties  in  wick  pore  size,  wick  permeability, 
surface  conditions,  liquid  inventory,  and  fluid  properties 
has  rendered  the  classical  deterministic  design  approaches 
inaccurate.  Probabilistic  design,  which  treats  the  design 
parameters  as  random  variables,  has  been  used  successfully  to 
describe  the  integrity  of  structural  components.1  Thus,  it 
is  logical  to  consider  the  use  of  probabilistic  techniques  in 
the  design  of  heat  pipes. 

3 

Haring  and  Greenkorn  used  a  statistical  distribution  to 
describe  the  pore  radius  in  a  porous  medium.  Figure  1,2  shows 
the  beta  distribution  which  was  used  to  model  the  uncertainty 
of  pore  radius  in  a  porous  medium.  Pore  radius  is  so  vari¬ 
able  that  deterministic  descriptions  cannot  be  made.  It  is 
virtually  impossible  to  randomly  choose  a  particular  pore  and 
accurately  calculate  its  radius.  The  probability  distribution 
enables  the  researcher  to  do  the  next  best  thing.  That  is, 
one  will  be  able  to  describe  a  range  of  radii  in  which  a  par¬ 
ticular  pore  radius  might  be  observed.  One  is  able  to  observe 
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Figure  1.2.  Beta  Distribution  Model  of  Pore 
Radius  in  a  Porous  Medium 


Figure  1.3.  Permeability  Distribution  of  100 
Mesh  Stainless  Steel  Screen 


the  range  of  variation  and  the  probabilities  of  residing  in 
certain  locations  of  the  range.  Also,  experiments  at  The 
University  of  New  Mexico  have  shown  the  wide  variability  of 
wick  permeability  data. ^ 

Figure  1.3  is  a  plot  of  a  continuous  distribution  approx¬ 
imation  of  experimental  permeability  data  and  illustrates  the 
high  variability  in  this  heat  pipe  design  variable.  The  mea¬ 
sured  permeability  of  wick  samples  taken  from  the  same  materi¬ 
al  is  not  exactly  the  same  even  though  great  care  may  be  taken 
to  uniformly  clean  and  assemble  each  sample  in  an  identical 
way.  Choosing  one  particular  value  of  wick  permeability  or 
ev*en  the  average  value  would  be  erroneous  in  light  of  the  high 
scatter  between  samples.  Therefore,  it  is  more  correct  to  de¬ 
scribe  wick  permeability  using  the  average  of  all  samples  as 
an  estimate  of  the  mean  permeability,  K,  and  the  standard  de¬ 
viation,  s^,  as  a  measure  of  scatter. 

Heat  pipe  operation  is  a  function  of  permeability,  sur¬ 
face  tension,  wetting  angle,  fluid  properties  and  geometrical 
configuration.®  If  any  of  these  parameters  are  distributed 
and,  therefore,  described  by  a  probability  distribution,  the 
resulting  heat  flux  capability  of  the  heat  pipe  will  be  dis¬ 
tributed  with  a  mean  heat  flux,  Q,  and  a  standard  deviation, 

S  ,  as  a  measure  of  variability. 

Plotting  of  typical  design  results,  as  in  Figure  1.4, 

c>  • 

shows  the  power  of  the,  probabilistic  design  method.  The  ‘5 
deterministic  method  gives  an  acceptable  answer  and  a 
factor  of  safety  while  the  probabilistic  method  indicates 
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a  finite  probability  of  failure.  One  minus  the  probability  * 
of  failure  (1-p)  is  the  probability  of  a  successful  design.  ' 
The  probability  of  obtaining  a  successful  design  is  usually 
termed  the  reliability 
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Figure  1.4.  Design  Methodology  Comparison 

The  usefulness  of  the  probabilistic  design  technique 
may  be  extended  to  heat  exchanger  design.  There  may  be  hun¬ 
dreds  of  heat  pipes  in  a  heat  exchanger  system,  each  having 
a  different  heat  transfer  capability.  Conventional  determin¬ 
istic  design  dictates  that  the  minimum  heat  transfer  capa¬ 
bility  be  used  as  characteristic  of  each  heat  pipe.  This  is 
essentially  a  "worst  case"  analysis.  Another  conventional 
design  technique  uses  the  mean  heat  transfer  rate  multiplied 
by  a  safety  factor  based  on  experience.  The  calculated  per¬ 
formance  of  the  heat  exchanger,  using  either  the  "worst  case" 
values  or  the  mean  values  with  a  safety  factor,  no  doubt  will 
give  a  design  safety  factor  which  works.  But  typically,  such 


a  design  is  very  conservative  and  the  overdesigned  heat  ex¬ 
changer  is  larger  and  more  expensive  than  necessary.  The 
probability  of  manufacturing  one  hundred  heat  pipes  that  have 
a  performance  decidedly  worse  than  average >  and,  therefore, “a 
"worst  case,"  may  be  very  small.  For  example,  if  the  proba¬ 
bility  of  occurrence  of  the  worst  case  is  .05,  then  the  prob¬ 
ability  of  manufacturing  a  heat  exchanger  with  100  "worst  case" 
heat  pipes  is  (.05)^^,  which  is  a  very  rare  event. 

The  purpose  of  this  research  project  is  to: 

1.  Measure  statistical  data  on  water  heat  pipes  using 
wire  mesh  wicks. 

2.  Develop  a  probabilistic  model  of  heat  pipe  operation. 

Based  on  the  literature,  it  appears  that  heat  pipe  design 

considerations  have  been  purely  deterministic.  The  publica- 
4 

tion  of  Holm'  indicates  some  realization  of  experimental  un¬ 
certainty  in  the  collection  of  data.  Holm  presents  his  data 
using  an  error  bound  but  does  not  mention  anything  about  re- 

5 

peatability  and  distribution  theory.  Phillips  presents  data 
on  nominal  pore  diameter  for  200  mesh  stainless  steel  screen. 

It  is  optimistic  to  hope  that  this  deterministic  pore  diameter 
value  is  representative  of  the  sample.  A  mean  pore  diameter, 
standard  deviation  and  distribution  will  better  represent  the 
physical  characteristics  of  this  wick.  Also,  it  has  been  ob¬ 
served  at  UNM  that  permeabilities  of  screen  samples  have  large 

standard  deviations  despite  many  efforts  to  obtain  uniform 

.  20 
samples. 


A  survey  of  the  literature  reveals  no  effort  to  account 
for  parameter  variability  in  the  field  of  heat  transfer.  Little 
information  is  available  on  the  statistical,  behavior  of  heat 
transfer  parameters  such  as  convection  coefficient,  conductiv¬ 
ity,  and  permeability.  The  author  has  not  yet  found  any  lit-' 
erature  which  applies  probabilistic  approaches  to  heat  trans¬ 
fer  design. 


CHAPTER  2 
THEORY 

To  develop  the  probabilistic  theory  for  heat  pipe  design, 
it  is  necessary  to  derive  the  deterministic  model  ca  which  the 
probabilistic  model  is  based.  The  classic  derivation  of  the 

O 

heat  pipe  design  equation  appears  in  the  literature.  In  this 
chapter,  the  classic  derivation  of  the  fully  saturated  wick 
model  is  described.  This  derivation  is  based  on  the  prin¬ 
ciples  of  conservation  of  momentum,  energy,  and  mass  for  a 
differential  element  in  the  wick.  The  fully  saturated  wick 
model  is'  later  modified  for  the  partially  saturated  wick  con¬ 
dition  and  will  be  used  as  the  basis  of  the  probabilistic  de¬ 
sign  equation. 

2.1  Deterministic  Model  of  Heat,  Pipe  Operation 

The  heat  pipe  to  be  considered  in  this1  work  is  cylindri¬ 
cal  in  shape  with  an  annular  wick,  as  shown  in  Figure  2,1. 

The  major  assumptions  used  in  the  development  of  the  deter¬ 
ministic  model  are? 

1.  The  system  is  treated  as  one-dimensional. 

2.  The  wick  is  fully  saturated  with  liquid. 

3.  The  heat  flux  is  uniform  over  the  evaporator  and 
condenser  sections. 

4.  The  liquid-vapor  interface  meniscus  can  be  charac¬ 
terized  by  one  radius  of  curvature. 

5.  The  liquid  and  vapor  are  at  the  same  temperature 

along  the  entire  length  of  the  pipe. 

6.  Wick  properties  are  isotropic. 
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Containment  Shell 


Figure  2.1.  Operating  Heat  Pipe 


7.  Differentials  raised  to  the  second  or  higher  powers 

are  neglected.  '  ij - 

8.  The  heat  pipe  is  wick  limited ,  that  is,  viscous 
pressure  drops  in  the  liquid  saturated  wick  are  so 
much  larger  than,  those  in  the  vapor  tfrat  the  vis¬ 
cous  pressure  drop  in  the  vapor,  is  negligible. 

9.  The  vapor  condensing  on  the  liquid-vapor  interface 
has  a  velocity  in  the  y  direction  only.  Therefore, 
there  is  no  contribution  to  momentum  changes  in  the 
x-  and  z  directions. 

A  differential  fluid  element  taken  in  the  condenser  sec¬ 
tion  of  the  heat  pipe  wick  may  be  used  for  the  analysis.  Ap¬ 
plying  conservation  of  mass  to  the  element  shown  in  Figure 
2.2,  the  relationship  between  liquid  and  vapor  flow  rates  is 


obtained 


+  mv  m£(x+dx) 

,«n> 


(2.1) 


Where 


=  <><ie(£wVu>! 

X 


(2.2) 


h  (x+dx)  =  pAe(Vw)(uft  +  3 GT  dx) 


(2.3) 


Combining  Equations  2.1,  2,2,  and  2.3,  and  solving  for  m 


V  "  PJ i  3jT  dxe(fwV 


(2.4) 


From  the  liquid-vapor  interface  of  the  differential  element 


o  u  Z  dx 
Fv  v  w 


(2.5) 


*where  tw  is  the  average  wick  thickness  as  would  be  measured 
by  a  micrometer. 
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Figure  2.2.  Momentum  Terms  for  a  Differential  Element  of  Wick 


Figure  2.3.  Energy  Terms  for  a  Differential  Element  of  Wick 


Application  of  conservation  of  momentum  in  the  x  direction 
to  the  element  of  Figure  2,2  yields 

SFx  =  Mx+dx  '  Mx  !2-6;> 


where  the  momentum  terms  are 


P  o  ?  — 

M  =  —  nQ  e  t  Z T 
x  g  Si  w  w 


P«  2  -  P p.  ^<uo> 

Mx+dx  =  r  U1  e  Vw  +  3T  -TE-  dX  £  Vw 

c  c 


(2.7) 


(2.8) 


where  no  contribution  is  made  to  momentum  in  the  x  direction 
by  the  vapor  as  it  enters  the  element  as  a  result  of  assump¬ 
tion  9 . 

The  force  terms  are  composed  of  the  capillary  pressure 
forces  as  described  by  the  Laplace-Young  equation 


Ap  =  o(i-  +  |~) 
R1  R2 


By  assumption  4,  this  equation  reduces  to 


(2.9) 


AP  = 


2a 

R 


(2.10) 


where  Ap  is  the  pressure  drop  across  the  liquid  vapor  inter¬ 
face  at  x.  The  pressure  forces  are 


F  =  (p  -  —-)  e  t  Z 
p  ktv  R'  w  w 
*x 


(2.11) 


^x+dx 


=  (pv  +  dpv  -  ^;dR 


)  e  t  z 
w  w 


(2.12) 
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Ol?  ' 


For  the  heat  pipe  with  an  annular  wick  with  large,  vapor  pas¬ 
sage,  a  low  density  vapor  and  low  vapor  flow  rate,  pv  is  con- 

<■ 

st ant  and  dpy  is  zero.  The  force  on  the  element  is  found 
by  the  summation  of  Equations  2.11  and  2.12  • 


Fpx  -  FP(x+dx)  -  -  e  Vw  20  a2  7  ~ 


(2.13) 


RdR  will  be  neglected  when  compared  to  R^ . 

Because  of  the  low  flow  rates  and  velocities  that  occur 
in  capillary  wicks,  the  flow  is  free  of  irtertial  effects,  is 
laminar,  is  described  by  Darcy's  law 


-  m,u. 
dp  _  & 

dx  Kp^, 


(2.14) 


and  viscous  forces  on  the  element  are  given  by 


F£  =gax  [e<Vw)l 


(2.15) 


The  net  viscous  force  on  the  element  is  obtained  by  combining 
Equations  2.2,  2.14,  and  2 . 15 


Ff  =  f-  V,  V*  dx 


(2.16) 


If  the  heat  pipe  operates  in  a  gravitational  field  (or 
acceleration  field) ,  another  force  term  must  be  added  which 
is  the  weight  component  of  the  fluid  element  along  the  direc¬ 
tion  of  flow 


F  =  2— 


P g  e  t  Z  cost  dx 


(2.17) 
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The  gravity  component  force  may  be  ulus  or  minus  depending  on 
the  orientation  of  the  heat  pipe.  In  this  analysis,  as  shown 
in  Figure  2.1,  the  heat  pipe  is  oriented  so  that  the  evapora¬ 
tor  is  above  the  condenser  and  the  gravity  force  opposes  the 
liquid  flow.  The  summation  of  forces  on  the  element  is 


ZFx  FP(x)  Fp(xtdx)  “f  ~g 


(2.18) 


Combining  Equations  2.6,  2.7,  2.8,  2.13,  2.16,  2.17,  and  2.18, 
we  obtain  the  differential  momentum  equation  for  the  fluid 
element 


.  2o  S|  -  |  p  u  dx  -  |-  p  cos*  dx 
R  yc 


"t  a<V2 
gc  ^ 


dx 


(2.19) 


Conservation  of  energy  must  also  be  satisfied.  Referring  to 
Figure  2.3,  neglecting  kinetic  energy  effects  the  energy  terms 
are 


Energy  in 


Q  =  m  h 
v  v  v 


dT 

-  k.(t  z  )~ 
£  w  w  dx 


Energy  out 

.  a  (m  h  ) 

®£(x+dx)  ^SpS,  +  dx  ^ 

rdT  d2T  -j 

-  Ww>  1st  +  -T-  *9 

ax 


Qq  (heat  removed  externally) 
Wg  (work  against  gravity) 

(2.20) 


From  assumption  5,  the  conduction  terms  are  negligible  and 
the  energy  equation  becomes, 


(2.21) 


^£(xidx)  '  Qo  ~  **g 


replacing  a  with  Sqaation  2.4 


dr 


da£ 

d3T  ^  e  V5* 


(2.22) 


cocbining  Sanations  2.21  and  2.22, 


dG£  dtn£h£) 


etzdx  et*  (he 
W  W  W  c? 


and  simplifying  the  derivative 


dx  **£^£  dr  u£^£  dr 


(2.23) 


^  °o  .  »q  0 

^  eVw  "  %  ** 


(2.24) 


Fros  assumption  5,  arial  temperature  gradients  - 


are  negligible  and  assuming  axial  pressure  gradients  are 
snail  for  the  water  heat  oioe 


£ 

—  ~  0 


(2.25) 


Combining  Equations  2.24  and  2.25 


pj\  -  h£)  g-i  ^ 

*Vwa*  etwzw*= 


where 


K  0  _ 

V7  =  +  — 2__  cos  Aet  Z  dx 
g  J  g_  v  w  w 


(2.26) 


(2.27) 


Q  =  +  q  dx  Z 
o  “  w 


(2.28) 


h  -  h.  =  h_ 
v  £  fg 


(2.29) 
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(2.30) 


Combining  Equations  2.26,  2.27,  ara*  2.28 


GO 


Pai^a 


p  __JL - SL_  JLJl  2—  cos*  =  0 

rg  dx  —  -■  a 


n  simplifying 


et. 


du 


;£ 


G 


eViafg 


_  5L_  u 

g  Jh_  £ 
“c  rc 


(2.31) 


'c  rg 

Solving  for  the  average  velocity  and  integrating  over  the 


length,  of  the  condenser,  the  energy  equation  be^otnes 


U£  = 


g  x  g  cos®  f 

-2-—  J  uzax 


hfgeV* 


g  dh, 
c  rg 


(2.32) 


Integrating  to  a  position  x  along  the  condenser  and  rearrang¬ 
ing  terns 


q  x  g  x 

£ -  =  (1  j.  -2—  _  -) 

u5n-  et  pff  yc  rg 


(2.33) 


‘rfg- 

The  gravity  tern  inside  the  bracket  is  much  smaller  than  1 
and  can  be  neglected,  resulting  in  the  following  simplified 
form  of  the  energy  equation 
q _ x 


U£  * 


h_  et  p0 
rg  w 


(2;34) 


Combining  Equations  2.34  and  2,19  to  form  the  energy-continuity- 
momentum  integral  equation 


/R 

V° 

■  r 


2a 


2£ 


dR 

R2 

2 


-  r  ^ 

Jn  K  h-  t 


u„ 


?  p< 


xdx 


,2  <£-2 
gc  fg  pJi&  fcw 


0  “fg^w 
xdx 


S-P£cos$dx 


(2.35) 


(at  x  =  0  R-*-co  at  x  =  x  R 

max 


mm 


=  r) 
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Sis  sass  slew  disfcribat5.ca  resulting  from  the  of 

uniform  heat,  flex  in  tire  evaporator  m^sneor  jg  shown 
in  Figure  2.4. 


Figure  2.4.  Mass  Flow  Distribution 
Equation  2.35  is  valid  for  the  condenser  region  only  and  new 

a  relation  in  the  evaporator  region  is  needed.  To  link  the 
energy  and  momentum  analysis  between  the  condenser  and  evapor¬ 


ator,  we  will  use  the  following  boundary  conditions 


u 


g  I, 

”c  c 


max. 


hfgVpi 


=  u 


G  L 

-e  e 


max. 


h_  t  ep. 
rg  w  £ 


(2.36) 


u£  = 


a  I, 

- 2=_£_  -  (i  -  x) 

*W^ep0  e  h,  t  ep. 


0<x<I» 


f  gw£ 


fg~w  2, 


(2.37) 


qe  = 


q  L 
c  c 


(2.38) 


g  L 
c  c 


k-  tep. 


•1  + 


(5£  -  £ 


n 


fa 


(2.39) 


Equation  2.39  is  the  link  between  the  evaporator  and  condenser 
velocities  and  can  be  combined  with  Equation  2.35,  resulting 
in  the  energy- continuity-momentum  integral  equation  for  the 
evaporator.  The  same  but  more  simplified  analysis  may  be 
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( 


applied  to  tis  aaiabstic  sectic-n  since  there  the  neat  flaz 
term  is  2ero  sad  the  only  forces  entering  into  the  equation 
are  field  friction  and  capillary  pressure  rise.  Addition  of 


the  three  mar.  ei>. tra-cc-atinui tg -energy  equations  yields 


psh  d  f^Si' 


P: 


\/2 


5,,  ^ y=l 

^  ■—  .is  “  1 - . 


V2  '  ■“« 


"i  !  j‘r 

rl\  e 

2 


-  M 

-C  I 

(2.40)  / 


_  Vfg0  feA  [2_  _  9&  A 

V;  \PJ  \re  -Sc  / 


A 


(2.41) 


Equation  2.41  is  the  basic  deterministic  maximum  heat  trans¬ 
fer  rate  equation  for  the  heat  pipe  considered  in  this  work. 


The  tens  {ph o/u)  is  ccmonly  referred  to  as  the  liquid  trans¬ 
port  factor,  which  is  a  group  of  fluid  properties.  The  term 
(Eh/L* )  is  a  group  of  wick  properties  where  I»*  is  the  affec¬ 
tive  flow  length  or  average  flow  length  of  a  given  liquid  par¬ 
ticle  in  the  wick.  1*  is  obtained  from  the  integration  of 
Darcy’s  law.  Equation  2.14,  according  to  the  mass  flow  dis¬ 
tribution  of  Figure  2.4.  The  mass  flow  equations  for  the 
heat  pipe  considered,  in  this  work  are  given  by 


m  (x)  =  m  j—  (0  <_  x  <  L  ) 

e  .  se  e 

ma(x)  =  m  (Le  <  x  <  Le  +  La) 

/  (x  -  L  -  L  )\ 

mc(x)  =  m  (l - - - — J  (Le  +  La  <  x  <  L)  (2.42) 


The  mechanics  of  the  integration  of  Darcy's  law  are  shown  in 
Eguation  2.43  resulting  in  the  expression  for  L'. 


19 


ne(x)  OX  -r 


lL(X>  dx 

cL 


E„(x)  dx 


(2.43) 


PjOL' 

■& 


(2.44) 


L  L 

I»*  =  2~  +  I»a  2“  (f°r  fully  saturated  wick)  (2.45) 

Equation  2.45  appears  as  part  of  Equation  2.40  as  a  re¬ 
sult  of  the  combination  and  integration  of  the  momentum- 
energy-ccntinuity  equation.  The  reason  for  the  reiteration 
of  the  derivation  of  Equation  2.45  is  that  we  will  modify 
this  integration  step  to  obtain  the  design  equation  for  the 
partially  saturated  wick  condition. 

2.2  Deterministic  Model  of  Heat  Pipe  Operation  with  Par¬ 
tially  Saturated  Wick 

One  of  the  major  assumptions  in  the  preced.ing  analysis 
was  assumption  2  which  stated  that  the  heat  r-’-pe  wick  was 
fully  saturated.  Many  times  this  is  not  the  case.  The 
liquid-vapor  interface  may  recede  into  the  wick  result¬ 
ing  in  performance  different  from  that  predicted  by  Equation 
2.40.  An  attempt  to  modify  Equation  2.40  for  the  partially 
saturated  wick  condition  follows. 

Modifying  assumptions  for  a  partially  saturated  wick 
analysis  are: 
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1.  Capillary  force  properties  are  non-uniform  across 
the  wick  as  the  result  of  the  variation  of  critical 
capillary  radius  in  tne  wick.  Critical  radius  will 
be  a  function  of  t;  therefore,  critical  radius  will 
decrease  in  the  direction  towards  the  heat  pipe  wall. 

2.  Permeability  is  assumed  uniform  r^ross  the  wick. 

3.  Uniform  heat  flux  in  the  evaporator  and  condenser 
is  maintained  daring  the  partially  saturated  (de- 
s at urate a)  mode  of  operation. 

4.  The  effect  of  the  desaturation  is  to  increase  vis¬ 
cous  liquid  flow  losses  in  the  wick  due  to  higher 
liquid  velocities  in  the  desaturated  wick.  The 
higher  liquid  velocities  are  required  to  maintain 
the  mass  flow  in  a  smaller  area. 

5.  A  heat  pipe  functioning  at  non-equilibrium  during 
the  desaturation  or  recession  process  will  re¬ 
establish  equilibrium  only  if  the  receding  liquid 
vapor  interface  encounters  a  smaller  pore  size  re¬ 
sulting  in  sufficient  additional  capillary  force. 

Figure  2.5  graphically  illustrates  the  assumptions  of 
the  recession  or  desaturation  mode  of  heat  pipe  operation. 
Figure  2.5(a)  shows  the  mass  flow  distribution  for  the  satu¬ 
rated  and  partially  saturated  mode  of  operation.  The  shape 
of  the  distribution  is  a  result  of  the  assumption  that  the 
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2.5.  Mass  Flow  Rate,  Liquid  Recession  Distribution 
and  Liquid  Pressure  Distribution  for  Fully 
Saturated  and  Partially  Saturated  Wicks 


heat  flux  in  the  condenser  and  evaporator  are  uniform.  Figure 
2.5(b)  shows  the  liquid- vapor  interface  for  the  fully  satu¬ 
rated  and  partially  saturated  modes  of  operation.  The  fully 
saturated  condition  of  Figure  2.5(b)  indicates  that  there  is 
sufficient  capillary  force  at  the  upper  t  level  in  the  wick 
-to  sustain  the  frictional  and  gravitational  losses  of  the  sys¬ 
tem.  If  this  balance  cannot  be  maintained  by  the  capillary 
forces,  liquid  will  be  depleted  due  to  evaporation  in  the 
evaporator  more  quickly  than  it  can  be  restored  by  the  capil¬ 
lary  forces  at  the  given  liquid- vapor  interface.  The  result 
is  a  non-equilibrium  condition  in  which  the  liquid-vapor  in¬ 
terface  recedes  into  the  wick.  Now  chat  the  fluid  is  receding, 
there  are  two  possible  outcomes: 

1.  The  fluid  will  recede  until  it  encounters  a  stronger 
capillary  force  which  will  reestablish  an  equilib¬ 
rium  condition. 

2.  The  fluid  will  recede  until  it  encounters  the  heat 
pipe  wall  resulting  in  the  dryout  and  failure  of 
the  heat  pipe. 

For  the  analysis  of  recession,  we  will  assume  that  the  liquid- 
vapor  interface  will  reestablish  an  equilibrium  condition  at 
some  level  in  the  wick.  Returning  to  Figure  2.5(b),  we  see  a 
hypothetical  liquid-vapor  interface  distribution  for  a  par¬ 
tially  saturated  wick.  The  t2  level  in  the  wick  is  assumed 
to  have  the  highest  capillary  force  capability  and,  therefore, 
the  smallest  critical  radius.  The  t^  level  has  a  higher  cap¬ 
illary  force  capability  than  the  -t  level  but  a  smaller  capa¬ 
bility  than  the  €2  level. 
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Referring  to  Figure  2.5  (b  and  c) ,  the  mechanics  of  the 
recession  will  be  described.  Starting  at  the  condenser  end 
of  the  heat  *• ipe ,  we  observe  the  liquid  pressure  and  liquid 
level  (position  of  vapor-liquid  interface) .  Koving  from  right 
to  left  in  the  condenser  region,  the  liquid  pressure  begines 
to  decrease  due  to  gravity  and  viscous  forces.  This  pressure 
loss  can  still  be  maintained,  however,  by  the  capillary  pres¬ 
sure  at  the  t  level  in  the  wick.  Moving  into  the  adiabatic 
w  ~ 

region,  the  pressure  loss  in  the  liquid  has  exceeded  the  capil¬ 
lary  pressure  capability  of  the  t  level  pores  and  the  fluid 
recedes.  The  t1  level  pores  haVe . sufficient  capillary  force 
to  maintain  the  liquid  pressure  loss  and  the  liquid  level  re¬ 
mains  at  t^.  Moving  into  the  evaporator  region  with  ever  in- 
creasing  liquid  pressure  loss,  it  is  observed  that  the  t1 
level  pores  cannot  sustain  any  more  pressure  loss  and  the  li.- 
quid  level  recedes  to  the  next  level,  t 2.  At  the  t2  level, 
capillary  forces  are  sufficient  to  maintain  the  viscous  and 
gravitational  losses  of  the  system.  The  t2  level  has  the 
highest  capillary  pressure  capability  and  any  additional  re¬ 
cession  will  result  in  burnout.  The  difficulty  in  this  anal¬ 
ysis  is  the  determination  of  the  capillary  properties  at  dif¬ 
ferent  levels  if  such  a  situation  exists.  In  this  work,  the 
above  theory  is  applied  to  wire  mesh  wicks  whose  capillary 
properties  can  be  determined  at  the  various  levels.  This 
may  be  a  difficult  matter  for  powder  meta?  or  other  types  of 


wicks. 


The  preceding  analysis  indicated  that  thd  liquid-svapor 
interface  recedes  into  the  wick,  resulting  i‘u,  a  decreased  flow 
area  and  increased  liquid  pressure  loss.,  i*he  only  difference 
this  introduces  into  the  combined  energy-continuity-raomentum 
equation,  2.35,  is  the  thickness  of  the  liquid  element.  All 
other  terms  of  Equation  2.35  must  be  maintained  at  their  pres¬ 
ent  values,  consistent  with  Figure  2.5(a),  except  the  liquid 
pressure  loss  term  given  by  Equation  2.1<i,  Darcy's  law,  which 
will  be  used  to  form  the  integral  equation  for  the  pressure 
loss  in  the  liquid.  For  the  recessed  liquid-vapor  interface. 
Equation  2.45  does  not  apply  because  the  liquid  flow  area  is 
not  constant  with  length.  This  is  shown  in  Figure  2.6 


Figure  2.6.  Liquid-Vapor  Interface  Distribution 


Figure  2.6  is  a  redrawing  of  Figure  2.4(b)  and  will  be  used 
to  determine  the  total  liquid  pressure  loss  for  the  recessed 
liquid  level  condition.  Application  of  Equations  2.14  and 
2.42  to  Figure  2.6  the  total  pressure  loss  is 
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Simplifying,  we  obtain 
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and  the  resulting  model  for  heat  pipe  operation  in  the  de- 
saturated  or  recessed  wick  condition  is 
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(2.49) 


The  only  difference  between  Equation  2.49  and  2.41  is  the  L'(A)r 

term  (Darcy  flow  length)  and  the  critical  radius,  r  .  In 

c 

Equation  2.49,  we  use  the  L' (A^.  to  account  for  the  additional 
pressure  loss  due  to  the  recession  and  rfr  to  account  for  the 
increased  capillary  force  at  the  fully  recessed  level,  there¬ 
fore,  level  t 2  in  Figure  2.5.  Figure  2.6  may  be  interpreted 
as  a  general  case  of  fluid  recession  and  does  not  imply  that 
there  are  always  three  distinct  levels  of  constant  pore  size, 
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r.  The  integration  of  Equation  2 . 46  may  be  applied  to  any 
amount  of  recession  "steps"  as  long  as  one  integrates  between 

the  discontinuities.  The  discontinuities  of  Figure  2.6  are 

'0 

at  l>e/  l»r,  and  Le  +  La* 

In  the  design  of  heat  pipes  operating  under  the  desatu- 
rated  wick  condition,  the  recession  lengths,  therefore,  Lfr 
and  Lr  of  Figure  2.6,  must  be  determines  to  jlhsert  into  Equa¬ 
tion  2.48.  Also,  one  must  know  the  variation  of  r  as  a  func- 

c 

tion  of  the  t  dimension  of  Figure  2.6.  The  rc  variation  may 
be  theoretically  hypothesized  or  experimentally  determined. 

In  this  work,  r  as  a  function  of  t  is  determined  experiment- 
ally.  Given  that  the  rc  variation  is  known,  one  can  calculate 
the  recession  lengths  by  plotting  the  pressure  distribution 
along  the  wick  and  observing  where  the  liquid  pressure  loss 
exceeds  the  capillary  pressure  rise  of  a  given  pore.  The 
pressure  loss  in  the  wick  is  a  combination  of  viscous  and 
gravity  effects  and  is  given  by 


m(x)uA  PA9H 
KA(x)  pj,  ~  c£  L 


(2. 50) 


Integration  of  Equation  2.50  results  in  the  pressure  distri 
bution  along  the  wick  of  a  heat  pipe  shown  in  Figure  2.7. 
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If  one  has  no  excess  liquid  in  a  pipe  and  neglects  the  effect 
of  the  diameter  of  the  pipe,  Equation  2.51  becomes 
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where  m (x) 


is  defined  by 


m  (x)  xdx  _ 

A(x)  "  gQ 

Equation  2.42. 


cos<*> 


(2.52) 


Equation  2.52  is  difficult  to  solve  in  closed  form  be¬ 
cause  the  area  term  is.  a  function  of  x.  A  trial  and  error 
solution  is  proposed  using  Equations  2.49,  2.50,  2.52,  and 
2.53. 


m_ 


max 


(2.53) 


First,  compute  Q^ax  using  Equation  2.49,  initially  setting  L*  (A)^_ 
equal  to  L'/A.  Second,,  compute  m  using  Equation  2.53.  Then, 
using  Equation  2.52,  start  plotting  p  ,  beginning  at  the  con- 

A 

denser  end  at  x  equal  o.  Increment,  using  small  x  until 
the  term  (p  -  p  )  exceeds  the  capillary  pressure  rise  of  the 
pores  at  the  inside  diameter  of  the  wick.  At  this  point,  set 
A(x)  equal  to  the  new  and  smaller  recessed  area  as  a  result 
of  the  fluid  recession  and  record  the  x  position.  Referring 
to  Figure  2.6,  the  new  area  will  be  Ar%  Continue  to  plot 
pressure  along  the  entire  length  of  the  pipe  in  like  manner, 
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always  observing  tbs  mass  flow  distribution  of  Figure  2.5  fa)  - 
Fron  this  pressure  plot,  ore  can  calculate  the  recessed  lengths 
or  the  wick  ana,  therefore,  calculate  L’  C~)  r-  She  new  L'  {?-) 
iron  Equation  2.48  is  inserted  into  Equation  2.49  and  the  pro¬ 
cedure  is  repeated  until  the  O  value  converges,  which  usu- 

js ax 

ally  takes  about  four  iterations.  A  computer  program  was 
written  to  perform  this  computation  and  is  discussed  in  Chan¬ 
ter  6. 

2.3  The  Probabilistic  Z-lodel  of  Heat  Pipe  Operation 

Many  of  the  design  parameters  of  Equation  2.49  are  ex¬ 
tremely  variable.  Experiments  show  that  the  wicking  proper¬ 
ties,  K,  A,  and  rc  may  vary  plus  or  minus  30%  of  the  absolute 
mean  value.  The  reasons  for  the  variability  or  uncertainty 
in  the  determination  of  these  variables  result  from  variability 
of  materials,  variability  of  manufacturing,  and,  to  some  ex¬ 
tent,  experimental  measurement  variability.  Generally,  the 
uncertainty  in  these  design  variables  is  too  large  to  be  ne¬ 
glected.  The  probabilistic  model  of  heat  pipe  operation  in¬ 
corporates  the  variability  of  the  design  parameters  so  that 
the  variability  of  the  design  result,  ,  can  be  determined. 

Given  that  the  design  variables  of  a  particular  system 
are  described  by  a  distribution  of  values  rather  than  a  single 
deterministic  point,  we  may  be  able  to  fit  certain  functions 
to  data  to  describe  the  variability  of  design  variables. 

These  functions  are  called  probability  distributions.  The 
permeability,  K,  may  be  described  by  a  certain  probability 

distribution  while  the  critical  radius,  r  ,  may  be  described 

c 

by  another. 
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1.  Q,  K,  r  ,  3,  2.* '{&>  ,  and  A  are  assnsfi  tc  bs  irJgegehsest 

S  2T 

random  variables  described  Isy  a  csntinsgas  probabil¬ 
ity  distribution. 

2.  Sis  ccatinucgg  probability  d?  striSatioas  will  be 
described  by  two  parameters,  tbs  mean  and  standard 
deviation  . 

3-  all  otter  variables  will  be  considered  deterministic 
since  their  variabilities  are  comparatively  snail. 

The  notation  for  a  given  randan  variable  using  these  assump¬ 
tions  will  be  (u  ,  o  )  where  zz  is  the  random  variable,  n 

zz  zz 

is  the  mean  and  a  is  the  standard  deviation.  The  oercea- 

zz 

bility  variate  pair  being  (u-,,  o^)  . 

The  probabilistic  design  equation  is  the  deterministic 
design  equation  with  the  random  variate  pairs  inserted.  The 
probabilistic  design  equation  is 
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(2.54) 

Equation  2.54  represents  the  functional  relationship  be¬ 
tween  the  random  variates  of  the  probabilistic  design  model. 

Many  times,  the  distribution  parameters  (the  mean  and 
standard  deviation  in  this  analysis)  must  be  estimated  from 
experimental  data.  The  consistent  and  unbiased  estimator  for 

the  mean  of  a  random  variable,  x,  is  defined  by  Miller  and 
17 

Freund  as 
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Since  all  the  rand  eg  variates  gill  fee  det-era.in.sd  using  ezperi- 
aegtal  data,  their  parameter  gill  be  estimated  and  Equation 
2.54  for  tee  case  of  partially  saturated  wick  be  cores 
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The  solution  of  Equation  2.57  is  the  random  variate  (Q,S^)  , 
which  may  be  translated  into  a  probability  distribution  func¬ 
tion  such  as  the  normal  distribution.  Figure  2.8  shows  a 
hypothetical  distribution  of  maximum  heat  transfer  rate  for 
a  given  heat  pipe  design.  The  designer  may  determine  the 
variability  of  the  design  and  decide  if  it  meets  or  exceeds 
the  specified  minimum  value,  Qg.  The  area  to  the  right  of 
Qs  indicates  the  probability  that  a  heat  pipe,  with  the  given 
variable  properties,  will  exceed  Q  .  The  design  engineer  can 
adjust  the  ^ax  distribution  by  changing  the  distribution  of 
the  heat  pipe  properties  so  that  a  very  small  portion  of  the 
distribution  lies  below  Q  in  the  failure  region. 

o 
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Figure  2.8.  Distribution  of  Heat  Flew  Capability 
of  a  Heat  Pipe  whose  Design  Param¬ 
eters  Are  Assured  Hasdoa 

Solution  techniques  sill  be  discussed  in  depth  in  Chap¬ 
ter  6  for  Equation  2.57. 

She  basic  argument  behind  the  probabilistic  design  ap¬ 
proach  as  that  we  can  obtain  a  quantitative  measure  of  the 
uncertainty  of  the  system  performing  as  required.  Many  times 
we  may  obtain  a  set  of  measurements  that,  when  inserted  into 
the  deterministic  equation  (Equation  2.49,  or  any  equation 
for  that  matter) ,  give  a  vastly  different  result  from  reality. 
She  conclusion  may  be  that  the  particular  deterministic  equa¬ 
tion  is  poor,  but  this  may  not  be  a  proper  conclusion.  One 
must  look  at  the  distribution  of  functional  solution  variables 
to  determine  the  validity  of  the  model.  It  may  be  highly 
possible  to  obtain  a  calculated  result  relatively  far  from 
the  observed  mean  result  yet  still  lie  in  the  range  of  the 
distribution. 


Chapters  6  and  7 -gill  analyze  the  probabilistic  design 

saethodologv  in  depth  end  comparisons  will  he  zxade  with  actnal 

test  data.  In  Chapter  3,  the  details  of  experimental  tech- 

niqnas  for  the  saasaresent  of  our  design  variables,  K,  a,  E>*, 

r  ,  and  S,  are  Presented, 
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The  heat,  pipe  assign  used  as  a  basis  for  obtaining  data 
on  a  typical  wire  mesh  wick  beat  pipe  is  shown  in  Figure  3.1. 
The  material  used  in  the  construction  of  the  wicks  and  pipe  was 
204  stainless  steel.  The  dimensions  shown  in  Figure  3.1  were 
chosen  as  typical  for  a  wire  mesh  heat  pipe  of  this  design  with 
an  a*-  wick.  The  active  length  of  the  heat  pipe  was  22  in¬ 

ches  and  the  evaporator  was  13.5  inches.  The  evaporator  was 
chosen  to  be  large  (a  significant  portion  of  the  pipe)  to  in¬ 
sure  low  radial  neat  transfer  and  therefore  little  chance  of 
vapor  blockage  in  the  wick  as  a  result  of  high  radial  heat 
fluxes.  The  inside  diameter  of  the  heat  pipe  was  .743  inches 
and  the  outside  diameter  was  13/16  inches.  Thermocouples 
were  placed  along  the  outside  surface  of  the  heat  pipe  to 
monitor  axial  temperature  loss.  The  end  caps  on  the  heat 
pipe  were  removable  to  facilitate  the  testing  of  many  differ¬ 
ent  wicks  without  the  additional  labor  and  expense  of  construc¬ 
ting  a  new  heat  pipe  for  each  wick.  Vapor  thermo couples  were 
attached  to  the  removable  end  caps  to  monitor  vapor  tempera¬ 
ture.  Four  different  sizes  of  wire  mesh  were  used  to  obtain 
four  variations  in  the  design  equation. 

Figure  3.2  shows  the  stainless  steel  heat  pipe  with  end 
caps.  The  thermocouple  wiring  harness  was  attached  and  bonding 
cement  was  applied  to  the  thermocouple  junctions.  Figure  3.3 
shows  the  same  heat  pipe  of  Figure  3.2  after  insulation,  heater 
wire,  and  calorimeter  have  been  installed.  Excessive  insulation 
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was  added  to  the  evaporator  to  ggxinize  the  input  heater  effec¬ 
tiveness.  Figure  3.4  shows  the  experimental  heat  pipe  mounted 
and  ready  for  testing.  The  heat  pipe  is  tested  with  the  evap¬ 
orator  higher  than  the  condenser  so  that  the  wick  limited  con¬ 
dition  can  be  reached  without  excessive  heat  flux.  The  appara¬ 
tus  used  to  measure  the  gravity  effect  is  shown  in  Figure  3.5. 

A  sliding  probe  two  feet  in  length  is  inserted  into  the  heat 
pipe  until  electrical  contact  is  made  with  the  excess  working 
fluid  in  the  bottom  of  the  pipe.  Contact  is  indicated  by  a 
reading  on  a  microammefcer.  The  length  of  the  probe,  Lp,  is 
measured  and  the  pressure  drop  due  to  gravity,  Ap^,  is  cal¬ 
culated  using  Equation  3.1.  Figure  3.6  shows  the  wire  mesh 
wick  being  inserted  into  a  heat  pipe.  The  wicks  were  manu¬ 
factured  to  fit  as  tight  as  possible  and  were  inserted  through 
the  condenser  end  and  pressed  tightly  into  the  evaporator. 
Figure  3.7  shows  the  heat  pipe  under  test.  Figure  3.24  is 
a  simplified  drawing  of  the  heat  pipe  operating  at  steady 

Ap£g  =  LPisin0  +  rpi 

state  and  temperature  recorder  monitoring  temperature  distri¬ 
butions  along  the  pipe  and  the  calorimeter  temperature  rise. 

Figure  3.9  shows  the  technique  used  to  experimentally 
determine  wick  permeability.  A  balloon  approximately  two 
feet  in  length  was  inserted  into  the  vapor  cavity  of  the 
heat  pipe  and  pressurized  to  50  psi.  This  pressure  was 
chosen  to  ensure  a  sufficient  force  to  press  the  balloon 
against  the  wick  as  shown  in  Figure  3.26.  The  objective  of 
the  measurement  is  to  force  fluid  through  the  wick  in  a 
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manner  similar  to  the  actual  operation  of  the  heat  pipe.  Figure 
3.26  shews  the  small  error  that  will  result  in  the  permeability 
test  due  to  the  inability  of  the  balloon  to  seal  off  a  small  area 


near  the  seam.  This  error  is  considered  insignificant  since  that 
seam  void  carries  a  small  amount  of  fluid  by  the  meniscus 
shown  in  Figure  3.27.  A  pressure  head  of  six  inches  of  water 
was  used  to  drive  the  fluid  through  the  wick.  The  liquid  flow 
velocity  resulting  from  this  driving  pressure  was  approximately 
.015  ft/sec  (.015  ft/sec  was  calculated  to  be  the  maximum  flow 
encountered  during  operation  of  these  heat  pipes  and  is  in  the 
Darcy  flow  regime)  and  the  pressure  loss  was  taken  across  the 
wick  structure  only.  Permeability  was  calculated  using  Equa¬ 
tion  2.14.  Figure  3.11  shows  the  preparation  for  a  permeability 
test  with  the  balloon  ready  to  be  pressurized. 


2  o=  pg  „ 

r  g 
e  ^o 
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The  apparatus  used  in  the  determination  of  capillary 
critical  radius  is  shown  in  Figure  3 . 14 .  Wick  samples  were 
pressed  between  two  0-rings  and  submerged  in  a  reservoir  of 
working  fluid.  The  sample  covered  a  2.25  inch  hole  which 
was  attached  to  a  2.25  inch  plastic  tube.  The  tube  acted  as 
a  fluid  container  so  that  the  liquid  head  could  be  supported 
by  the  capillary  forces  of  the  wick  sample  structure.  The 
fluid  level  was  then  lowered  until  the  wick  sample  could  no 
longer  support  the  fluid  inside  the  plastic  cylinder  by  sur¬ 
face  tension  forces.  At  that  instant  the  capillary  rise 
height  was  measured.  Critical  radius,  r£,  was  determined 
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using  Equation  3.2.  Figure  3.12  shows  the  capillary  rise 
height  apparatus  broken  down  and  in  the  testing  configuration. 
Figure  3.13  shows  the  measurement  of  wire  mesh  thickness  after 
the  capillary  rise  test.  The  wire  mesh  thickness  data  will  be 
used  in  Chapter  5  to  determine  liquid  recession  depths. 

Figure  3.15  shows  the  apparatus  used  to  manufacture  the 
wire  mesh. wicks.  The  wrapping  mandrel  twenty-five  inches  in 
length  is  mounted  between  spring  loaded  jaws  which  press 
against  the  wire  mesh  as  it  is  wrapped.  Figure  3.16  shows  the 
initial  phase  of  manufacture.  A  piece  of  wire  mesh  is  cut  to 
size  and  a  retainer  rod  is  spot  welded  at  the  edge.  The  wire 
mesh  is  inserted  into  the  mandrel  and  placed  into  the  appara¬ 
tus.  The  wire  mesh  is  then  wrapped  on  the  mandrel  with  the 
spring  loaded  jaws  pressing  against  the  wire  mesh  to  produce 
a  tight  fit  as  shown  in  Figure  3.17.  Figures  3.18  through 
3.23  reiterate  this  sequence  in  detail.  Figure  3.25  shows 
the  final  result  of  the  wick  manufacture.  The  1/4  inch  over¬ 
lap  at  the  edge  of  the  wick  is  necessary  to  insure  strong  spot 
welds.  The  resultant  structure  is  quite  strong  and  incompres¬ 
sible  which  gave  accurate  and  repeatable  readings  in  the  per¬ 
meability  test. 

The  procedure  for  testing  is  given  in  the  appendix. 
Chapter  4  will  deal  with  properties  of  heat  pipe  wicks  mea¬ 
sured  using  the  apparatus  just  described. 


Figure  3.1.  Experimental  Heat  Pipe  Apparatus 


Figure  3.6. 
Inserting  Wire 
Mesh  Wick  into 
Heat  Pipe 


Figure  3.7. 
Testing  Maximum 
Heat  Transfer 
Rate 


Constant  Head 


Figure  3.9.  Apparatus  for  the  Determination  of  Wick  Permeability 
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Figure  3.18.  Initial  Phase  of  Wire  Mesh 
Wick  Manufacture 
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Figure  3.19.  Apparatus  Ready  for  Wrapping 
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Figure  3.21.  The  Wrapping  Is  Initiated  by 

Turning  the  Collet  and  Lowering 
the  Compression  Jaws 
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Additional  Liquid  Flow  Paths 
in  Wirs  Mesh  Keat  Pipes 


CHAPTER  4 

HEAT  PIPE  WICK  PROPERTIES 

The  properties  of  heat  pipe.wicking  materials  are  ex¬ 
tremely  variable.  Wick  properties  such  as  porosity,  critical 
radius,  wetting  angles.,  and  cross  sectional  area  are  not 
identical  for  the  same  type  of  wicks  manufactured  from  the 
same  materials  under  carefully  controlled  uniform  conditions. 

In  order  to  measure  the  extent  of  such  variability  a  total  of 
50  .square  weave  wire  mesh  wicks  were  tested  in  the  following 
categories 

Type  of  Wick  Number  of  Wicks  Tested 

100  Mesh  304  Stainless  Steel  5 

(tight  wrap) 

100  Mesh  304  Stainless  Steel  30 

(moderately  tight  wrap) 

200  Mesh  304  Stainless  Steel  5 

50  Mesh  304  Stainless  Steel  5 

100  Mesh  Copper  5 

Five  properties  were  measured  for  each  wick.  These  properties 
were  permeability  (K) ,  critical  radius  (r£) ,  wick  cross  sec¬ 
tional  area  (A) ,  porosity  (e) ,  and  tortuosity  (b)  . 

In  this  chapter  tabular  values  of  the  experimental  data 
used  to  determine  the  heat  pipe  wick  properties  are  presented. 
The  reason  for  presentation  of  this  data  is  to  show  the  extreme 
variability  of  the  readings  and  to  show  how  one  incorporates 
these  readings  into  estimates  of  the  random  variates  used  in 
design. 


55 


4.1  Permeability 


Permeability  measurements  for  each  of  the  50  wire  mesh 
wicks  tested  are  shown  in  Figures  4.1  through  4.5.  Each  mean 
value  and  standard  deviation  was  generated  using  Equations 
2.55  and  2.56  using  the  number  of  replicated  points  which  is 
designated  as  the  sample  size.  The  individual  means  and  stan¬ 
dard  deviations  result  from  experimental  error  in  observing 
the  data.  The  overall  mean  is  the  mean  of  all  the  individual 
means  and  the  overall  standard  deviation  is  calculated  using 
Equation  2.56  for  all  the  data  points  in  the  data  set.  (Note: 
Each  data  set  is  data  for  a  specific  type  of  wick. )  The  over¬ 
all  mean  and  overall  standard  deviation  are  used  as  the  para¬ 
meters  of  the  random  variate  pair  to  describe  variability  of 
the  manufacturing  process.  According  to  Data  Set  1,  the  per¬ 
meability  variate,  (K,S^) ,  for  the  100  mesh  stainless  steel 

"9  -9 

wicks,  is  (5.243x10  ,  1.096x10  ).  Now  that  the  distribution 

parameter  estimates  are  known,  it  is  helpful  to  assume  a  proba¬ 
bility  distribution  for  the  given  random  variable.  Since 
these  data  are  repeated  readings,  one  generally  chooses  the 
normal  distribution  to  model  variability  of  the  design  vari¬ 
able.  All  of  the  data  were  tested  for  normality  and  did  not 

17 

reject  the  Kolomogorov-Smirnov  Significance  Test  at  the  99% 
significance  level;  thus,  the  data  are  considered  normally 
distributed. 

Factors  contributing  to  the  uncertainty  in  prediction  of 
the  permeability  are: 
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1.  Variability  of  wire  diameter  in  the  mesh  and  between 

mesh  lots 

2.  Weave  manufacture  variability 

3.  Wrapping  manufacture  variability 

a)  Weights  of  wire  mesh  wicks  are  not  identical  in¬ 
dicating  more  wife  mesh  on  certain  wicks. 

b)  Wrap  compression  is  not  completely  consistent  due 
to  change  alignment  of  the  two  screen  layers. 

c)  Tightness  of  wick  fit  in  the  pipe  may  vary  re¬ 
sulting  in  a  variable  seam  channel  flow  path. 

One  may  notice  the  significant  difference  between  the  perme¬ 
ability  in  data  sets  1  and  2  shown  in  Tables  4.1  and  4.2 
Although  each  of  the  wicks  had  identical  mesh  sizes,  the  wicks 
of  data  set  1  were  wrapped  to  fit  very  tightly  in  the  heat 

pipe.  This  tight  fit  compressed  the  screen  structure  which 

/ 

resulted  in  a  low  permeability.  This  demonstrates  the  sensi¬ 
tivity  of  the  permeability  measurement,  to  compression  o£  the 
wick  layers  against  the  pipe  wall. 

The  uncertainty  in  the  prediction  of  heat  flow  capability 
for  a  heat  pipe  with  the  above  permeability  characteristic  is 
directly  proportional  to  the  variability  of  the  permeability. 

4  <,2  Critical  Radius 

Capillary  critical  radius  readings  are  affected  by 

1.  Wetting  angle 

2.  Width  of  wire  mesh  openings 

3.  Compression  of  the  two  or  multiple  layers  of  wire  mesh 
The  data  of  Table  4.6  shows  the  uncertainty  in  the  measurement 
of  capillary  pressure  capability  of  one  and  two  layers  of  wire 
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mesh.  Here  capillary  pressure  capability  is  the  pressure  rise 
measurement  using  the  apparatus  of  Figure  3.14. 

The  data  of  Table  4.6  are  taken  from  data  set  1  and  2* 
since  100  mesh  stainless  steel  screen  was  used  in  both  sets. 

The  equality  of  standard  deviations  between  the  two  layer  and 
one  layer  data  indicates  the  consistency  of  the  capillary  pres¬ 
sure  measuring  device.  All  of  the  single  layer  capillary  pres¬ 
sure  measurements  proved  to  be  statistically  different  from 

17 

the  two  layer  measurements  according  to  a  hypothesis  test. 

It  should  be  pointed  out  that  the  difference  between  the 
one  and  two  layer  capillary  pressure  capabilities  may  vary  de¬ 
pending  on  the  degree  of  compression  and  distortion  of  the  two 
layers.  If  the  two  layers  of  wire  mesh  were  compressed  tightly 
together,  a  rather  large  difference  in  capillary  pressure 
capability  would  be  expected  between  the  compressed  two  layer 
structure  and  the  single  layer  of  wire  mesh  due  to  the  smaller 
capillary  pores  formed  between  the  two  layers.  If  the  two 
layer  structure  was  loose,  very  little  difference  would  be 
expected.  The  wick  structures  mentioned  in  Chapter  3  are 
compressed  sufficiently  to  cause  intermeshing  of  the  multiple 
layers  of  wire  mesh. 

The  difference  in  capillary  pressure  capability  of  one 
and  two  layers  of  the  same  wire  mesh  is  an  important  factor 
in  the  recession  theory  analysis  of  these  types  of  wicks. 

During  the  two  layer  capillary  pressure  tests,  the  liquid 
receded  into  the  two  layer  structure  and  encountered  higher 
capillary  force  resulting  in  the  difference  between  one  and 
*Feference  Appendix  A 
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two  layer  readings.  This  will  be  discussed  in  detail  in 
Chapter  5. 

Table  4.7  shows  small  samples  of  data  from  other  tests 
and  in  every  case  the  one  and  two  layer  capillary  pressure 
capabilities  proved  to  be  different ,  statistically. 

Tables  4.8,  4.9,  4.10,  4.11,  and  4.12  show  the  critical 
radius  values,  r^r ,  for  all  the  data. 

Uncertainty  of  the  widths  of  wire  mesh  openings  (due  to  the 
manufacturing)  are  considered  the  prime  variabilities  in  the 
measurement  of  this  property.  Wetting  angle  was  assumed  con¬ 
stant  throughout  the  structure  since  the  wicks  were  cleaned 
uniformly. 

4.3  Wick  Cross  Sectional  Area 

The  wick  cross  sectional  area  distributions  proved  to  be 
very  consistent  with  low  standard  deviations.  Tables  4.13, 
4.14,  4.15,  4.16,  and  4.17  show  cross  sectional  area  data  for 
all  of  the  area  tests.  The  low  standard  deviations  indicate 
a  consistent  intermeshing  of  the  multiple  layers  of  wire  mesh. 
Uncertainties  in  the  measurement  of  this  property  arise  from 

a)  Differences  in  wire  mesh  lots 

b)  Measurement  variations  for  the  wick  inside  diameter 
(one  must  avoid  additional  compression  of  the  wick 
structure  while  making  the  measurement) 

4  A  Porosity  and  Tortuosity 

Porosity  and  tortuosity  for  all  50  wicks  were  calculated 
from  measured  values  according  to  the  following  definitions 
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Void  Vol  ,  Weight  of  Wick  ...... 

P°roslty  e  -  Total  vol  =  1  -  P7TSTU5 —  (4-1) 


er, 


tortuosity  b  =  — 


K 


(4.2) 


The  porosity  distributions  proved  to  be  almost  identical 
from  wick  to  wick.  Tables  4.18,  4.19,  4.20,  4-. 21,  and  4.22 
show  the  porosity  and  tortuosity  data.  The  tortuosity  data 
were  almost  independent  of  wire  mesh  size  except  for  data  set 
1  which  had  a  very  low  permeability.  One  might  expect  higher 
tortuosity  for  these  wicks  but  looking  at  the  variables  for 
the  calculation  of  b  we  are  able  to  see  the  opposite.  The  rg 
used  in  the  calculation  was  the  value  of  the  two  layer  wick,  (rfr) 
which  was  about  20%  smaller  th^n  that  for  the  single  layer. 

The  permeability  for  this  wick  design  is  higher  than  other 
designs  because  of  the  additional  flow  through  the  seam  flow 
channel  (see  Figure  3.26).  This  tends  to  lower  the  tortuosity, 
b,  somewhat  as  compared  to  other  wire  mesh  structures  since  the 
permeability  seems  to  have  dominated  the  tortuosity  calcula¬ 
tion.  The  standard  deviations  for  e  and  b  were  calculated  from 
the  input  variabilities  into  the  equations  using  the  algebra 

7 

of  moments  method.  The  uncertainty  in  the  calculation  of  b 
was  rather  high  as  a  result  of  the  high  uncertainties  of  the 
input  properties  e,  rQ ,  and  K. 


4.5  Summary  of  Wick  Property  Data 

Table  4.23  summarizes  the  uncertainties  in  predicting 
properties  of  wire  mesh  wicks  of  this  particular  design. 
Many  of  the  properties  may  vary  considerably  depending  on 
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the  type-  ole  manufacture  methods  used.  Since  all  of  the  data 
distributions  are  considered''  normal,  one  can  quote  the 
mean  property  value  plus  or  minus  three  standard  deviations 
and  be  assured  that  99.7%  of  the  property  values  will  lie 
within  these  bounds.  (For  small  sample  sizes,  tolerances 
should  be  placed  on  this  probability.)  Figures  4.1  -  4.5 
show  discrete  probability  plots  of  the  various  design  variables 
which  were  assumed  normally  distributed. 


61 


TABLE  4.1 


PERMEABILITY  DISTRIBUTION  FROM 
DATA  SET  1 


(100  mesh,  stainless 
steel,  304  2  layer) 


Mean  K  (ft2) 

Standard  Deviation  K  (ft2) 

Sample  Size 

3 . 584xl0**9 

1.926xl0~10 

30 

5.244xl0"9 

2.685xl0“10 

5.246xl0~9 

2.874xl0"10 

5.407X10"9 

1.495xl0"10 

6 . 654xl0~9 

42 . 158xl0-10 

s 

/ 

Overall  Mean  5.243x10 

-9 

Overall  Standard  Deviation  1.096x10 
*Note:  These  standard  deviations  are  a  result  of  experimental 
error  in  the  determination  of  each  individual  K  reading.  The 
sample  sizes  are  the  number  of  replications  per  K  reading. 


TABLE  4.2 

PERMEABILITY  DISTRIBUTION  FROM 
DATA  SET  2 


(100  mesh  stainless 
steel,  304  2  layer 
wrap  mandral  .001" 
smaller  diameter  than 
Data  Set  1) 


Mean  K  (ft2) 

2 

Standard  Deviation  K  (ft  ) 

Sample  Size 

5.951  xlO"9 

4.467  xl0~10 

10 

6.080  xlO"9 

3.364  xlO-10 

| 

6.081  xlO**9 

2.297  xlO"10 

6.204  xlO"9 

2.575  XlO"10 

7 

62 


6.234x10 

6.250x10 

6.326x10 

6.412x10 

6.427x10 

6.566x10 


-9 


-9 


-9 


-9 


-9 


6,698x10 

6.850x10 

6.909x10 

6.953x10 

7.184x10 


-9 

-9 

-9 


-9 

-9 

-9 


7.191x10 

7.323x10 

7.463x10 

7.474x10 

7.479x10 

7.482x10 

7.775x10 

8.116x10 


-9 

-9 

-9 

-9 


-9 

-9 

-9 

-9 


1.154x10 

1.214x10 


-9 

-9 


1.246x10 


-9 

-9 


1.337x10 
5 .  338xl0*"9 
3.379x10"9 


1.701x10 


-9 


TABLE  4.2  continued 

.-10 


2.833x10' 

4.356x10 

2.176x10 


-10 


-10 


1.968x10 


-10 

! 

.-10 


3.795x10' 
3.085X10"10 
3.641x10 


3.748x10 

4.862x10 

3.122x10 

3.601x10 

3.192x10 

2.834x10 


-10 

-10 

-10 


-10 

-10 

-10 

-10 


1.381x10 

5.334x10 

2.634x10 

4.464x10 

3.342x10 

4.518x10 

2.877x10 

4.302x10 

2.674x10 


-10 

-10 

-10 


-10 

-10 

-10 

-10 

-10 

-10 


3.730x10 

2.142x10 

5.015x10 


-10 

-10 

-10 


3.439x10 


-10 

.-9 


Overall  Mean  7.193x10  sq.  ft. 
Overall  Standard  Deviation  0.841x10  sq. 


ft. 
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TABLE  4.3 


PERMEABILITY  DISTRIBUTION  FROM 
DATA  SET  3 


Mean  K  (ft2) 

Standard  Deviation 

2.584xl0-9 

1.483X10-10 

2.798xl0"9 

9.819xl0“10 

3 . 074xl0-9 

2.845X10-10 

3.171X10-9 

2.149xl0"10 

3 . 796xl0-9 

6.949X10"10 

(200  mesh  stainless 
steel,  316  3  layer) 


Sample  Size 
10 


w 


-9 


Overall  Mean  3.084x10 
Overall  Standard  Deviation  0.460x10 


sq.  ft. 
-9 


sq.  ft. 


TABLE  4.4 

PERMEABILITY  DISTRIBUTION  FROM 
DATA  SET  4 


Mean  K  (ft2) 

Standard  Deviation 

2.739xl0~8 

-9 

1.339  x10 

2.822X10-8 

-9 

8.494x10 

3 . 209xl0~8 

1 . 625xl0-9 

3 . 230  xlO-8 

9 . 071 xlO-9 

3 . 290  xlO-8 

1 . 625  xlO-9 

(50  mesh  stainless 
steel,  304  2  layer) 


Sample  Size 
10 


W 


Overall  Mean  3.056x10  8  sq.  ft. 

_o 

Overall  Standard  Deviation  0.260x10  sq.  ft. 
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TABLE  4,5 

PERMEABILITY  DISTRIBUTION  FROM 
DATA  SET  5 

(100  mesh  Cu  2  layer) 


Mean  K  (ft^) 

Standard  Deviation  K  (ft; 

Sample 

Size 

6 . 031xl0~9 

4.131xl0"10 

10 

-9 

6.400x10 

3.475xl0~10 

6.42ixl0~9 

4 . 344xl0~10 

6 . 577xl0~9 

3 .lllxlO-10 

6.737xl0-9 

2.586xl0"10 

7 

Overall  Mean  6.433xl0~9  sq.  ft. 

Overall  Standard  Deviation  0.263xl0~9  sq.  ft. 


TABLE  4.6 

CAPILLARY  PRESSURE  CAPABILITY  OF  ONE  LAYER 
AND  TWO  LAYERS  OF  WIRE  MESH 
(100  mesh  stainless  steel  30.4,  35  readings 
taken  from  data  set  #1  and  #2) 


One  layer 
capillary 
pressure 
capability 
(inches  Il^OJ/h^ 

5.000 
4  .625 

5.250 
5.000 

4.625 

4.750 

4.625 

4.750 
4.75Q 

4.500 

5.250 
4.000 

4.750 

5.250 

5.500 
4.675 

4.500 

4.250 
4.125 

4.625 
5.000 

4.750 

5.364 

4.750 

5.250 
5.0C0 

4.875 

4.875 
4.362 
4.375 

4.500 
4.50C 
5.000 

4.750 

4.500 


Two  layer 
capillary 
pressure 
capability,  h ^ 

5.625 

5.250 

5.875 

5.500 

5.  COG 

5.250 

5.500 
5.899- 

5.750 

5.750 
6.250- 

5.625 
5,350- 

6.  P.50- 
6.000 
5.S0C 
5.25C- 

4.625 

5.364 

4.750 

5.625 
5  r  875 

5.750 

5.500 

5.875 

5.750 

5.500 

5.500 

5.625 

5.875 

5.364 

5 . 364 

5.500 

5.500 
5.36& 


Difference 

data 


0  .625 
0.625 
0  .625 
0.500 
0.375 
0.500 
0.875 
!  '  149 
1  .000 
l  .250 
1  .000 
1  .625 
0  .600 
1  .000 
0  .500 
0  .824 
0.750 
0.375 
1  .239 
0.125 
0.625 
1 .125 
0  .335 
0  .750 
0.625 
0 .750' 
0.625' 
0.625- 
1  .262- 
1  .500* 
0.864- 
0.864' 
0.500' 
0.7501 
0.864 


Mean  4.761 


5.552 


0.790 


Stan-  0,349  0.350  0.332 

dard 

deviation 

(all  distributions 

Normal) 
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TABLE  4.7' 

CAPILLARY  PRESSURE  CAPABILITY  OP  ONE  LAYER 
AND  TWO  LAYERS  OP  WIRE  MESH 
(data  taken  from  data  sets  3,4,5) 


One  Layer 
capillary 
pressure 
capability 
,nches  H2O)  ,h^ 

Two  layer 
capillary 
pressure 
capability 

h2 

Difference 

data 

h2  -  hl 

(Data  set  3, 

200  mesh 

stainless 

steel, 

3  layer) 

(note:  for 
data  set  3, 
this  is  3 
layer  capability) 

8.000 

9.750 

1.750 

7.000 

9.500 

2.500 

7.500 

9.000 

1.500 

9.250 

10.250 

1.000 

7.350 

8.750 

1.400 

7 . 819 (Overall  Mean)  9.450 


_ (Overall  Stan- 
dard  Deviation) 


0.597 


1.630 

0.556 


(Data  set  4, 

50  mesh 

stainless 

steel, 

2  layer) 

2.364 

2.875 

0.510 

2.364 

2.750 

0.385 

2.000 

2.625 

0.625 

1.875 

2.625 

0.750 
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TABLE  4.7 

(continued) 

2.250  2.625 

0.375 

2.171 (Overall  Mean)  2.700 

0  i  529 

n  (Overall  Stan-  n  9 

dard  Deviation 

0.160 

(Data  set  5, 

100  niesh 

Cm,-  2  layer) 

5.250 

5.875 

0.625 

5.250 

5.875 

0.625 

5.164 

5.750 

0.585 

5  i  250 

5.875 

0.625 

5.250 

5.875 

0.625 

5. 233 (Overall  Mean) 

5,850 

0.617 

n  mn  (Overall  Stan- 
*  dard  Devia- 

0,056 

0.018 

tion) 


TABLE  4.8 


CRITICAL  RADIUS  DISTRIBUTION 


FROM  DATA  SET  1 


Mean  r  _  inches) 
rr 


4  055x10 
4.247x10' 


4.337x10 


4.809x10 

4.906xi0~3 


(100  mesh  stainless 
steel,  two  layer) 


Standard  Deviation (inches) 
_ •  fr _ 

l.OOOxlO-4 


Sample  Size 


Overall  Mean  4.470x3  0  inches 
Overall  Standard  Deviation  0.369xl0-:1  inches 


TABLE  4.9 


CRITICAL  RADIUS  DISTRIBUTION 


FROM  DATA  SET  2 


Mean  r^r( inches) 


(100  mesh  stainless 
steel,  two  layer) 

Standard  Deviation  rfj,{ inches)  Sample  Size 


3.639x10 

3.645x10 


1.000x10 


3.786x10 

3.858xl0~3 

3. 870xl0_3 

3.877X10"3 

3.878xl0“3 

3.887xl0~3 
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TABLE  4 1,9  (continued) 


3.958xl0"3 
3.961xl0"3 
3. 9 69x10“ 3 
3.969xl0“3 
5. 051xl0~3 
5.051xl0"3 
5.053xl0"3 
4. 139x10" 3 


5.140x10 


-3 


4.140x10 

4.140x10 


-3 


r>~3 


-3 


4.140x10 

4.140x10 

4.141x10 

4.143x10 

4.245x10 

4.257x10 

4.257x10 

4.265x10 

4.337x10 

4.337X10 


4, 539x10"  3 


-3 


-3 


-3 


-3 


-3 


-3 


-3 


-3 


-3 


V 


-3 


Overall  Mean  4.069x10  inches 

-3 

Overall  Standard  Deviation  0.214x10  inches 

Ncte:  All  r^r  reac.ingc  are  those  of  the  multiple  layers  of 
screen. 


70 


TABLE  4.10 
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TABLE  4.12 


Mean  r^r 

3.888x10 

3,888x10' 

3.888x10 

3.888x10 

3.971x10 


(inches) 


-3 

-3 

-3 

-3 

-3 


CRITICAL  RADIUS  DISTRIBUTION 
FROM  DATA  SET  5 


(100  mesh  copper, 
two  layer) 


Standard  Deviation  r^r (inches) 


Sample  Size 


1.000x10 


-4 


V 


Overall  Mean  3.904x10 


-3 


Overall  Standard  Deviation  1.000x10 


-4 


TABLE  4.13 


WICK  CROSS  SECTIONAL  AREA  DISTRIBUTION 

(Data  from  data  set  1,  100 
mess  stainless  steel,  tv7o  layer) 


Mean  A (sq.ft.) 

Standard  Deviation  A 

Sample  Size 

2.285xl0~4 

5 . 040xl0~5 

16 

2.354xl0“4 

6.289xl0-5 

2.375xl0-4 

5. 250xl0-5 

2.445X10"4 

4.790xl0~5 

2.503X10"4 

1.104xl0**5 

7 

-4 

Overall  Mean  2.392x10  sq. 

ft. . 

-5 

Overall  Standard  Deviation  0.848x10  sq.  ft. 


TABLE  4.14 

WICK  CROSS  SECTIONAL  AREA  DISTRIBUTION 

(Data  from  data  set  2,  100 
mesh  stainless  steel,  two  layer) 

Mean  A(sq.ft.)  Standard  Deviation  A (sq.ft)  Sample  Size 


2.236xl0~4 

6.790X10”6 

-4 

2.261x10 

7.510X10-6 

2. 270xl0-4 

4.900X10-6 

2. 294xl0-4 

2.294X10"6 

2. 324xl0-4 

5.370X10-6 

2.334xl0*~4 

6.270xl0-6 

2. 336xl0-4 

1.124X10"3 

2.337xl0~4 

6.280X10"6 

2.353X10"4 

7 . 900xl0~6 

2.359X10'4 

7.890xl0-6 

2.363xl0”4 

7.780xl0~6 

2.373X10"4 

6.270xl0~6 

2.382X10"4 

1.236xl0"5 

2.387X10*’4 

1 . OOlxlO-5 

2. 387xl0~4 

6.00010“6 

2.392xl0~4 

9.78010"6 

2. 392xl0~4 

6.880xl0~6 

2. 401xl0-4 

8 . 750xl0~6 

2. 402xl0-4 

9.200xl0”6 

2.416xl0~4 

1.022xl0~5 

2.416X10-4 

1.134xl0~5 

2.421X10"4 

1.300x10  5 

2. 426xl0~4 

8.500X10-6 

2.426xl0"4 

1.169X10"5 

2. 441xl0-4 

1.117X10"5 

2. 469xl0-4 

1.674X10-5 

2. 475xl0~4 

1.479X10-5 

2. 555xl0-4 

1.190X10”5 

2.558xl0"4 

1.189X10”5 

2.584xl0-4 

6 . 540x10*  3 

Overall  Mean  2.393x10  3  sq.  ft. 

Overall  Standard  Deviation  0.816x10  5  sq.  ft. 
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TABLE  4.15 


WICK  CROSS  SECTIONAL  AREA  DISTRIBUTION 

(Data  from  data  set  3,  200  mesh, 
stainless  steel,  3  layer) 


Mean  A(sq.  ft.) 

Standard  Deviation  A 

Sample  Size 

1.579xl0~4 

5.390x10  6 

1 

6 

1.614xl0~4 

1.294X10-5 

1.677xl0"4 

6.440X10-6 

,  „-4 

,  -6 

1.682x10 

7.860XX0 

r 

7 

-4 

-6 

Y 

1.688x10 

8.570X10 

-4 

Overall  Mfan  1.648x10  sq.  ft. 

Overall  Standard  Deviation  0.487x10  5  sq.  ft. 


TABLE  4.16 

WICK  CROSS  SECTIONAL  AREA  DISTRIBUTION 


(Data  from  data  set  4,  50  mesh 
stainless  steel,  2  layer) 


Mean  S(sq.  ft.) 

Standard  Deviation  A 

Sample  Size 

5.203xl0“5 

2.091xl0~6 

16 

5. 295xl0-5 

l.lOOxlO-5 

5.314xl0-5 

1.890X10-5 

5.406X10-5 

1.860X10-6 

5.438xl0-5 

l.OOlxlO-5 

7 

-4 

Overall  Mean  5.331x10  sq.  ft. 

-5 

Overall  Standard  Deviation  0.936x10  sq.  ft. 
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TABLE  4.17 


WICK  CROSS  SECTIONAL  AREA  DISTRIBUTION 


(Data  from  data  set  5,  100  mesh 
Cu  2  layer) 


Mean  A(sq,  ft.) 

Standard  Deviation  A 

Sample  Size 

2,759xl0”4 

1.719xl0**5 

16 

2.836X10"4 

1. 474xlo""5 

2.893X10"4 

1. 037xl0”5 

2.939xl0"4 

1. 053xl0”5 

3.092xlQ-4 

9 . 650xl0~6 

7 

-4 

Overall  Mean  2.903x10  sq.  ft. 

-4 

Overall  Standard  Deviation  0.124x10  sq.  ft. 
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TABLE  4.18 


POROSITY  AND  TORTUOSITY  DISTRIBUTIONS 


(Data  from  data  set  1,  100 
mesh  stainless  steel,  two 
layer) 


Mean 

Standard 

Mean 

Standard 

Porosity 

Deviation 

Tortuosity 

Devistion 

(dimensionless) 

Porosity 

Dimensionless)  Tortuosity 

0.582 

0.0^24 

18.553 

3.244 

0.577 

0.113 

18.944 

2.827 

0.591 

0.0898 

14.113 

2.376 

0.614 

0.0180 

19.561 

1.403 

0.602 

0.0782 

14.322 

2.016 

Overall 

Mean 

0.593 

16.098 

Overall 

Standard 

Deviation 

0.0359 

2.727 

TABLE  4.19 

POROSITY  AND  TORTUOSITY  DISTRIBUTIONS 


Mean 

Porosity 

Standard 

Deviation 

Porosity 

0.596 

0.0249 

0.636 

0.0107 

0.611 

0.0190 

0.583 

0.0210 

0.619 

0.0186 

0.606 

0.0273 

(Data  from  data  set  2,  100 
mesh  stainless  steel,  two 
layer) 


Mean 

Tortuosity 

Standard 

Deviation 

Tortuosity 

10.257 

0.975 

8.472 

0.557 

8.229 

0.687 

8.794 

0.740 

9.836 

0.752 

8.927 

0.888 
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TABU  4.19 

0.607 

0.0189 

0.597 

0.0166 

0.594 

0.0148 

0.604 

0.0151 

0.578 

0.0168 

0.569 

0.0132 

0.582 

0.0233 

0.584 

0.0186 

0.585 

0.0168 

0.577 

0.0137 

0.574 

0.0154 

0.591 

0.0206 

0.580 

0.0226 

0.587 

0.0203 

0.573 

0.0158 

0.578 

0.0130 

0.563 

0.0114 

0.585 

0.0184 

0.584 

0.0122 

0.570 

0.0157 

0.552 

0.0143 

0.581 

0.0127 

0.578 

0.0116 

0.589 

0.0179 

Overall 

Mean 

0.588 

Overall 

Standard 

Deviation 

0.0170 

(continued) 


8.157 

0.663 

10.688 

0.772 

9.721 

0.652 

8.757 

0.661 

10.323 

0.690 

9.984 

0.610 

10.519 

0.984 

9.106 

0.742 

11.644 

0.833 

9.136 

0.645 

9.618 

0.710 

9.300 

0.707 

9.225 

0.795 

10.312 

0.866 

10.303 

0.708 

7.924 

0.525 

8.116 

0.473 

10.375 

0.863 

9.885 

0.685 

1.704 

0.551 

10.841 

0.670 

8.314 

0.572 

9.219 

0.511 

9.482 

0.670 

9.473 

0.937 
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TABLE  4.20 


POROSITY  AND  TORTUOSITY  DISTRIBUTIONS 


(Data  set  3,  200  mesh 
stainless  steel,  three 
layer) 


Mean 

Porosity 

Standard 

Deviation 

Porosity 

Mean 

Tortuosity 

Standard 

Deviation 

Tortuosity 

0.631 

0.0142 

10.074 

0.600 

0.622  . 

0.00588 

9.221 

0.463 

0.641 

0.00553 

10.186 

0.562 

0.645 

0.00844 

11.941 

0.361 

0.630 

0.00950 

10.184 

0.287 

Overall 
Mean  0 , 

.634 

10.321 

Overall  0 

Standard 

Deviation 

.00933 

0.992 

TABLE  4.21 

POROSITY  AND 

TORTUOSITY  DISTRIBUTIONS 

(Data  set  4 ,  50  mesh 
stainless  steel,  two 
layer) 

Mean  P 
Porosity 

Standard 

Deviation 

Porosity 

Mean 

Tortuosity  ' 

Standard 

Deviation 

Tortuosity 

0.617 

0.0312 

7.617 

1.035 

0.614 

C.0145 

9.468 

0.985 

0.656 

0.0183 

7.957 

1.598 

0.640 

0.0149 

7.428 

2.701 

0.653 

0.0171 

9.696 

1.025 

Overall 
Mean  0 

.636 

8.433 

Overall  0 

Standard 

Deviation 

.0197 

1.068 
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TABUS  4.22 

POROSITY  AND  TORTUOSITY  DISTRIBUTIONS 


(Data  set  5,  100  mesh,  Cu 
2  layer) 


Mean  Standard 

Porosity  Deviation 

Porosity 

Mean 

Tortuosity 

Standard 

Deviation 

Tortuosity 

0.648  0.0122 

10.105 

0.676 

0.600  .  0.0256 

10.450 

1.000 

0.634  0.0142 

10.358 

0.911 

0.626  0.0202 

10.277 

0.837 

0.623  0.0147 

10.373 

0.757 

Overall  0.626 

Mean 

10.313 

Overall  0.0175 

Standard 

Deviation 

0.313 

to  ft  ff) 
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CHAPTER  5 

RECESSION  IN  WIRE  MESH  WICKS 

In  Chapter  4  the  basic  determination  of  wire  mesh  wick 
properties  was  described.  In  this  chapter  the  application 
of  the  recession  theory  of  Chapter  i  to  the  wire  mesh  wicks 
will  be  presented.  The  wire  mesh  wicks  t~  be  analysed  will 
be  the  two  layer  square  weave  type  wici,J  that  were  used  in 

the  experimental  measurements. 

7 

Many  descriptions  of  wire  mesh  wick  characteristics  have 
modeled  the  surface  tension  phenomenon  as  shown  in  Figure  5.1. 
The  screen  layers  are  spaced  one  screen  opening  apart,  and 
the  radius  is  equal  to  one  half  the  screen  opening  size.  Li¬ 
quid  recession  is  shown  in  Figure  5.2(a).  This  desaturation 
process  is  unstable  since  the  capillary  forces  have  not  in¬ 
creased  due  to  the  recession  (Figures  5.2(a)  and  (b) ) .  Using 
this  physical  model,  the  wick  would  dry  out  at  relatively  low 
heat  transfer  rate,  yet  experimentation  shows  this  not  to  be 
the  case.  Many  times  fluid  recession  occurs  down  to  the  se¬ 
cond  layer  and  remains  stable,  and  no  dry  out  is  observed  at 
higher  than  expected  heat  transfer  rates.  Thus,  a  more  ac¬ 
curate  description  of  the  wire  mesh  capillary  structure  is 
apparently  required  to  help  explain  these  observations. 

Figure  5.3  shows  a  sketch  of  plain  square  weave  wire  cloth. 
Figure  5.4  shows  a  sketch  of  two  layers  of  wire  cloth  compressed 
together  as  occurs  in  the  construction  of  a  heat  pipe  wick. 
Generally,  when  wicks  are  wrapped,  the  two  layers  do  not  align 
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perfectly  as  suggested  in  Figure  5.1.  Normally  the  wrapping 
results  in  the  structure  shown  in  Figure  5.5.  The  high  points 
of  the  bottom  layer  of  wire  mesh  tend  to  fit  into  openings  of 
the  top  layer,  resulting  in  a  structure  with  a  thickness  less 
than  the  sums  of  the  thicknesses  of  the  individual  layers. 
Figures  5.6,  5.7,  and  5.8  show  a  cross  section  view,  A-A,  of 
the  compressed  doube  layer  of  wire  mesh  previously  shown  in 
Figure  5.4. 


Figure  5.4.  Compressed  Double 
Layer  of  Wire  Mesh 


Figure  5.5.  Proposed  Wire  Mesh  Capillary  Model 


Figure  5.7.  Two  Layers  of  Screen 'Wick  with  Liquid  Recessed  to  the  First  Layer 


■aye 


In  Figure  5.5,  the  simplified  drawing  of  the  intermeshing 

of  the  two  layers  shows  that  a  smaller  capillary  radius  exists 

at  the  interface  of  the  two  layers.  Figures  5.6,  5.7,  and  5.8 

are  detailed  drawings  of  Figure  5.5.  The  mean  thickness  of  a 

double  layer  of  100  mesh  stainless  steel  was  .0160  inches  with 

a  standard  deviation  of  .00107  inches  based  on  72  readings. 

In  Figure  5.6,  the  fully  saturated  wick  is  shown  with  a  large 

capillary  radius  at  the  t  level.  In  Figure  5.7,  liquid  is 

w 

shown  recessed  to  the  first  layer  of  wire  mesh.  This  recessed 
level  will  be  termed  the  first  layer  capillary  pressure  capa¬ 
bility  of  the  wick  and,  when  the  liquid-vapor  interface  re¬ 
cedes  to  that  level,  the  fluid  experiences  a  capillary  force 
characteristic  of  the  first  layer.  The  distance  from  the 
wall  of  the  first  layer  is  called  the  recessed  depth,  t^,  and 
is  graphically  determined  from  Figure  5.7.  The  graphical  de¬ 
termination  of  t^  is  performed  by  measuring  the  distance  from 
the  pipe  wall  to  the  hypothesised  liquid  level  on  a  large 
cross  sectional  drawing  of  the  wick.  For  the  first  layer, 
t^  is  three  fourths  the  thickness  of  the  wick,  t  ,  with  the 
same  standard  deviation.  Therefore,  for  a  100  mesh  stainless 
steel  wick,  the  liquid  experiences  a  .capillary  force  equiva¬ 
lent  to  the  first  layer  capability  at  a  distance  .00120  inches 
from  the  wall  with  a  standard  deviation  of  .00107  inches.  The 
first  layer  capillary  pressure  capability  was  determined  ex¬ 
perimentally  and  is  tabulated  in  Tables  4.6  and  4.7  of  Chap¬ 
ter  4. 
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If  additional  capillary  force  is  required  of  the  wick,  the 
liquid  level  will  receed  to  the  configuration  shown  in  Figure 
5.8.  The  liquid- vapor  interface  has  now  encountered  a  smaller 
capillary  radius,  r2.  This  configuration  is  termed  the  two 
layer  capillary  pressure  capability,  the  data  for  this  liquid 
level  tabulated  in  Tables  4.6  and  4.7.  The  recessed  depth  of 
this  level  of  recession,  t2,  is  again  determined  graphically 
and  averages  59%  of  the  wick  thickness  with  the  standard  der 
viation  of  the  wick  thickness.  For  the  100  mesh  stainless 
steel  wick  the  liquid  level  will  encounter  a  capillary  pres¬ 
sure  equal  to  the  two  layer  capillary' pressure  capability 
.0094  inches  from  the  wall  of  the  heat  pipe  with  a  standard 
deviation  of  00107  inches.  If  additional  capillary  force 
is  required  of  the  wick,  no  smaller  rg  will  be  encountered 
during  the  recession  and  the  wick  will  dry  out. 

To  add  to  the  information  presented  in  Figures  4.7  and 
4.8  we  will  refer  to  Figures  5.9  and  5,10  which  are  photo¬ 
graphs  of  the  wick  structures  tested  in  this  work.  The  photo¬ 
graph  of  the  100  mesh  two  layer  wick  (12. 5x)  of  Figure  5.9 
shows  the  intermeshing  of  the  two  layers  of  wire  mesh  consis¬ 
tent  with  Figures  5.6  and  d.7.  Figure  5,12  shows  a  top  view 
of  the  intermeshing  and  illustrates  how  the  high  points  of 
the  lower  layer  tend  to  fill  the  openings  of  the  top  layer 
of  wire  mesh.  Figure  5.11  is  a  sequential  recession  of  fluid 
in  a  100  mesh  stainless  steel  two  layer  wick.  High  intensity 
lamps  were  used  to  observe  the  liquid  vapor  interface  as  it 
recedes.  The  liquid  receded  to  a  location  between  the  two 
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200-mesh  stainless  steel 
3-layer  wick  2X 


200-mesh  stainless  steel 
3-layer  wick  12. 5X 


100-mesh  stainless  steel 
2-layer  wick  2X 


100-mesh  stainless  steel 
2-layer  wick  12. 5X 


Figure  5.9. 


Wire  mesh  structures  for  200-  and  100- 
mesh  stainless  steel  wicks 
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50-mesh  stainless  steel 
2-layer  wick  2X 


50-mesh  stainless  steel 
2-layer  wick  12. 5X 


100-mesh  copper 
2-layer  wick  2X 


Figure  5.10.  Wiremesh  structures  for  50-mesh  stainless 
steel  and  100-mesh  copper  wicks 


Fully  saturated 
wick  at  low  watt¬ 
age 
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Wattage  increased, 
recession  begins 


Recession 

Increases 


Liquid  has  re¬ 
recessed  to 
minimum  capil¬ 
lary  radius 


Burnout  begins 
randomly 


Wick  function 
terminated  at 
complete  burnout 


Figure  5.11.  Sequential  observation  of  liquid  reces¬ 
sion  and  final  burnout  for  100-mesh 
stainless  steel  2-layer  wick 


Figure  5.12. 
Comparison  of 
Randon  Inter¬ 
meshing  of  Two 
Layers  of  100- 
Mesn  Stainless 
Steel  Screen 
Compressed  To- 


layers  where  maximum  capillary  forces  were  encountered.  When 
txie  heat  transfer  rate  was  increased  further,  the  liquid-vapor 
interface  receded  to  dry  out.  Maximum  capillary  forces  were 
encountered  at  the  interface  of  the  two  screen  layers  which 
is  illustrated  in  Figure  5.8.  Tables  4 . 6  and  4.7  of  Chapter 
4  indicate  a  significant  capillary  force  difference  between 
one  and  two  layers  of  wire  mesh. 

Figure  5.14  shows  the  three  possible  configurations  of 
the  liquid-vapor  interface  that  will  be  observed  in  this 
analysis.  Figure  5.14  (a)  shows  the  fully  saturated  wick 
with  large  capillary  radii  at  the  t  level.  The  capillary 
pressure  capability  of  this  level  is  small  and  will  be  ne¬ 
glected.  This  configuration  of  the  liquid  vapor  interface 
is  assumed  for  the  condenser  and  adiabatic  sections  of  the 
heat  pipe  considered  here. 

The  partially  recessed  configuration  of  Figure  5.14  (b) , 
where  the  wick  is  saturated  at  a  level  where  tn  =  ,75t  is 
assumed  to  exist  from  x  =  Lfr  to  x  =  Lg.  The  configuration 
of  Figure  5.14  where  the  liquid  is  fully  recessed  to  the  level 

where  t«  =  .  59t  is  assumed  to  exist  from  x  =  0  to  x  =  L-  . 

2  w  fr 

The  reason  for  this  detailed  analysis  is  to  determine  the 

actual  saturation  distribution  as  shown  in  Figure  5.13  (b) , 

of  the  wick  so  that  L' (A)r  may  be  calculated.  Now  that  the 

values  of  r  at  various  levels  are  specified,  we  can  calculate 
e 

jj'  (A)  using  Equation  2.52.  Applying  Equation  2.14  to  the  evap¬ 
orator  section  of  our  heat  pipe  we  obtain  the  expression  for 
liquid  pressure  loss  and 
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Ap* =  4  J 

X, 

m(x)  ,  ym 

A [ST  PKA 

(-a  ♦  $) 
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+  La  + 


(5.2) 


L' (A)  is  the  expression  that  accounts  for  the  reduction  of  flow 
area  due  to  recession.  Equation  5.2  can  now  be  used  in  con¬ 
junction  with  Equations  2.45,  2.53,  and  5.1  in  the  iterative 
solution  discussed  in  Chapter  2. 

This  concludes  discussion  of  the  deterministic  recession 
model.  Chapter  6  will  outline  computation  techniques  for  the 
solution  of  Equation  2.52  using  probabilistic  techniques.  The 
analysis  presented  in  this  chapter  is  considered  valid  for  any 
size  square  weave  wire  mesh  wrapped  according  to  the  procedure 
shown  in  Chapter  3.  Although  this  particular  analysis  was  done 
with  two  layer  wicks,  the  effects  of  multiple  layers  (greater 
than  two)  could  be  incorporated  into  this  analysis.  For  in¬ 
stance,  this  analysis  was  used  successfully,  on  a  three  layer 
200  mesh  wick  assuming  no  recession  into  the  third  layer. 
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CHAPTER  6 


COMPUTATIONAL  METHODOLOGY 


6.1  Methodology 

In  this  chapter  the  foundation  is  developed  for  the  solu¬ 
tion  of  Equation  2.57,  the  probabilistic  model  of  heat  pipe 
operation.  Simple  deterministic  equations  are  used  as  ex¬ 
amples  of  the  formulation  of  a  probabilistic  design  model. 

The  probabilistic  model  uses  the  basic  equation  structure 
of  the  deterministic  model.  Probabilistic  methodology  is 
applied  to  account  for  the  variability  in  the  deterministic 
model. 

To  illustrate  the  solution  to  the  probabilistic  design 
equation,  simple  mathematical  models  will  be  used. 

Many  physical  phenomena  may  be  described  by  mathematical 
models.  A  simple  example  is  shown  below. 

Q3  =  Qi  +  Q2  (6.1) 

Q  =  Flow  Rate 

Assuming  that  this  deterministic  model  describes  the  physical 
phenomena  adequately,  we  may  form  a  probabilistic  model  by 
transforming  the  deterministic  variables  into  random  variables 
which  are  governed  by  some  statistical  distribution 


^Q3Pl'Q3P2' *  *  *  Q3Pn^  ~  ^1^>1,^1^>2 f  *  *  *^lPn^ 

+  (Q2P1'Q2P2' ' * ,Q2Pn* 


Q-j_fQ2#Q3  =  Random  variables 

QP  =  Random  variable  distribution  parameters 

Ql,Q2  =  Independent 

For  our  particular  analysis,  virtually  all  distributions  of 
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physical  properties  are  two  parameter.  Therefore*  the  - 
tribution  is  determined  solely  by  the  mean  and  standard  devi¬ 
ation.  We  may  now  transform  our  general  probabilistic  model 
to  that  of  simple  two-parameter  distributions  and  obtain 

<63'V  =  <VSQl>  +  (q2'sq2>  (6-3) 

where  (Q,Sq)  is  a  random  variate  pair. 

The  reason  for  changing  the  deterministic  variable  to  a 
random  variable  is  that  most  physical  variables 
may  not  be  known  precisely.  If  there  is 

uncertainty  in  describing  these  input  variables,  there  will 
be  uncertainty  in  the  functional  result.  Describing  the  un¬ 
certainty  of  the  functional  result  may  be  the  best  way  to 
describe  the  natural  phenomenon. 

There  are  three  techniques  for  finding  functional  vari¬ 
ability  and  they  are  discussed  as  follows: 
a)  Simulation:  We  will  use  the  flow  rate  example  to  illus¬ 
trate  the  procedure.  We  are  given  the  following 

,  Q.  =  Plow  Rate  ~  ZZ*(Q, ,S  ) 

°3  =  Q1  +  Q?  1  <6  *> 

Q2  =  Flow  Rate  ~  zz2 (^2 (^2 ,SQ2) 

This  solution  technique,  uses  a  random  number  generator  to 
supply  a  ZZ1  distributed  flow  rate,  and  a  ZZ2  distributed 
flow  rate,  Q2.  These  individual  random  numbers  are  inserted 
into  the  above  equation  and  the  random  result,  Q3,  is  recorded. 
This  procedure  is  repeated  many  times,  and  the  Q3  parameters 
and  distribution  can  be  determined  using  Equations  2.55  and 
2,56  and  a  significance  test. 

*  ZZi  designates  a  particular  two  parameter  distribution  such 
as  the  normal  distribution. 
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Figure  6.1.  Simulation  of  Functional  Variability 


The  simulation  algorithm  is  as  follows 


Do  down  to. and  including  5,  n  times 
Call  ZZ^  Distributed  Random  Q1 
Call  ZZ2  Distributed  Random  Q2 
Random  Q3  =  Random  Q-^  +  Random  Q2 
Store  Random  Q3 
5  Continue 

Call  Subroutine  Parameter  Estimation 

Call  Significance  Test 

End 

From  this  technique  we  can  determine  the  functional  mean, 
variability  (standard  deviation)  and  distribution, 
b)  Partial  Derivative  Method:  This  method  uses  a  Taylor  ap¬ 
proximation  of  functional  variability  and  is  derived  in  most 
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elementary  statistical  texts.  To  find  functional  variability 
using  this  method  we  proceed  as  follows: 


Given 


P  =  f x2...  xn) 

E  (F)  *  F  =  f  (E  (xx)  ,E  (x2) . .  .E  (xn) ) 
All  X  Independent 


F  «  f (xlfx2. ..xn) 


(6.5) 

(6.6) 

(6.7) 


(the  mean  function  result  is  the  computation  of  all  the 
variables  means  in  the  given  equation) 

The  functional  variability  (standard  deviation)  is  derived 

from  the  Taylor  approximation  as, 

and  for  our  example 


Q3  -  Qx  -  Q2 

SQ3  "  <SQl2  +  SQ22)1/2 


Note  that  this  method  yields  a  good  approximation  of  the  func¬ 
tional  mean  and  standard  deviation  but  indicates  nothing  of 
the  functional  distribution. 

c)  Algebra  of  moments:  This  method  uses  the  algebra  of  moments 
or  expectations  to  determine  functional  variability  exactly. 
Unfortunately,  algebra  of  expectations  becomes  difficult  when 
functions  contain  many  random  variables  and  powers  of  random 
variables.  Because  of  this  difficulty  this  method  is  used  only 
on  simple  functions  and  these  results  are  shown  in  Table  6.1. 
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TABLE  6.1 


MEANS  AND  STANDARD  DEVIATIONS  FOR 
SIMPLE  FUNCTIONS  USING  ALGEBRA  OF  EXPECTATIONS 


Function 

Mean 

Standard  Deviation 

a=b+c 

a=b+c 

sa='sb2+Sc2>V2 

a=b-c 

a=b-c 

Sa=(Sb2+Sc2)1/2 

a=bc 

a=bc 

-2  2-2  2 
sa=(b  sc  +C  Sb 

+s/s  2: 
b  c 

a=b/c 

a=b/c 

rc2sb2+s2s 2 

^  _  b  c 

1/2 

aL  c4 

,2 

a=b 

“  u2 

a=b 

S^(45\2+2Sb4 

}  l/2 

b,c  independent 

The  algebra  of  expections  gives  the  functional  mean  and 
the  standard  deviation  but  does  not  indicate  the  func¬ 
tional  distribution.  For  example,  the  algebra  of  expec¬ 
tations  gives  the  following  exact  result  which  is  in¬ 
cidentally  the  same  as  the  partial  derivative  result. 


°3  ”  Q1  +  Q2 


/  2  2  N  1/2 

(  Q1  °2  ) 


Examples  of  Functional  Uncertainty 


The  three  sample  functions  we  will  investigate  are* 

1)  Q3  =  Qx  +  Q2 

2}  A  =  B  -  C 
3)  V  =  Q/A 

The  first  function,  +  Q2,  appears  in  Table  6.2.  The  stan¬ 
dard  deviation  of  the  function  was  calculated  using  simulation, 
partial  derivative  and  exact  techniques  at  various  standard 
deviations  and  sample  sizes.  Looking  at  the  range  of  input 
variable  standard  deviations ,  one  can  conclude  that  this  func¬ 
tion  is  well  behaved.  The  functional  standard  deviations  are 
actually  smaller,  percentage  wise,  than  the  input  variable 
standard  deviations.  Even  at  large  input  standard  deviations 
this  holds  true  and  all  techniques  give  the  same  resultant 
standard  deviation. 

The  second  function,  A  =  B  -  C,  is  not  so  well  behaved. 
Table  6.3  shows  what  can  happen  if  the  input  variables,  B  &  C, 
are  nearly  the  same  in  magnitude  with  high  standard  deviations. 
This  is  the  classic  numerical  problem  of  subtracting  two  num¬ 
bers  of  almost  equal  magnitude.  Notice  that  at  moderate  input 
standard  deviations  (say  9%  of  the  input  variable  mean)  the 
functional  standard  deviation  grows  to  225%  of  the  functional 
result.  The  simulation  mean  becomes  unstable  at  these  high 
standard  deviations  as  one  might  expect. 

*  The  reason  for  choosing  these  simple  functions  is  that  they 
appear  in  many  engineering  applications  and  also  form  the 
basic  structure  of  Equation  2.57. 
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*(XX  =  Q3  =  Qx  +  Q2,  QrN(10,  SQ  )  and  Q2  ~N(5,SQ  ) 

1  2 

**In  this  case  the  partial  derivative  and  algebra  of  moments  method  are  identical. 


TABLE  6.3 


STANDARD  DEVIATION  GROWTH  OF  THE  FUNCTION 
YY*  USING  SIMULATION,  PARTIAL  DERIVATIVE 
AND  ALGEBRA  OF  MOMENTS  METHOD 


Random 

Variable 

Standard 

Deviation 

(in  %  of  ■ 

mean) 

Mean 

Simulation 

Stdv. 

.1 

50.0098 

1.3278 

1 

48.6319 

15.3096 

2 

51.6134 

28.9453 

3 

48.6156 

43.4557 

4 

47.0834 

53.3973 

5 

48.5659 

70.2449 

6 

53.0261 

83.8806 

7 

44.2171 

99.3238 

8 

51.4417 

107.2407 

9 

53.8297 

131.5490 

*(YY  =  A  =  B 

-  C,  B  - 

.  N(iooo,sb: 

.  Partial 

Derivative  & 
Algebra  of 

Sample  Moments  Method** 


Size 

Mean 

Stdv. 

300 

50.00 

1.3793 

300 

50.00 

15.1724 

300 

50.00 

28.9655 

300 

50.00 

42.7586 

300 

50.00 

56.5571 

300 

50 . 00 

70.3448 

300 

50.00 

70.9831 

300 

50.00 

97.9310 

300 

50.00 

111.7231 

300 

50. 00 

125.5173 

,  C  ~  N(950,SC)) 


**  In  this  case  the  partial  derivative  and  algebra  of  moments 
method  ave  identical. 
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case  the  partial  derivative  and  algebra  of  moments  method  are 


The  third  function,  .V  =  Q/A,  presents  a  problem  when  de¬ 
nominator  random  variables  approach  zero.  Table  6.4  shows  the 
standard  deviation  growth  of  this  function.  Notice  that  in 
this  case  the  partial  derivative  and  exact  method  diverge  from 
the  simulation  method.  The  simulation  is  affected  considerably 
by  random  choices  of  A  near  zero  while  the  other  techniques 
are  not.  This  divergence  may  be  a  good  indicator  of  numeri¬ 
cal  problems  in  this  function. 

The  purpose  of  these  examples  is  to  point  out  the  numeri¬ 
cal  problems  associated  with  the  computation  of  thei  functional 
standard  deviation.  Caution  must  be  observed  in  this  calcula¬ 
tion  since  the  physical  phenomenon  variability  may  differ, 
significantly,  from  the  calculation. 

6.2  Application  to  Heat  Pipe  Design 

In  the  next  several  pages  the  three  techniques  are  pre¬ 
sented  for  finding  functional  variability  as  applied  to 
Equation  2.57. 


a)  Simulation 

The  simulation  technique  is  the  most  useful  because  the 
distribution  of  the  heat  transfer  rate  can  be  determined  as 
well  as  the  mean  and  variability.  The  mean  heat  transfer 
rate  is 


5  =  .£ 


(estimate  of  mean  heat  transfer 


- ,  <Ki)  (  i  ! 

■ N*  ^  ter ' 


V  / 


(6.10) 


,  ,f*2a 

“i  u. 


The  q^  values  are  obtained  by  randomly  choosing  a  value  for 
each  of  the  random  variables  (r^r,H,K, A,L' (A)^).  The  random  val¬ 
ues  (K^,  rf  ^i^i^1  (A) r^)  are  chosen  by  entering  the  cumulative 
distribution ' functions  of  each  of  these  parameters  using  a  random 
number  generation  technique.  The  procedure  for  generating  a 
random  variable  from  a  particular  distribution  is  shown  in 
Figure  6.1.  The  cumulative  distribution  functions  are  deter¬ 
mined  from  experimental  data  using  Equations  2.55  and  2.56. 

The  estimate  of  the  heat  transfer  rate  standard  deviation  is 
obtained  from  Equation  2.56  and  is  rewritten  in  Equation  6.11. 


n  (q±  -  Q) 2 

S  =  £  - j -  estimate  of  variability  (6.11) 

Q  1=1  n”  (standard  deviation) 

The  parent  distribution  of  q^  is  hypothesized  using  a  goodness 

of  fit  test.  This  technique  becomes  accurate  for  large 

and  moderate  coefficients  of  variation. 

b)  Partial  Derivative  Method 

The  mean  heat  transfer  rate  is  determined  by  introducing 
the  means  of  the  random  variables  into  the  deterministic  design 
equation  and  solving  for  Q. 
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(6.12) 


Q  =  n£  r 


K' 


u 


l» 1  (A)  -■  fr 


p^g(H) 

'  v~ 


) 


The  standard  deviation  is  given  by  an  approximation  using  a 


Taylor  expansion 


(6.13) 


(6.14) 


The  distribution  of  Q  cannot  be  determined  through  significance 
tests  and  must  be  assumed  or  calculated  in  closed  form. 


c)  Algebra  of  Moments 


1  .• 


This  technique,  introduced  by  Haugen,  is  an  exact  tech¬ 
nique  for  determining  Q,Sg.  The  mean  is  calculated  as  in  the 
partial  derivative  method 

(all  random  variables  are  independent) 

The  standard  deviation  is  formulated  ixsing  algebra  of  the  first 
and  second  moments.  For  the  functions 


2,1  =  XY 


110 


sz,  /x^s  +  y2s  ^  +  s  2s  2 
l  V  y  x  x  y 


Z2  =  X  -  Y 


These  formulae  may  be  used  in  succession  to  obtain  S  .  The 

q 

distribution  is  determined  in  the  same  way  as  the  partial 
derivative  method. 

The  computer  program  block  diagram  is  shown  in  Figure 
6.1,  and 'the  actual  listing  appears  in  the  appendix.  The 
computer  program  is  simply  the  iteration  procedure  mentioned 
in  Chapter  2  using  the  three  solution  techniques  of  this 
chapter . 


START 


K,Sk 

A,  S 

r£  ,  S 
fr  rfr 

L'  , 

rr  ,sr 

H ,  S 

r  rr 

Radius  of  vapor 
passage,  S 

rv 

Heat  pipe  dimensions 
Fluid  density 
Fluid  latent  heat 
Fluid  surface  tension 
Fluid  viscosity 


WRITE 

All  the  above 

parameters 

CALCULATE  N, 


DO  down  to  and  including 
statement  100  M  times 


DO  down  to  and  including 
statement  40  n  times 


► 


* 


Figure  6.2.  (continued) 
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Figure  6,2P  (continued) 
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CHAPTER  7 


THE  PROBABILISTIC  DESIGN  EQUATION 


The  probabilistic  model  was  mentioned  in  Chapter  5  and 
is  written  below  tor  the  recessed  condition 


<K.-Sk}  /  2 

(AS  )  =  N  - i. - 

Q  *  (L*  (A)  'ST  ,  \ *  r  ..  ,S  ' 

'  'r  L  (A;r  :r'  rfr 


(7.1) 

The  assumptions  are  the  same  as  the  deterministic  equation 
with  the  addition  of  the.  following: 

a)  A,  K,  r (A)r  and  H  are  random  variables  and  are 

described  by  simple  two  parameter  continuous 
distributions . 

b)  Variabilities  of  all  other  parameters  are  considered 
small  and  are  therefore  deterministic 

c)  All  random  variables  are  independent 

d)  Coefficients  of  variation  (x/s  )  are  large  (>7)  for 
denominator  random  variables. 


7 . 1  Comparison  of  Probabilistic  Model  with  Experimental  Data 
To  test  the  validity  of  the  probabilistic  model,  fifty 
heat  pipe  wicks  were  manufactured,  tested,  and  their  proper¬ 
ties  recorded  in  the  Append  .x  and  analyzed  in  Chapter  5.  The 
following  analysis  will  examine  data  set  2  in  detail  since 
this  data  set  is  comp'osed  of  thirty  tests  and  carries  the 
most  significance.  The  analysis  for  the  other  data  sets  are 
identical  and  will  be  mentioned  throughout  the  discussion. 

For  each  of  the  fifty  heat  pipe  wicks  tested,  K,  A,  L* ,  and 
r^r  were  recorded.  A  steady  state  heat  transfer  rate  was 


115 


established  and  the  evaporator  section  was  raised  until  the 
wick  began  to  dry  out.  The  steady-state  Maximum  heat  trans¬ 
fer  rate,  Qob'  was  recorded  along  with  the  gravity  effect,  H, 
at  dry  out.  The  experimental  values  of  K,  A,  L1,  rfr,  and  H 
were  inserted  into  the  recession  design  equation,  Equation  2.57, 
and  Qca  was  calculated.  The  deviation  of  the  observed  experi¬ 
mental  mean  heat  transfer  rate  from  the  calculated  mean  heat 
transfer  rate,  Qcq  -  QQb,  was  calculated  for  each  test.  Table 

7.1  shows  the  mean  and  standard  deviation  of  Q  -  o  .  .  The 

ca  ob 

mean  difference,  Q^a  "  Qob'  was  small  in  all  cases  when  com¬ 
pared  to  the  standard  deviation.  We  will  test  this  hypothesis 

17 

using  the  following  Student’s  T  Statistic.. 

Tcalc.  -  1 . P  I*1” V2  -  T  “/2  <»-«  0.2) 


(7.2) 


Q  -  Q  . 
n  ca .  ob . 

d  =  y  — - - - 

i=i  n 

If  the  calculated  T  value,  Tcalc.,  is  less  than  the  critical  T 
value  at  the  a  =  .005  significance  level,  we  will  accept  the 
hypothesis  that  there  is  no  significant  difference  between  Q 
and  Qot>.  Table  7.1  shows  the  results  of  the  hypothesis  test. 

At  every  heat  transfer  rate,  the  mean  observed  heat  transfer 
rate  was  not  significantly  different  from  the  mean  calculated 
heat  transfer  rate.  This  indicates  that  Equation  2.57  is  a 
good  predictor  of  the  mean  heat  transfer  rate.  Figures  7.7 
through  7.11  show  plots  of  the  data  from  Table  7.2  and  a  com¬ 
parison  of  the  mean  maximum  calculated  heat  transfer  rate  with 
experimental  data  for  each  of  the  data  sets.  Figures  7.1 
through  7.3  show  the  distribution  of  Q„  -  Q  ,  ,  which  is  called 
the  difference  distribution.  Appendix  A  shows  values  of  Q 

Ca 

and  Sq  for  all  the  heat  pipes  tested. 
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Probability  Probability  Probability 


Q  -  Q  ,  Watts 
ca  ob 
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Figure  7.3.  Difference  Distribution  for  Data 
Set  2  (100-130  watts)  Using  Reces 
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These  values  were  calculated  from  Equation  2.57  using  the 

simulation  technique.  200  random  variables  were  calculated 

for  each  simulation  and  the  resultant  Q  distribution  did 

ca 

not  reject  the  Kolmogorov  Smirnov  significance  test  at  95% 
significance.  The  distribution  of  Q  was  hypothesized  as 

CcL 

normal.  The  uncertainty  bound  on  the  predicted  heat  trans¬ 
fer  rates  for  each  pipe  will  be  plus  or  minus  three  standard 
deviations  from  the  mean  since  these  limits  will  encompass 
99.7%  of  the  Q  Normally  distributed  random  variables.  If 
the  observed  heat  transfer  rate  lies  in  this  region  of  plus 
or  minus  three  standard  deviations  of  Q  ,  the  probabilistic 
design  equation  is  credited  as  having  predicted  the  occurrence 
of  the  experimental  result,  Q  ^  If  Qob  falls  outside  the 
three  standard  deviation  bounds,  the  probabilistic  design 
equation  will  be  considered  inadequate  in  the  prediction  of 
the  occurrence  of  Q^.  Referring  to  the  data  of  Appendix  B, 
each  of  the  135  observed  heat  transfer  rates,  Qq^,  were  within 

the  three  standard  deviations  of  Q  .  If  we  assume  that  the 

ca 

three  standard  deviation  bound  on  the  variability  of  the 
calculated  heat  transfer  rate  to  be  correct,  then  the  proba¬ 
bility  of  observing  Q  ^  within  these  bounds  will  be  .9972. 

From  135  tests,  all  the  observed  heat  transfer  rates  were 
within  these  bounds.  The  probability  of  observing  135  ob¬ 
served  heat  transfer  rates  within  the  three  standard  devia- 

135 

txon  bound  and  0  outside  the  bound  is  (.9972)  or  .668. 

This  indicates  a  high  probability  that  the  calculated  vari¬ 
ability  of  the  observed  heat  transfer  rate,  S^,  is  correct 


No.  of  Occurrences  of  Q  .  No.  of  Occurrences  of  Q  .  No.  of  Occurreno 3&  of  Q 


Figure  7.4. 
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Standard  Deviations  from  the  Mean 


Figure  7.5. 
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Figure  7.6. 


Distribution  ox  GCa 
about  Qq'd  for  100- 
Mesh  Stainless  Steely 
Two  Layer  Wicks  at  a 
Mean  Wattage  of  50.566 


Distribution  of 
about  Qca  for  100- 
Mesh  Stainless  Steel, 
Two  Layer  Wicks  at  a 
Mean  Wattage  of  81.333 


Distribution  of  Q0k 
about  Qca  for  100- 
Mesh  Stainless  Steel, 
Two  Layer  Wicks  at  a 
Mean  Wattage  of  116.799 
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©  Experimental  Data 


Figure  7.9.  Comparison  of  the  Mean,  Maximum  Heat 
Transfer  Rate  with  Experimental  Data 
from  Data  Set  3  (200  mesh  stainless 
steel  3  layer  wicks) 


Maximum  Heat  Transfer  Rate,  Q  (watts) 

'  max 


Figure  7.10.  Comparison  of  the  Mean,  Maximum  Heat 
Transfer  Rate  with  Experimental  Data 
from  Data  Set  4  (50-mesh  stainless 
steel  two  layer  wicks) 


according  to  data  from  135  tests.  Figures  7.4,  7.5,  and  7.6 

show  the  distribution  of  the  observed  heat  transfer  rate,  Q  ,  , 

'  *cb' 

about  the  mean  calculated  heat  transfer  rate,  Q  for  data 

cs 

set  2.  The  distribution  of  Figure  7.4  did  not  reject  the 
Komogrov-Smirnov  significance  test  when  tested  for  Normality 
and  neither  did  the  distributions  of  Figures  7.5  and  7.6. 
However,  as  heat  transfer  rate  increased,  the  occurrence  of 
Qob  aPPeared  more  likely  at  the  tails  of  the  distributions 
and  the  resulting  distributions  became  flatter  than  expected. 

The  consistency  of  predicting  the  range  of  occurrence  of 
^ob  us^n9  Equation  2.56  is  evident.  However,  some  basic 
numerical  problems  may  result  in  applications  to  low  heat 
transfer  rates  at  high  gravity  effect  and  also  to  low  heat 
transfer  rates  at  high  permeabilities.  Upon  examination  of 
our  mathematical  model  below,  we  can  readily  see  how  problems 


Q  =  N, 


‘  K  / 
L*  (A)r  \l 


p„gH> 


arise.  Shown  below  are  two  sample  calculations  using  data 
from  data  set  2.  One  calculation  has  a  high  gravity  head 


Q  =  N, 


Q  =  N, 


—  f  24 ‘  _  62.4  (2.96)  \ (Low 

r  '4.069  x  10~3  4.4  x  10“3(12) /  Gravity 

Loss) 

-  (5900  -  3500) 


(7.3) 


Q  =  N  yr-^TT-  f 5900  - — V4- 

1  '  4.4  x  10“3(12)  /  Gravity 


(7.4) 


LOSS) 


=  NJI  lT(aT  ’  (5900  ■  5000> 

V  >  r 


loss  while  the  other  has  a  low  gravity  loss.  Notice  at  the 
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high  gravity  head  loss  the  subtraction  of  the  two  rather  large 
numbers.  This  shows  the  classic  numerical  problem  of  loss  of 
significance  when  subtracting  two  nearly  identical  numbers. 

Such  a  problem  was  discussed  in  Chapter  6  and  if  one  refers 
to  Table  6.3  of  Chapter  6  the  sensitivity  of  this  computation 
to  error  at  relatively  small  standard  deviations  can  be  realized. 
Input  standard  deviations  for  the  surface  tension  term  are  on 
the  order  of  5%  and  the  gravity  terms  are  approximately  1%. 

As  shown  on  the  previous  page,  these  particular  values  of  stan¬ 
dard  deviation  lead  to  problems  at  low  heat  flow  since  the  sur¬ 
face  tension  and  gravity  terms  are  on  the  order  of  magnitude  of 
5000.  If  one  refers  to  Chapter  6,  Table  6.3,  we  can  observe  a 
reliable  prediction  of  the  mean  of  the  difference  function. 
Despite  the  high  resultant  standard  deviation,  the  error  in  the 
mean  (calculated  by  simulation)  is  small  in  comparison  (about 
10%) .  We  may  conclude  that  although  low  heat  flow  (near  wick 
dry  out)  standard  deviation  calculations  may  be  unreliable, 
the  mean  heat  flow  calculation  is  quite  good  as  indicated  by 
the  data  discussed  earlier  in  this  chapter.  If  calculations 
must  be  made  in  the  low  heat  flux  range  (less  than  50  watts 
and  near  burn  out) ,  the  resultant  calculated  standard  devia¬ 
tions  of  heat  flow  will  probably  be  conservative.  Another 
numerical  problem  results  at  high  permeabilities  and  low  crit¬ 
ical  heat  transfer  rates.  This  is  exemplified  in  data  set  4 
of  the  appendix  (50  mesh  stainless  steel  wicks)  .  The  permer- 
ability  was  large  and  the  capillary  pumping  term  was  small 
resulting  in  a  loss  of  significance  in  the  difference  term. 
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This  resulted  in  extremely  large  calculations  of  standard  de¬ 
viation  and  seemed  too  conservative  as  indicated  by  the  data. 

We  have  been  examining  the  prediction  of  heat  transfer 
rate  knowing  the  properties  of  the  particular  wick  tested. 

If  one  wishes  to  design  a  heat  pipe  using  the  manufacturing 
technique  of  Chapter  3,  then  the  wick  properties  of  Table  4.23 
should  be  used.  These  property  values  give  the  range  of  occur¬ 
rence  of  wick  properties  for  a  particular  heat  pipe  manufac¬ 
tured.  Table  7.2  shows  the  mean  calculated  wattage  and  stan¬ 
dard  deviation  at  various  values  of  H,  using  the  three  tech¬ 
niques  mentioned  in  Chapter  6.  The  wattage  values  of  Table  7.2 
give  the  range  of  maximum  heat  transfer  rate  one  might  expect 
if  a  heat  pipe  was  designed  and  operated  at  the  particular 
gravity  effect,  H.  Notice  that  the  calculated  standard  devia¬ 
tions  are  quite  high  due  to  the  high  standard  deviations  of 
the  input  properties  of  Table  4.23.  The  three  techniques  used 
to  solve  the  probabilistic  design  equation  were  quite  close. 

The  simulated  results  were  slightly  higher  than  the  rest  due 
to  high  variability  in  the  denominator  random  variables.  Now 
that  the  variability  of  the  manufacturing  process  has  been  des¬ 
cribed,  its  applicability  to  design  will  be  discussed. 

7.2  Results  of  the  Probabilistic  Design 

Assume  that  our  design  specification  required  that  a  par¬ 
ticular  heat  pipe  of  this  design  operate  at  75  watts  at  a  grav¬ 
ity  head  of  2.96"  H20.  Referring  to  data  set  2,  the  results  of 
the  probabilistic  model  (H  =  2,96"  H20)  yields 
Mean  heat  flow  125.707  watts 
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os 


.Notation:  Qcaicr  Sim-Calculated  wattage  using  simulation,  SgC,  sim  -  standard  deviation 
of  calculated  wattage  using  simulation  Qcalc/  P.D.&  Exact  -  Calculated  wattage 
using  partial  derivative  and  exact  techniques . (algebra  of  moments  technique). 


Standard  deviation  24.678  watts 


Distribution  Normal  (passes  Kolmogorov 

significance  test) 

Assuming  the  estimates  of  the  mean  and  standard  deviation  to 
be  perfect,  the  probability  of  failure  for  this  heat  pipe  at 
75  watts  is 


125.707-75 

- 247578"' 


=  2.06  =  .98 


1  -  .98  =  probability  of  failure 

=  .02 


Since  30  samples  were  required  for  the  estimation  of  the 
mean  and  standard  deviation  of  the  heat  flows  these  estimates 
.  are  by  no  means  perfect.  Due. to  the  randomness  of  these  es¬ 
timates  we  will  quote  a  maximum  probability  of  failure  at 
some  confidence  level.  Using  tolerance  theory,  we  may  con¬ 
clude  that  the  upper  tolerance  limit  of  the  probability  of 
failure  for  this  heat  pipe  design  is  .10  at  a  confidence  of 
99%,  or  that  we  are  99%  confident  that  the  probability  of 
failure  does  not  exceed  .10.  These  tolerance  limit  values 
were  computed  from  a  table  based  on  a  given  mean,  standard 
deviation  and  sample  size.  This  result  indicates  that  the 
maximum  amount  of  failures  one  might  expect  from  100  units  of 
this  design  is  10  units.  If  this  probability  of  failure  is 
unacceptable,  we  can  decrease  the  probability  of  failure  by 
changing  design  parameters  such  that  the  meam  heat  flow  is 
larger.  For  example,  we  can  enlarge  the  wick  area  and  one 
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might  calculate  a  heat  flow  of  150.00  watts  at  a  standard 
deviation  of  25  watts  with  a  maximum  probability  of  failure 
of  .025  or  175  watts  at  a  standard  deviation  of  30  watts 
with  a  .005  probability  of  failure.  The  difference  between 
the  mean  calculation  and  the  specification  will  indicate  the 
risk  one  takes  in  designing  a  certain  heat  flow.  The  design 
distribution  can  be  manipulated  to  minimize  the  risk  of  pro¬ 
ducing  below  specification.  The  reliability  of  this  design 
is  an  initial  start  up  reliability.  An  increase  in  failure 
probability  may  occur  due  to  gas  generation  in  the  heat  pipe. 
The  resulting  gas  may  block  the  condenser. 


7.3  Heat  Pipe  Temperature  Drop 

The  heat  pipe  thermal  analog  circuit  appears  in  Figure 

7.12.  Applying  this  analysis  to  our  low  mass  flow  water  heat 

pipes,  we  are  allowed  to  neglect  the  evaporation  temperature 

loss  (T  -  T  ) f  the  vapor  flow  loss  (T  -  T  )  and  the  con- 
ew  ev  ev  cv 

densation  temperature  loss  (T  -  T  ) .  These  conclusions  were 

cv  cw 

a  result  of  inserting  the  maximum  expected  mass  flow  rate’,  ( .  7 
lbm/hr)  into  the  pressure  loss  equation,  7.5  (derived  from 
kinetic  theory)  and  the  Clausius  equation,  7.6. 


®  ■  5^rL  "here  R  ’  Constant1 


M 


M  =  Molecular  Weight 


AT  =  T  AP/phfgJ 


(7.5) 


(7.6) 


Figure  7..  12  shows  the  general  thermal  analog  circuit  and  our 
thermal  circuit  is  reduced  to 


T  T  . 

eo  A  ei 

a. .  — /vs/v  —  * - v/V*V“ 


T  T  . 

ew  ci 

— * — vW — 3 — s/\/\r 


co 
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ep 
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ew 


cw 


R 


cp 
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Heat  Addition  Q  I-Ioat  Removal  Q 
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where 


ep 


=  £n  ro^ri 


.  x .  /r 

R.  =  £n  i7  v 

^  2rk  l 

'  w  c 


xi/rv 

Rew  ~  £n  2iE 


w^E 


r  /r . 
Rcp  =  £n  2tiK_L 


p  c 


(7.8) 


(7.9) 

(7.10) 


This  model  should  adequately  describe  temperature  loss  assum¬ 
ing  the  fully  wetted  wick.  For  the  fully  recessed  situation 
near  burn  outf  the  evaporator  wick  resistance  is  modified  as 
follows 


^PEW  x  rfp,i/rprw  x 

^RW  =  *n  V(ro  +  1/4(ri  -  ro>)/2T,k«<Le  -  Lfr> 

"fRW  =  *n  V(ro  +  ’41  (ri  '  ro)/2"kl(Lfr) 

It  is  assumed  that  the  fully  and  partially  recessed  areas  of 
the  evaporator  wick  section  combine  in  parallel  to  form  a 
total  resistance. 

The  thermal  conductivity  of  the  wick-fluid  matrix  is  dif¬ 
ficult  to  determine,  but  may  be  approximated  by 

kw  =  +  (1  -  e)ks  (7.11) 

assuming  that  the  wick  is  well  bonded  to  the  pipe  wall  (little 
wick-pipe  thermal  resistance) .  This  equation  is  inaccurate  in 
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our  particular  case  because  of  the  poor  pipe-wick  contact  and 
oxide  coatings  which  yield  a  high  thermal  resistance  as  shewn  in 
Figure  7.13  „  it  is  therefore  assumed  that  the  wick  contact 
resistance  with  the  wall  is  on  the  same  order  of  magnitude  as 
that  of  the  fluid  and  thermal  conductivity  of  the  wick  is 
essentially  that  of  the  liquid. 


Figure  7.13.  Wick  Contact  with  Pipe  Wall 


As  shown  in  Table  7.3  the  stainless  steel  wick  tempera¬ 
ture  loss  predictions  (using  conductivity  of  the  fluid  as  the 
wick  conductivity)  were  always  slightly .higher  than  the  ob¬ 
served  temperature  losses.  The  stainless  steel  wick  has 
sufficient  rigidity  to  allow  some  parallel  heat  paths  to  exist. 
However,  the  above  assumption  (conductivity  of  the  wick  equal 
to  conductivity  of  the  fluid)  yields  a  conservative  result  and 
is  more  indicative  of  reality  than  the  predictions  of  Equation 
7.11. 

The  copper  wicks,  on  the  other  hand,  seemed  to  be  described 
accurately  by  the  fluid  conductivity  assumption.  It  is  theorized 
that  the  copper  wick  lacks  the  rigidity  for  even  mediocre  thermal 
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TABLE  7.3 

TEMPERATURE  LOSS  PREDICTION  (Recession.  Model) 


Data 

Set 

Wattage 

Obs. 

DT 

°F 

S:  Obs. 

DT 

Op 

Calc. 

DT 

op 

S  Calc. 

DT 

°F 

1 

42-59 

8.309 

3.432 

8.2S7 

* 

3.177. 

1 

70-90 

9.363 

3.293 

15.046 

5,046 

1 

.  100-130 

12.799 

3.271 

21.939 

8.386 

2 

42-59 

6.799 

1.788 

7.163 

2.608 

2 

70-90 

9.033 

1.564 

12.438 

3.132 

2 

100-130 

12.099 

2.411 

18.835 

3.803 

3 

20-30 

2.000 

0,707 

2.582 

1.235 

3 

60-80 

4.799 

1.923 

8.070 

1.743 

4 

50-60 

7.399 

1.341 

16 . 3£2 

21.767 

4 

100-130 

14.600 

2.191 

41.899 

23.400 

5 

45-55 

7.400 

1.516 

7.789 

1.638 

5 

90-120 

16.400 

1.949 

18.980 

2.352 

Note:  a)  Obs  DT  -  mean  observed  temperature  loss 

b)  S  obs  DT  -  standard  deviation  of  observed 
temperature  loss 

c)  Calc.  DT  -  calculated  mean  temperature  loss 
using  recession  model 

d)  S  calc.  DT  -  standard  deviation  of  the  cal¬ 
culated  temperature  loss 


contact  with  the  pipe  wall  and*  that  essentially  no.  parai-iex 
.heat  flow  paths-  exist. 

Figure  7,.  14  shows  a  typical  wire  mesh  heat  pipe'  tempera¬ 
ture  distribution.  Refer  to  the  appendix  data  sets  for  vapor 
temperature  and  temperature  loss. 


CHAPTER  8 


CONCLUSIONS 

8 . 1  Probabilistic  Design 

The  probabilistic  approach  to  design,  is  illustrated  in 
this  work.  For  any  given  deterministic  design  equation, 

A  =  F(x^,  x2, . .  .^5 

that  describes  the  behavior  of  some  physical  phenomenon,  we 
may  transform  into  the  probabilistic  domain  where 

RA  -  F (R^  ,  R^  Rx  ) 

1  /'2  n 

R  indicates  the  particular  variable  to  be 
random  and  described  by  some  probability 
function 

The  random  result  of  the  probabilistic  design  equation,  RA, 
indicates  the  uncertainty  of  the  design  quantity  to  be  pre¬ 
dicted  as  a  function  of  the  input  parameter  uncertainties. 

The  value  of  the  probabilistic  design  is  that  it  eliminates 
the  requirement  for  safety  factors  (which  are  arbitrarily 
conceived  and  may  not  describe  the  situation  at  hand)  in 
favor  of  a  quantitative  measure  of  the  range  of  a  design  re¬ 
sults.  One  can  therefore  indicate  the  probability  of  manu¬ 
facturing  a  successful  design  hence  predict  it’s  reliability. 

Three  methods  for  the  solution  of  the  probabilistic  design 
equations  are  mentioned.  They  are 

a)  Simulation 

b)  Algebra  of  moments 

c)  Taylor  approximation 
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These  methods  work  well  in  prediction  of  most  design  function 
variabilities  but  caution  must  be  exercised  in  the  handling 
of  certain  functions  such  as  RA-RB  and  1/RA  at  large  coeffi¬ 
cients  of  variation. 

8.2  Heat  Pipe  Design 

The  recession  model  is  quite  reliable  in  the  prediction 
of  the  mean  heat  flow  capability  of  the  wire  mesh  wick  heat 
pipe.  The  prediction  of  heat  flow  standard  deviations  was 
good  at  high  heat  flows  (therefore  low  gravity  components) . 

At  low  heat  flows  the  magnitude  of  the  gravity  and  surface 
tension  terms  approach  each  other  to  form  the  classical 
numerical  problem  mentioned  in  Chapter  6.  At  these  low  heat 
flows,  prediction  of  the  standard  deviation  is  conservative. 
The  recession  model  is  by  no  means  perfect  and  the  descrepancy 
of  standard  deviation  predictions  at  low  heat  flows  is  in¬ 
herent  to  this  particular  model. 

Uncertainty  in  the  description  of  heat  pipe  wicking  prop¬ 
erties  is  discussed.  A  "rule  of  thumb"  for  predicting  under- 
tainties  in  heat  pipe  properties  is  listed  below  in  Table  8.1. 
The  variability  of  these  properties  are  a  result  of  experimen¬ 
tal  and  manufacturing  variability.  It  should  be  mentioned 
that  the  property  data  generated  in  this  work  describes  the 
wicks  manufactured  by  the  process  mentioned  in  Chapter  3. 

Other  types  of  wire  mesh  designs  may  have  significantly  dif¬ 
ferent  properties. 
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TABLE  8.1 

UNCERTAINTY  OF  WICK  PROPERTIES 


AS  A  %  OF  THE 

Property 

Permeability,  K 

Critical  radius,  re 

Wick  cross  sectional 
area,  A 

Porosity,  e 

Tortuosity,  b 


MEAN  PROPERTY 

Standard  Deviation 
in  %  of  Mean 

10% 

5% 

3% 

3% 

10% 


* 


i 

l 
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APPENDIX  A 


BASIC  RAW  DATA 

DATA  NOMENCLATURE 


The  following  is  an  explanation  of  data  format. 


K 

Wick  permeability  as  described  by  Darcy 

(Ft2) 

SK 

Standard  deviation  of  permer±>ility 

(Ft2) 

RE 

Effective  capillary  critical  radius  based 
on  two  layers  of  screen 

(inches) 

SR 

Standard  deviation  of  critical  radius 

(inches) 

A 

Flow  area  of  the  wick 

(Ft2) 

SA 

Standard  Deviation  of  flow  area 

(Ft2) 

I  LAYER 

Capillary  pressure  capability  of  one 
layer  of  wire  mesh 

(in. 

h2o) 

II  LAYER 

Capillary  pressure  capability  of  two 
layers  of  wire  mesh 

(in. 

h2°) 

SUB  TEST 

FORMAT 

VAPOR 

TEMP 

Heat  pipe  vapor  temperature  at  equi¬ 
librium 

( °F) 

H 

Total  gravity  head  loss  experienced 
by  heat  pipe 

(in. 

h2o) 

SH 

Standard  deviation  head  loss 

(in. 

h2o) 

ALPP 

Darcy  effective  flow  length  (L1)  for 
non  recession  model 

if  _  yL’m 

paSp 

(Ft) 

SRALPP 

Standard  deviation  of  Darcy  effective 
flow  length  non  recessed  model 

(Ft) 

Note: 


Each  of  the  standard  deviations  are  a  result  of 
experimental  error  in  determination  of  the  parti 
cular  properties 


/V3 


Al 


ALPPR 


SRALPR 

RL 


SRL 

CALC 
WATTS  R 

SWR 

CALC 

WATTS 

SW 

OBS 

WATTS 

DT  CALC 

SDT 

OBS  DT 


Darcy  effective  flow  length  for  recession 
model 

Standard  Deviation  for  Darcy  effective 
flow  length  recession  model 

Distance  from  the  evaporator  edge  to 
the  beginning  of  the  fully  recessed 
region  of  the  wick 


Standard  deviation  for  RL 

Heat  pipe  wattage  capability  using  the 
recession  model  Eq(  ) 

Standard  deviation  of  CALC  WATTS  R 

Heat  pipe  wattage  capability  using  non 
recessed  model  Eq(  ) 

Standard  deviation  of  CALC  WATTS 

Experimentally  observed  wattage  capability 
of  a  heat  pipe  with  the  above  mentioned 
properties 

Calculated  temperature  drop  from  evapora¬ 
tor  surface  to  condenser  surface.  All 
vapor  temperature  drops  were  negligible 

Standard  deviation  of  DT  CALC 

Experimentally  observed  temperature  drop 
along  a  heat  pipe  with  the  above  mentioned 
properties 


A2 


(Ft) 

(Ft) 

(Ft) 

(Ft) 

(Watts) 

(Watts) 

(Watts) 

(Watts) 

(Watts) 

( °F) 
(°F) 

(°F) 


test  data  sheet  for  data  set  #1 


TEST  # 


17 


18 


19 


K 

SK 

RE 

SR 

A 

SA 

I  LAYER 

II  LAYER 


3.5840  E-09 
1.9260  E-10 
4.0550  E-03 
1.0000  E-04 
2.2830  E-04 
5.0400  e-06 
5.0000 
5.6250 


6.6540  e-09 
.2.1580  E-10 
4.8090  E-03 
1.0000  E-04 
2.3540  e-04 
6.2900  e-G6 
4.6250 
4.7500 


5.2460  2-09 
2.8740  E-10 
4.2470  E-03 
1.0000  E-04 
2.3750  E-04 
5.2000  E-Ofc 
4.1250 
5.3640 


(42-59W) 
SUB  TEST) 


VAPOR  TEMP 

H 

SH 

ALPP 

SRALPP 

ALPPR 

SRALPR 

RL 

SRL 

CALC  WATTS  R 
SWR 

CALC  WATTS 
SW 

OBS  WATTS 
DT  CALC 
SDT 

OBS  DT 


140 

140 

146 

3.0600 

3.2700 

3.7400 

0.0417 

0.0417 

0.0417 

1.0400 

1.0400 

1.0400 

0.0300 

0.0300 

0.0300 

1.3260 

1.3380 

1.3430 

0.0651 

0.0769 

0.0700 

1.0410 

1.2080 

0.7910 

0.3210 

0.4500 

0.2890 

48.1120 

49.0530 

42.4140 

5.3890 

7.9690 

5.6810 

60.0870  . 

60.8100 

54.0170 

5.7640 

8.4360 

6.3500 

55.0000 

52.0000 

53.0000 

8.3830 

8.6460 

7.1740 

1.0700 

1.6870 

1.0530 

10.0000 

7.0000 

8.0000 

(CONTINUED) 


A3 


TEST  DATA  SHEET  CONTINUED  FOR  DATA  SET  #  1 


(70  -  90W 
SUB  TEST) 


VAPOR  TEMP 

H 

SH 

ALPP 

SRALPP 

ALPPR 

SRALPR 

RL 

SRL 

CALC  WATTS  R 
SWR 

CALC  WATTS 
SW 

OBS  WATTS 
DT  CALC 
SDT 
OBS  DT 


165 

2.500GT 

0.0417 

1.0400 

0.0300 

1.3380 

0.0748 

1.1240 

0.4270 

67.0410 

7.2360 

83.7070 

8.0990 

78.0000 

11.9080 

1.4090 

15.0000 


166 

2.2500 

0.0417 

1.0400 

0.0300 

1.3600 

0.0965 

1.2080. 

0.5760 

101.8750 

14.4850 

127.7850 

18.4170 

84.0000 

18.2870 

3.0880 

9.0000 


182 

2.7600 

0.0417 

1.0400 

0.0300 

1.3440 

0.0722 

0.8330 

0.2830 

85.6470 

9.2860 

100.4470 

12.4160 

80.0000 

14.7720 

1.7190 

8.0 


(100  -  130) W 

SUB  TEST) 

• 

VAPOR  TEMP 

198 

185 

204 

H 

1.8000 

1.8700 

2.0800 

SH 

0.0417 

0.0417 

0.0417 

ALPP 

1.0400 

1.0400 

1.0400 

SRALPP 

0.0300 

0.0300 

0.0300 

ALPPR 

1.4140 

1.4280 

1.3720 

SRALPR 

0.1300 

0.1490 

0.0696 

RL 

1.2080 

1.2080 

0.8540 

SRL 

0.7490 

0.8100 

0.3270 

CALC  WATTS  R 

96.3000 

136.1640 

122.2110 

SWR 

12.0020 

22.0540 

12.0210 

CALC  WATTS 

118.9330 

169.8960 

155.1590 

SW 

10.6320 

26.3460 

16.6560 

OBS  WATTS 

110.0000 

110.0000 

100.0000 

DT  CALC 

17.1850 

24.3230 

20.9000 

SDT 

2.2210 

4.2000 

2.0280 

OBS  DT 

18.0000 

13.0000 

10.0000 

i4k 


TEST  DATA  SHEET  FOR  DATA  SET  #1 


TEST  2 

20 

21 

K 

5. .2440 

E-09 

5.4870 

E-09 

SK 

2.6850 

E-10 

1.4950 

E-10 

RE 

4.9060 

E-03 

4.3660 

E-03 

SR 

1.0000 

E-04 

1.0000 

E-04 

A 

2.5030 

E-04 

2.4450 

E-04 

SA 

1.1040 

E-05 

4.7900 

E-05 

I  LAYER 

4.2500 

4.5000 

II  LAYER 

4.6250 

5.2500 

(42 -59$) 
SUB  TEST) 


VAPOR  TEMP 

H 

SH 

ALPP 

SRALPP  ' 

ALPPR 

SRALPR 

RL 

SRL 

CALC  WATTS  R 
SWR 

CALC  WATTS 
SW 

OBS  WATTS 
DT  CALC 
SDT 

OBS  DT 


138 

2.9600 

0.0417 

1.0400 

0.0300 

1.3340 

0.0645 

1.0410 

0.3080 

48.1020 

5.3550 

60.1170 

5.6620 

43.0000 

7.5660 

0.9680 

6.0000 


147 

3.7500 

0.0417 

1.0400 

0.0300 

1.3260 

0.0648 

0.9990 

0.3070 

41.9380 

4.9170 

52.9920 

5.7160 

43.0000 

5.4070 

0.7900' 

3.0000 


(CONTINUED) 


A5 


TEST  DATA  SHEET  CONTINUED  FOR  DATA  SET 


(70  -  90W 
SUB  TEST) 


VAPOR  TEMP 

H 

SH 

ALPP 

SRALPP 

ALPPR 

SRALPR 

RL 

SRL 

CALC  WATTS 
SWR 

CALC  WATTS 
SW 

OBS  WATTS 
DT  CALC 
SDT 
OBS  DT 


178 
2.3900* 
0.0417 
1.0400 
0.0300 
1.3630 
0.0749 
1.1660 
0.4370 
R  78.5490 
8.9850 
97.8670 
10.9480 
87.0000 
12.6340 
1.6100 
8.0000 


157 

2.7300 

0.0417 

1,0400 

0.03.00 

1.3310 

0.0680 

0.9990 

0.3140 

83.5200 

7.0820 

104,2500 

7.6330 

73,0000 

12.9560 

1.1420 

7.0000 


(100  -  130) W 
SUB  TEST) 


VAPOR  TEMP 

H 

SH 

ALPP 

SRALPP 

ALPPR 

SRALPR 

RL 

SRL 

CALC  WATTS 
SWR 

CALC  WATTS 
SW 

OBS  WATTS 
DT  CALC 
SDT 

OBS  DT 


192 

1.6000 
0.0417 
1.0400 
0.0300 
1.4110 
0.1300 
1.1660 
0.7490 
R 120.1130 
14.5570 
149.8650 
15.3920 
128.0000 
18.2170 
1.7020 
13.0000 


192 

1.9800 

0.0417 

1.0400 

0.0300 

1.3490 

0.0707 

1.0830 

0.4390 

129.5810 

10.555 

161.2640 

9.3140 

104.0000 

'20.0970 

1.5970 

10.0000 


TEST  DATA  SHEET  FOR  DATA  SET  #2 


TEST  # 

22 

23 

24 

K  ■ 

6.909 

E-09 

7.479 

E-09 

6.850 

E-09 

SK 

4.862 

E-10 

2.634 

e-io 

3.748 

E-10 

RE 

4.139 

E-03 

3.786 

E-03'- 

3.645 

E-03 

SR 

1.000 

E-04 

1.000 

E-04 

1.000 

E-04 

A 

2.475 

E-04 

2.584 

E-04 

2.555 

E-04 

SA 

1.479 

E-05 

6.540 

E-05 

1.190 

E-05 

I  LAYER 

4.625 

5.500 

5.250 

II  LAYER 

5.500 

6.000 

6.250 

(42-59W) 
SUB  TEST) 


VAPOR  TEMP  156.000 

156.000 

155.000 

H 

43.300 

4.740 

4.650 

SH 

0.042 

0.042 

70.042 

ALPP 

1.040 

1.040 

1.040 

SRALPP 

0.030 

0.030 

0.030 

ALPPR 

1.324 

1.317 

1.323 

SRALPR 

.064 

0.067 

0.064 

RL 

.999 

1.187 

0.999 

SRL 

.2846 

0.270 

0.273 

CALC  WATTS 

K36.840 

44.557 

56.165 

SWR 

7.502 

8.735 

9.430 

CALC  WATTS 

46.516 

56.569 

71.074 

SW 

8.743 

10.500 

10.562 

OBS  WATTS 

46.090 

48.000 

48.500 

DT  CALC 

5.671 

7.185 

8.922 

SDT 

1.272 

1.470 

1.640 

OBS  DT 

5.000 

5.000 

7.000 

(CONTINUED) 

til 


A7 


TEST' 

DATA  SHEET 

CONTINUED  TOR  -DATA  SET  # 

2 

(70  ~  90W 

SUB  TEST) 

VAPOR  TEMP 

192.000 

169.000 

189.000 

H 

3.840‘ 

3. '760, 

3.890 

SH 

0 . 042 

0.042 

0.042 

ALPP 

1.040 

1.040 

1.040 

SRALPP 

0.030 

0.030 

0,030 

ALPPR 

1.327 

1.329 

1.330 

SRALPR 

0.067 

0.069 

0.068 

RL 

1.104 

1.208 

1.083 

SRL 

0.340 

0.374 

0.351 

CALC  WATTS  R 

62.268 

104.462 

101.620 

SWR 

10.270 

11.761 

13.471 

CALC  WATTS 

77.506 

129.862 

126.742 

SW 

13.294 

15.189 

17.383 

OBS  WATTS 

75.000 

82,000 

84.000 

PT  CALC 

9.894 

17.290 

16 .571 

SDT 

1.760 

2.024 

2.362 

OBS  DT 

10.000 

8.000 

8.000 

(100  -  130) W 

SUB  TEST) 

VAPOR  TEMP 

198.000 

194.000 

187.000 

H 

2.850 

3.130 

2.990 

SH 

0.042 

0.042 

0.042 

ALt‘P 

1.040 

1.040 

1.040 

SRALPP 

0.030 

0.030 

0.030 

ALPPR 

1.357 

1.348 

1.351 

SRALPR 

0.077 

0.069 

0.073 

RL 

1.083 

1.208 

1.020 

SRL 

0.404 

0.453 

0.430 

CALC  WATTS  R 

127.166 

151.798 

158.986 

SWR 

15.435 

14.060 

16.870 

CALC  WATTS 

158.640 

191.804 

199.124 

SW 

19.603 

18.773 

21.956 

OBS  WATTS 

101.000 

110.000 

114.000 

DT  CALC 

19.936 

25.203 

25.577 

SDT 

2.492 

2.299 

2.743 

OBS  DT 

12.000 

12.000 

10.000 

tsO 


A8 


TEST  DATA  SHEET  FOR  DATA  SET  $2 


TEST  # 

25 

26 

•27 

K  ■ 

8.379 

E'r09 

7.184 

E-09 

6.250 

E-OS 

SK 

5.015 

E-10 

3.601 

E-10 

4.357 

E-10 

RE 

4.265 

E-03 

4.053 

E-03 

3.639 

E-03 

SR 

1.000 

E-04 

1.000 

E-04 

1.000 

E-04 

A 

2.336 

E-04 

2.558 

E-04 

2.469 

E-04 

SA 

I  LAYER 

II  LAYER 

1.124 

4.750 

5.350 

E-05 

1.189 

4.000 

5.675 

E-05 

1.674 

5.250 

6.250 

E-05 

(42-59W) 
SUB  TEST) 


VAPOR  TEMP 

159.000 

156.000 

140.000 

H 

4.300 

4.160 

4.730 

SH 

0.042 

0.042 

0.042 

ALPP 

1.040 

1.040 

1.040 

SRALPP 

0.030 

0.030 

0.030 

ALPPR 

1.329 

1.391 

1.347 

SRALPR 

0.057 

0.101 

0.075 

RL 

1.062 

0.729 

0.916 

SRL 

0.279 

0.332 

0.272 

CALC  WATTS  R 

43.536 

49.895 

45.583 

SWR 

7.900 

8.098 

8  „  3  9  6' 

CALC  WATTS 

53.920 

68.235 

58.079 

SW 

10.754 

8.432 

9.692 

OBS  WATTS 

51.000 

51.000 

57.000 

DT  CALC 

6.433 

7.833 

6.992 

SDT 

1.114 

1.421 

1.461 

OBS  DT 

6.000 

5.000 

6.000 

A9 


(CONTINUED) 


TEST  DATA  SHEET  CONTINUED.  FOR  DATA  SET  #  ? 


(70  -  90W 
SUB  TEST) 


VAPOR. TEMP 

H 

SH 

ALPP 

SRALPP 

ALPPR 

SRALPR 

RL 

SRL 

CALC  WATTS  R 
SWR 

CALC  WATTS 
SW 

OBS  WATTS 
DT  CALC 
SDT 
OBS  DT 


166.000 

186.000 

168.000 

3.970* 

3.760 

4.080 

0.042 

0.042 

0.042 

■1.040 

1.040 

1.040 

0.030 

0.030 

0.030 

1.323 

1.382 

1.396 

0.0646 

0.105 

0.089 

1.083 

0.749 

0.958 

0.2191 

0.330 

0.325 

62.783 

76.951 

76.970 

9.587 

11.413 

11.977 

78.854 

99.996 

96.957 

10.926 

14.798 

15.126 

82.000 

85.000 

91.000 

9.243 

12.118 

12.107 

1.562 

1.861 

2.084 

11.000 

8.000 

8.000 

(100  -  130) W 
SUB  TEST) 


VAPOR  TEMP 

H 

SH 

ALPP 

SRALPP 

ALPPR 

SRALPR 

RL 

SRL 

CALC  WATTS  R 
SWR 

CALC  WATTS 
SW 

OBS  WATTS  ' 
DT  CALC 
SDT 

OBS  DT 


198.000 

192.000 

185  000 

3.100 

3.000 

3.640 

0.042 

0.042 

0.042 

1.040 

1.040 

1.040 

0.030 

0.030 

0.030 

1.342 

1,385 

1.392 

0.078 

0.101 

0.113 

1.166 

0.729 

0.958 

0.464 

0.332 

0.332 

123.787 

126.594 

91.236 

14.974 

16.363 

12.654 

154.466 

167.072 

114.119 

19.230 

24.198 

14.898 

114.000 

116.000 

118.000 

18.670 

20.174 

14.153 

2.480 

3.136 

2.008 

12.000 

11.000 

12.000 

TEST  DATA  SHEET  FOR  DATA  SET  #2 


TEST  # 

29 

30 

31 

K  • 

8.116 

E-09 

6.081 

E-09 

6.326 

E-09 

SK 

4.052 

E-10 

2.641 

E-10 

2.126 

E-10 

RE 

3.961 

E-03 

3,958 

E-03 

3.859 

E-C3 

SR 

1.000 

E-04 

1.000 

E-04 

1.000 

E-04 

A  • 

2.441 

E-04 

2.402 

E-04 

2.359 

E-04 

SA 

1.117 

E-05 

9.200 

E-6 

7.890 

E-06 

I  LAYER 

4.500 

4.750 

4.750 

II  LAYER 

5.750 

5.750 

5.900 

(42-59W) 
SUB  TEST) 


VAPOR  TEMP 

166.000 

155.000 

175.000 

H 

4.460 

4.120 

4.180 

SH 

0.042 

0.042 

0.042 

ALPP 

1.040 

1.040 

1.040 

SRALPP 

0.030 

0.030 

0.030 

ALP  PR 

1.332 

1.329 

1.328 

SRALPR 

0.065 

0.065 

0.065 

RL 

0.895 

0.958 

0.958 

SRL 

0.271 

0.279 

0.289 

CALC  WATTS  R 

46.369 

49.950 

54.853 

SWR 

9.131 

6.951 

7.775 

CALC  WATTS 

59.364 

63.329 

69.723 

SW 

10.683 

7.892 

8.965 

OBS  WATTS 

54.000 

58.000 

53.000 

DT  CALC 

7.011 

7.501 

8.089 

SDT 

1.478 

1.133 

1.226 

OBS  DT 

6.000 

3.000 

10.000 

(CONTINUED) 

1 63 


All 


TEST  DATA  SHEET  CONTINUED  FOR  DATA  SET  #  2 


(70  -  90W 

SUB  TEST) 

VAPOR  TEMP 

192.000 

H 

3.970- 

SH 

0.042 

ALPP 

1.040 

SRALPP 

0.030 

ALPPR 

1.339 

SRALPR 

0.073 

RL 

0.937 

SRI. 

0.310 

CALC  WATTS  R 

77.526 

SWR 

11.672 

•CALC  WATTS 

97.116 

SW 

15.573 

OBS  WATTS 

79.000 

DT  CALC 

12.030 

SDT 

1.859 

OBS  DT 

8.000 

(100  -  130) W 

SUB  TEST) 

VAPOR  TEMP 

200.000 

H 

3.180 

SH 

0.042 

ALPP 

1.040 

SRALPP 

0.030 

ALPPR 

1.352 

SRALPR 

0.069 

RL 

0.917 

SRL 

0.358 

CALC  WATTS  R 

136.161 

SWR 

14.388 

CALC  WATTS 

171.234 

SW 

19.563 

OBS  WATTS 

118.000 

DT  CALC 

20.837 

SDT 

2.201 

OBS  DT 

1.500 

167.000 

194.000 

3.640 

3.540 

0.042 

0.042 

1.040 

1.040 

0.030 

0.030 

1.341 

1.339 

0.075 

0.074 

0.958 

0.979 

0.337 

0.344 

72.873 

88.706 

8.557 

9.884 

91.406 

111.176 

11.531 

13.286 

93.000 

77.000 

11.130 

13.328 

1.351 

1.525 

5.000 

12.000 

198.000 

206.000 

2.770 

2.900 

0.042 

0.042 

1.040 

1.040 

0.030 

0.030 

1.352 

1.348 

0.072 

0.069 

1.102 

0.999 

0.435 

0.398 

122.527 

130.109 

11.518 

11.607 

153.882 

163.779 

16.004 

16.587 

130.000 

114.000 

18.605 

19.395 

1.751 

1.710 

10.000 

13.000 

TEST  DATA  SHEET  FOR  DATA  SET  #2 


TEST  # 

33 

34 

35 

K  • 

8.214 

E-09 

7.324 

E-09 

8.154 

E-09 

SK 

4.302 

E-10 

2.834 

Er-10 

2.877 

E-10 

RE 

4.140 

E-03 

4.337 

E-03 

4.539 

E-03 

SR 

1.000 

E-04 

1.000 

E-04 

1.000 

E-04 

A 

2.426 

E-04 

2.294 

E-04 

2.334 

E-04 

SA 

8.500 

E-06 

8.470 

E-06 

6.270 

E-06 

I  LAYER 

4.625 

4.750. 

4.625 

II  LAYER 

5.500 

5.250 

5.000 

(42-59W) 

SUB  TEST) 

VAPOR  TEMP 

175 

165 

179 

H 

4.210 

3.990 

4.210 

SH 

0.042 

0.042 

0.042 

ALPP 

1.040 

1.040 

1.040 

SRALPP 

0.030 

0.030 

0.030 

ALPPR 

1.318 

1.317 

1.319 

SRALPR 

0.062 

0.068 

0.062 

RL 

1.062 

1.187 

1.021 

SRL 

0.238 

0.259 

0.233 

CALC  WATTS  R 

47.540 

40.836 

51.748 

SWR 

8.899 

6.668 

8.755. 

CALC  WATTS 

60.319 

51.572 

75.837 

SW 

10.298 

7.681 

10,232 

OBS  WATTS 

48.500 

48.000 

50.000 

DT  CALC 

7.203 

5.937 

7.551 

SDT 

1.419 

1.032 

1.331 

OBS  DT 

5.000 

8.000 

8.000 

(CONTINUED) 


A13 


TEST  DATA  SHEET  CONTINUED  FOR  DATA  SET  #  2 


(70  -  90W 

SUB  TEST.) 

VAPOR  TEMP 

190 

182 

180 

H 

3.830' 

3.660 

3 ;  680 

SH 

0.042 

0.042 

0.042 

ALPP 

1.040 

1.040 

1.040 

SRALPP 

0.030 

0.030 

0.030 

ALPPR 

1.324 

1.322 

1.329 

SRALPR 

0.065 

0.064 

0.066 

RL 

1.083 

1.26  8 

0.999 

SRL 

0.263 

0.288 

0.257 

CALC  WATTS  R 

66.836 

58.804 

86.713 

SWR 

9.729 

7.904 

10.211 

CALC  WATTS 

82.935 

72.906 

108.077 

SW 

12.166 

9.899 

12.433 

OBS  WATTS 

80.000 

77.000 

77.000 

DT  CALC 

10.388 

-8.740 

12.910 

SDT 

1.522 

1.191 

1.528 

OBS  DT 

8.000 

8.000 

8.000 

(100  -  130) W 

SUB  TEST) 

VAPOR  TEMP 

202 

194 

189 

H 

3.250 

2.970 

3.18 

SH 

0.042 

0.042 

0.042 

ALPP 

1.040 

1.040 

1.040 

SRALPP 

0.030 

0.030 

0.030 

ALPPR 

1.331 

1.335 

1.333 

SRALPR 

0.058 

0.061 

0.050 

RL 

1.104 

1.208. 

1.042 

SRL 

0.304 

0.359 

0.290 

CALC  WATTS  R 

115.933 

101.329 

123.185 

SWR 

12.139 

9.497 

11.304 

CALC  WATTS 

143.664 

125.615 

153.237 

SW 

14.569 

11.017 

13.306 

OBS  WATTS 

118.000 

109,000 

114.000 

DT  CALC 

17.809 

14.895 

18.203 

SDT 

1.813 

1.388 

1.607 

OBS  DT 

13.000 

10.000 

10.000 

A14 


TEST  DATA  SHEET  FOR  DATA  SET  #2 


TEST  # 

K  - 

SK 

RE 

SR 

A 

SA 

I  LAYER 

II  LAYER 


(42-59W) 
SUB  TEST) 


VAPOR  TEMP 

H 

SH 

ALPP 

SRALPP 

ALPPR 

SRALPR 

RL 

SRL 

CALC  WAT^S 
SWR 

CALC  WATTS 
SW 

OBS  WATTS 
DT  CALC 
SDT 

OBS  DT 


36 

37 

7.475 

E-09 

6.698 

5.334 

E-10 

3.146 

4.141 

E-03 

3.870 

1.000 

E-04 

1.000 

2.421 

E-04 

2.416 

1.300 

E-05 

1.022 

5.000 

5.250 

5.500 

5.875 

38 


E-09 

6.566 

E-09 

Et10 

3.085 

E-10 

E-03 

4.337 

E-03 

E-04 

1.000 

E-04 

E-04 

2.401 

E-04 

E-05 

8.750 

E-06 

4.675 

5.250 


172.000 

165.000 

180.000 

4.220 

4.170 

3.730 

0.042 

0.042 

0.042 

1.040 

1.040 

1.040 

0.030 

0.030 

0.030 

1.316 

1.317 

1.318 

0.067 

0.062 

0.063 

1.208 

1.166 

1.166 

0.251 

0.257 

0.275 

43.123 

61.611 

50.960 

8.549 

8.793 

6.935 

54.369 

77.503 

64.387 

9.742 

9.927 

8.150 

47.000 

54,000 

53.000 

6.593 

9.413 

7.740 

1.418 

1.471 

1.140 

5.000 

7.000 

7.000 

(CONTINUED) 


TEST  DATA  SHEET  CONTINUED  FOR  DATA  SET  #2 


{70  -  90W 
SUB  TEST) 


VAPOR  TEMP 

194.000 

175.000 

195.000 

H 

3.380' 

3.380 

3.460 

SH 

0.042 

0.042 

0.042 

ALPP 

1.040 

1.040 

1.040 

SRALPP 

0.030 

0.030 

0.030 

ALPPR 

1.327 

1.327 

1.324 

SRALPR 

0.067 

0.067 

0.065 

RL 

1.208 

1.187 

1.208 

SRL 

0.332 

0.332 

0.305 

CALC  WATTS  R 

95.714 

101.988 

65.902 

SWR 

12.006 

11,386 

7.876 

CALC  WATTS 

119.515 

127.300 

82.085 

SW 

15.745 

13.903 

10,315 

OBS  WATTS 

75.000 

77.000 

76.000 

DT  CALC 

14.956 

15.862 

10.206 

SDT 

2.086 

1.850 

1.249 

OBS  DT 

10.000 

9,000 

12.000 

(100  -  130) W 

SUB  TEST) 

VAPOR  TEMP 

195.000 

205.000 

205.000 

H 

2.940 

2.650 

2.86 

SH 

0.042 

0.042 

0.042 

ALPP 

1.040 

1.040 

1 .040 

SRALPP 

0.030 

0,030 

0.030 

ALPPR 

1.338 

1.346 

1.336 

SRALPR 

0.063 

0.068 

0.061 

RL 

1.208 

1.208 

1.208 

SRL 

0.392 

0.443 

0.373 

CALC  WATTS  R 

126.022 

152.780 

102.459 

SWR 

14.558 

14.981 

8.918 

CALC  WATTS 

155.837 

189.397 

127.379 

SW 

17.116 

15.879 

11.376 

OBS  WATTS 

110.000 

119.000 

118.000 

DT  CALC 

19.456 

23.582 

15.716 

SDT 

2.273 

2.344 

1.377 

OBS  DT 

10.000 

14.000 

18.000 

t 

A16 


TEST  DATA  SHEET  FOR  DATA  SET  #2 


TEST  # 

K  • 

39 

7^192 

E-09 

'  40 

5.951 

E-09 

41 

6.235 

Et09 

SK 

3.193 

E-10 

3.363 

E-10 

2.832 

E-10 

RE 

4.051 

E-03 

3.878 

E-03 

3.877 

E-03 

SR 

1.000 

E-04 

1.000 

E-04 

1.000 

E-04 

A 

2.392 

E-04 

2.416 

E-04 

2.363 

E-04 

SA 

6.880 

E-06 

1.134 

E-05 

7.780 

E-06 

I  LAYER 

II  LAYER 

5.000 

5.625 

5.250 

5.875 

4.750 

5.875 

(42-59W) 

SUB  TEST) 

VAPOR  TEMP 

175.000 

160.000 

174.000 

H 

4.010 

4.560 

4.310 

SH 

0.042 

0.042 

0.042 

ALPP 

1.040 

1.040 

1.040 

SRALPP 

0.030 

0.030 

0.030 

ALPPR 

1.317 

1.315 

1.322 

SRALPR 

0.063 

0.062 

0.062 

RL 

1.187 

1.146 

0.979 

SRL 

0.260 

0.228 

0.225 

CALC  WATTS  R 

58.780 

34.876 

46.413 

SWR 

8.036 

6.665 

7.214/ 

CALC  WATTS 

74.435 

44.341 

59.278 

SW 

9.54  8 

7.920 

8.667 

OBS  WATTS 

53.000 

50,000 

51.000 

DT  CALC 

8.892 

5.299 

6.853 

SDT 

1.287 

1.089 

1.130 

OBS  DT 

8.000 

8.000 

8.000 

(CONTINUED) 


TEST  DATA  SHEET  CONTINUED  FOR  DATA  SET  #2 


i 


b 

- 


{70  -  90W 
SUB  TEST) 


VAPOR  TEMP 

187.000 

175.000 

180.000 

H 

4.010’ 

3.900 

3.630 

SH 

0.042 

0.042 

0.0*42 

ALPP 

i.040 

1.040 

1.040 

SRALPP 

0.030 

0.030 

0.030 

ALPPR 

1.328 

1.321 

1.332 

SRALPR 

0.067 

0.064 

0.067  . 

RL 

1.208 

1.187 

0.958 

SRL 

0.341 

0.278 

0.253 

CALC  WATTS  R 

108.078 

65.546 

80.464 

SWR 

10.049 

8.549 

8.619 

CALC  WATTS 

134.972 

81.772 

101.040 

SW 

13.158 

11.076 

11.275 

OBS  WATTS 

84.000 

73.000 

78,000 

DT  CALC 

16.639 

10.205 

12.094 

SDT 

1.597 

1.384 

1.327 

OBS  DT 

10.000 

10.000 

10.000 

(10G  -  130) W 
SUB  TEST) 


VAPOR  TEMP 

195.000 

189.000 

190.000 

H 

2.770 

3.160 

3.120 

SH 

0.042 

0.042 

0.042 

ALPP 

1.040 

1.040 

1.040 

SRALPP 

0.030 

0.030 

0.030 

ALPPR 

1.340 

1.334 

1.337 

SRALPR 

0.063 

0.060 

0.059 

RL 

1.208 

1.208 

0.950 

SRL 

0.402 

0.357 

0.273 

CALC  WATTS  R 

y 

140.222 

104.590 

108.881 

SWR 

10.791 

10.281 

9.333 

CALC  WATTS 

174.624 

129.736 

136.060 

SW 

13.838 

12.980 

11.979 

OBS  WATTS 

115.000 

120.000 

125.000 

DT  CALC 

21.447 

16.122 

16.206 

SDT 

1.665 

1.616 

1.366 

OBS  DT 

14.000 

18.000 

12.000 

ito 


A18 


TEST  DATA  SHEET  FOR  DATA  SET  #  2 


TEST  # 

42 

43 

44 

K  - 

6.954 

E-09 

7.483 

E-09 

-  6.081 

E-09 

SK 

3.122 

E-IQ 

4. -465. 

E~i  0 

2.298 

E-10 

RE 

3.969 

E-03 

4.140 

E-03 

3.969 

E-03 

SR 

1.000 

E-04 

1.000 

E-04 

1.000 

E-04 

A 

2.426 

E-04 

2.382 

Et-04 

2.353 

E-04 

SA 

1.169 

E-05 

1.236 

E-06 

7.960 

E-06 

I  LAYER 

.  5.365 

4.750 

5.000 

II  LAYER 

5.750 

5.500 

5.750 

(42-59W) 

■ 

SUB  TEST) 

VAPOR  TEMP 

165.000 

185.000 

176.000 

H 

4.440 

4.210 

4.170 

SH 

0.042 

0.042 

01042 

ALPP 

1.040 

1.040 

1.040 

SRALPP 

0.030 

0.030 

0.030 

ALPPR 

1.318 

1.323 

1.322 

SRALPR 

0.062 

0.059 

0.059 

RL 

1.041 

0.895 

0.894 

SRL 

0.226 

0.229 

0.224 

CALC  WATTS  R 

40.253 

40.468 

45.587 

SWR 

7.542 

8.214 

6.703' 

CALC  WATTS 

51.350 

51.983 

58.509 

SW 

9.047 

9.800 

8.129 

OBS  WATTS 

53.000 

45,000 

45.000 

DT  CALC 

6.107 

5.986 

6.683 

SDT 

1.236 

1.313 

1.049 

OBS  DT 

10.00,0 

8.000 

5.000 

(CONTINUED) 

!&(  . 

A19 


TEST  DATA  SHEET  CONTINUED  FOR  DATA  SET  #2 


(70  -  90W 

SUB  TEST) 

VAPOR  TEMP 

185.000 

192,000 

198.000 

H 

3.810' 

3.490 

3.420 

SH 

0.042 

0.042 

0.042 

ALPP 

1.040 

1.040 

1.040 

SRALPP 

0.030 

0,030 

0.030 

ALPPR 

1.327 

1.343 

1.341 

SRALPR 

0.065 

0.077 

0.760 

RL 

1.020 

0.874 

0.895 

SRI, 

0.256 

0.256 

0.255 

CALC  WATTS  R 

75.292 

85.704 

83.753 

SWR 

9.646 

10.709 

8.474 

CALC  WATTS 

94.157 

108.543 

105.741 

SW 

12.601 

14.156 

11.220 

OBS  WATTS 

86.000 

88.000 

72.000 

DT  CALC 

11.656 

13.420 

12.483 

SDT 

1.561 

1.702 

1.305 

OBS  DT 

8.000 

10.000 

8.000 

(100  -  130) W 

SUB  TEST) 

VAPOR  TEMP 

190.000 

205.000 

198.000 

H 

3.140 

3.000 

3.070 

SH 

0.042 

0.042 

0.040 

ALPP 

1.040 

1.040 

1.040 

SRALPP 

0.030 

0.030 

0.030 

ALPPR 

1.319 

1.333 

1.329 

SRALPR 

0.061 

0.064 

0.063 

RL 

1.020 

0.854 

0.874 

SRL 

0.285 

0.264 

0.254 

CALC  WATTS  R 

116.538 

118.852 

102.627 

SWR 

10.946 

12.124 

8.546 

CALC  WATTS 

145.170 

149.571 

129.286 

SW 

13.833 

15.507 

11.046 

OBS  WATTS 

114.000 

121.000 

120.000 

DT  CALC 

17.811 

17.700 

15.137 

SDT 

1.715 

1.846 

1.245 

OBS  DT 

8.1)00 

'  13; 000 

'  10 . 000 

TEST  DATA  SHEET  FOR  DATA  SET  #2 


TEST  # 

45 

46 

47 

K  • 

8.702  E-09 

8.338  E-09 

6.427  E-09 

SK 

3.439  E-10 

2.742  E-10 

3.796  E-10 

RE 

4.143  E-03 

4.140  E-03 

4.051  E-03. 

SR 

1.000  E-04 

1.000  E-04 

1.000  E-04 

A 

2.337  E-; 04 

2.270  E-04 

2.387  E-04 

SA 

6.280 .E-06 

4.900  E-06 

1.001  E-05 

I  LAYER 

4.875 

4.875 

4.050 

II  LAYER 

5.500 

5.500 

5.625 

(42-59W) 

SUB  TEST) 

VAPOR  TEMP 

168.000 

170.000 

179.000 

H 

4.078 

4.280 

4.140 

SH 

0.042 

0.042 

0.042 

ALPP 

1.040 

1.040 

1.040 

SRALPP 

0.036 

0.030 

0.030 

ALPPR 

1.317 

1.325 

1.346 

SRALPR 

0.058 

0.060 

0.064 

RL 

0.937 

0.854 

0.708 

SRL 

0.232 

0.221 

0.215 

CALC  WATTS  R 

65.462 

40.704 

44.146 

SWR 

9.656 

7.538 

7.281 

CALC  WATTS 

83.766 

52.739 

57.578 

SW 

11.397 

9.220 

8.816 

OBS  WATTS 

51.000 

50.000 

52.000 

DT  CALC 

9.530 

5.725 

6.495 

SDT 

1.446 

1.092 

1.157 

OBS  DT 

6.000 

.  10.000 

6.000 

(CONTINUED) 


!(&■ 


•3 


A21 


TEST  DATA  SHEET  CONTINUED  FOR  DATA  SET  #  2 


(70  -  90W 
SUB  TEST) 


VAPOR  TEMP 

185.000 

200.000 

190.000 

H 

3.460' 

3.58G 

3.480 

SH 

0.042 

0.042 

0.042 

ALPP 

1.040 

1.040 

1.040 

SRALPP 

0.030 

0.030 

0.030 

ALPPR 

1.337 

1.347 

1.367 

SRALPR 

0.075 

0.077 

0.082 

RL 

0.937 

0.833 

0.687 

SRL 

0.261 

0.244 

0.228 

CALC  WATTS  R 

99.720 

84.430 

78.652 

SWR 

10.273 

9.422 

9.171 

CALC  WATTS 

125.386 

106.740 

101.497 

SW 

13.540 

12,510 

12.328 

OBS  WATTS 

81.000 

83.000 

83.000 

DT  CALC 

14.706 

12.143 

11.786 

SDT 

1.542 

1.354 

1.433 

OBS  DT 

8.000 

10.000 

10.000 

(100  -  130) W 
SUB  TEST) 


VAPOR  TEMP 

190.000 

205.000 

190.000 

H 

2.970 

2.680 

2.910 

SH 

0.042 

0.042 

0.042 

ALPP 

1.040 

1.040 

1.040 

SRALPP 

0.030 

0,030 

0.030 

ALPPR 

1.326 

1.337 

1.356 

SRALPR 

0.067 

0.066 

0.0679 

RL 

0.916 

0.812 

0.666 

SRL 

0.270 

0.260 

0.227 

CALC  WATTS  R 

135.919 

149.236 

109.447 

SWR 

10.797 

10.381 

10.220 

CALC  WATTS 

170.892 

189.573 

140.463 

SW 

13.992 

13.625 

13.463 

OBS  WATTS 

112.000 

130.000 

125.000 

DT  CALC 

19.941 

21.276 

16.187 

SDT 

1,538 

1.423 

1.512 

OBS  DT  . 

12.000 

13.001 

14.000 

/(s?f 


A22 


TEST  DATA  SHEET  CONTINUED  FOR  DATA  SET  #  2 


(70  -  90W 
SUB  TEST) 


VAPOR  TEMP 

189.000 

180.000 

190.000 

H 

3.580* 

2.950 

3.310 

SH 

0.042 

0.042 

0.042 

ALPP 

1.040 

1.040 

1.040 

SRALPP 

0.030 

0.030 

0.030 

ALPPR 

1.349 

1.352 

1.348 

SRALPP. 

0.077 

0.080 

0.078 

RL 

0.812 

0.812 

0.833 

SRL 

0.233 

0.260 

0.254 

CALC  WATTS  R 

82.475 

113.993 

71.195 

SWR 

8.437 

9.536 

6,899 

CALC  WATTS 

104.588 

145.198 

90.239 

SW 

11.175 

12.906 

9.238 

OBS  WATTS 

80.000 

89.000 

80.000 

DT  CALC 

12.401 

16.294 

10.092 

SDT 

1.282 

1  435 

0.998 

OBS  DT 

10.000 

10.000 

10.000 

(100  -  130) W 
SUB  TEST) 


VAPOR  TEMP 

192,000 

192.000 

190.000 

H 

3.040 

2.700 

2.540 

SH 

0.042 

0.042 

0.042 

ALPP 

1.040 

1.040 

1.040 

SRALPP 

0.030 

0.030 

0.030 

ALPPR 

1.337 

1.337 

1.339 

SRALPR 

0.064 

0.066 

0.066 

RL 

0.791 

0.812 

0.792 

SRL 

0.236 

0.262 

C.261 

CALC  WATTS  R 

111.588 

133.755 

110.377 

SWR 

8.870 

9.888 

7.484 

CALC  WATTS 

141.818 

169.661 

140.466 

SW 

11.636 

12.848 

9.801 

OBS  WATTS 

121.000 

110.000  - 

110.000 

DT  CALC 

16.592 

18.949 

15.487 

SDT  ' 

1.270 

1.402 

1.040 

OBS  DT  • 

12.000 

10.000 

15.000 

!b$' 


A23 


TEST  DATA  SHEET 

FOR  DATA  SET  #2 

TEST  # 

48 

49 

50 

K  • 

€.204 

E-09 

8.246 

E-09 

6.412 

E-09 

SK 

2.576 

E-10 

2.674 

E-10 

1.968 

E-10 

RE 

3.887 

E-03 

4.257 

E-03 

4.257 

E-03 

SR 

1.000 

E-04 

1.000 

E-04 

1.010 

E-04 

A 

2.387 

E-04 

2.261 

E-04 

2.236 

E-04 

SA 

6.000 

E-06 

7.510 

E-06 

6.290 

E-06 

I  LAYER 

4.875 

4.500 

4.500 

II  LAYER 

5.875 

5.365 

5.365 

(42-59W) 

SUB  TEST) 

VAPOR  TEMP 

176.000 

174.000 

176. 00G 

H 

4.210 

3.780 

3.870 

SH 

0.042 

0.042 

0.042 

ALPP 

1.040 

1.040 

1.040 

SRALPP 

0.030 

0.030 

0.030 

ALPPR 

1.331 

1.330 

1.330 

SRALPR 

0.061 

0.062 

0.062 

RL 

0.812 

0.833 

0.830 

SRL 

0.216 

0.235 

0.233 

CALC  WATTS  R 

54.535 

6’  .828 

43.227 

SWR 

7.240 

7.782 

5.839 

CALC  WATTS 

70.427 

79.545 

55.724 

SW 

8.826 

9.459 

7.127 

OBS  WATTS 

49.000 

58.000 

43.000 

DT  CALC 

8.067 

8,728 

6.031 

SDT 

1.127 

1.188 

0.862 

OBS  DT 

6.000 

10.000 

7.000 

(CONTINUED) 

TEST  DATA  SHEET  FOR  DATA  SET  #2 


TEST  # 

51 

52 

53 

1 

K  ■ 

8.337 

E-09 

7.463 

E-09 

7.775 

E-09 

SK 

3.730 

E-10 

1.382 

E-10 

3.343 

E-10 

RE 

4.140 

E-03 

4.140 

E-03 

4.245 

E-03 

SR 

1.000 

E-04 

1.000 

E-04 

1.000 

E-04 

A 

2.373 

E-04 

2.324 

E-04 

2.392 

E-04 

SA 

6.270 

E-06 

5.320 

E-06 

9.780 

E-06 

I  LAYER 

5.000 

4.750 

4.500 

II  LAYER 

5.500 

5.500 

5.365 

(42-59W) 
SUB  TEST) 


VAPOR  TEMP 

168.000 

184.000 

185.000 

H 

4.160 

4.020 

3.980 

SH 

0.042 

0.042 

0.042 

ALPP 

1.040 

1.040 

1.040 

SRALPP 

0.030 

0.030 

0.030 

ALP  PR 

1.315 

1.323 

1.330 

SRALPR 

0.057 

0.059 

0.061 

RL 

0.999 

0.896 

0.833 

SRL 

0.235 

0.232 

0.233 

CALC  WATTS  R 

50.830 

50.886 

49.416 

SWR 

8.118 

7.406 

8.126 

CALC  WATTS 

64.889 

65.495 

63.716 

SW 

9.722 

9.112 

9.813 

OBS  WATTS 

48.000 

54.000 

49.000 

1 

| 

DT  CALC 

7.529 

7.362 

7.378 

t 

SDT 

1.255 

1.116 

1.297 

OBS  DT 

8.000 

6.000 

5.000 

(CONTINUED) 

.  til 


A25 


TEST  DATA  SHEET  CONTINUED  FOR  DATA  SET  #2 


(70  -  90W 

SUB  TEST) 

VAPOR  TEMP 

177.000 

192.000 

194.000 

H 

3.610* 

3.560 

3.660 

SH 

0.042 

0.042 

0.042 

ALPP 

1.040 

1.040 

1.040 

SRALPP 

0.030 

0.030 

0.030 

ALPPR 

1.330 

1.343 

1.347 

SRALPR 

0.073 

0.076 

0.077 

RL 

0.979 

0.874 

0.833 

SRL 

0.260 

0..  247 

0.245 

CALC  WATTS  R 

86.563 

79.212 

69.306 

SWR 

9.649 

8.123 

9.610 

CALC  WATTS 

108.494 

99.799 

87.654 

SW 

12.643 

10.782 

12.113 

OBS  WATTS 

86.000 

82.000 

85.000 

DT  CALC 

13.060 

11.655 

10.487 

SDT 

1.466 

1.202 

1.416 

OBS  DT 

10.000 

6.000 

8.000 

(100  -  130) W 

SUB  TEST) 

VAPOR  TEMP 

189.000 

202.000 

202.000 

H 

2.860 

2.990 

2.990 

SH 

0.042 

0.042 

0.042, 

ALPP 

1.040 

1.040 

1.040 

SRALPP 

0.030 

0.030 

0.030 

ALPPR 

1.324 

1.332 

1.337 

SRALPR 

0.063 

0.064 

0.065 

RL 

0.950 

0.854 

0.812 

SRL 

0.287 

0.256 

0.257 

CALC  WATTS  R 

140.225 

116.434 

116.188 

SWR 

11,067 

8.609 

10.409 

CALC  WATTS 

175.901 

147.187 

146.990 

SW 

14.331 

11.205 

13.472 

OBS  WATTS 

18.000 

125.000 

125.000 

DT  CALC 

20.936 

16.966 

17.341 

SDT- 

1.606 

1.209 

1.559 

OBS  DT  . 

10.000 

10.000 

10.000 

U>4 


A26 


TEST  DATA  SHEET  FOR  DATA  SET  #3 


TEST  # 

K  - 

SK 

RE 

SR 

A 

SA 

I  LAYER 

II  LAYER 


54 

3.074  E-09 
2.845  E-10 
2.336  E-03 
1.000  E-04 
1.614  E-04 
1.294  E-05 
8.000 
9.750 


55 

2.584  E-09 
1.483  E-10 
2.396  E-03 
1.000  E-04 
1.579  E-04 
5.390  E-06 
7.000 
9.500 


56 


3.796  E-09 
6.949  E-10 
2.525  E-03 
1.000  E-04 
1.688  E-04 
8.570  E-06 


(  20-30W) 
SUB  TEST) 


VAPOR  TEMP 

H 

SH 

ALPP 

SRALPP 

ALPPR 

SRALPR 

RI, 

SRL 

CALC  WATTS  R 
SWR 

CALC  WATTS 
SW 

OBS  WATTS 
DT  CALC 
SDT 

OBS  DT 


152.000 

146.000 

164.000 

7.230 

7.100 

7.010 

0.042 

0.042 

0.042 

1.040 

1.040 

1.046 

0.030 

0.030 

0.030 

1.227 

1.239 

1.225 

0.055 

0.056 

0.055 

0.937 

0.729 

0.999 

0,180 

0.148 

0.193 

26.665 

20.406 

25.080 

6.529 

4.382 

8.379 

32.035 

24.953 

30.116 

7.399 

5.118 

9.191 

26.500 

25.000 

25.400 

2.878 

2.136 

2.779 

0.779 

0.479 

0.976 

2.000 

1.000 

2.000 

(CONTINUED) 


A27 


TEST  DATA  SHEET  CONTINUED  FOR  DATA  SET  #  3 


(70  -  90W 
SUB  TEST) 


VAPOR  TEMP 

H 

SH 

ALPP 

SRALPP 

ALPPR 

SRALPR 

RL 

SBXj 

CALC  WATTS  R 
SWR 

CALC  WATTS 
SW 

OBS  WATTS 
DT  CALC 
SDT 
OBS  DT 


150.000 

2.860 

0.042 

1.040 

0.030 

1.257 

0.075 

0.854 

0.377 

92.933 

13.026 

110.930 

15.865 

78.000 

10.231 

1.588 

4.000 


192.000 

3.620 

0.042 

1.040 

0.030 

1.255 

0.069 

0.749 

0.216 

74.076 

7.632 

88.276 

9.227 

70.000 

7.952 

0.839 

3.000 


200.000 

3.980 

0.047 

1.040 

0.030 

1.255 

0.074 

1.083 

0.498 

99.119 

20.465 

117.492 

23.810 

75.000 

11.474 

2.384 

4.000 


no 


TEST  DATA  SHEET  FOR  DATA  SET 


TEST  # 

K  ■ 

SK 

RE 

SR 

A 

SA 

X  LAYER 

XI  LAYER 


57 

2.798  E-09 
9.819  E-ll 
2.162  E-03 
1.000  E-04 
1.677  E-04 
6.440  E-06 
9.250 
10.250 


58 

3.171  E-09 
2.146  L-10 
2.603  E-03 
1.000  E-04 
1.682  E-04 
7.860  E-06 
7.350 
8.750 


(20-30W) 
SUB  TEST) 


VAPOR  TEMP 

H 

SH 

ALPP 

SRALPP 

ALP  PR 

SRALPR 

RL 

SRL 

CALC  WATTS  R 
SWR 

CALC  WATTS 
SW 

OBS  WftTTS 
DT  CALC 
SDT 

OBS  DT 


162.000 

7.360 

0.042 

1.040 

0.030 

1.222 

0.055 

1.100 

0.181 

34.337 

6.742 

41.373 

7.892 

29.000 

3.842 

0.797 

3.000 


175.000 

6.850 

0.042 

1.040 

0.030 

1.222 

0.055 

1.062 

0.070 

18.559 

5.653 

22.590 

6.496 

24.000 

2.060 

0.659 

2.000 


(CONTINUED) 


m 


TEST  DATA  SHEET  CONTINUED  FOR  DATA  SET  #3 


(60  -  80W) 

SUB  TEST) 

VAPOR  TEMP 

203.000 

194.000 

H 

4.350 

4.050 

SH 

0.042 

0.042 

ALPP 

1.040 

1.046 

SRALPP 

0.030 

0.030 

ALPPR 

1.244 

1.242 

SRALPR 

0.067 

0.067 

RL 

1.208 

1.083 

SRL 

0.376 

0.314 

CALC  WATTS  R 

92.092 

75.657 

SWR 

9.862 

9.050 

CALC  WATTS 

107.954 

89.044 

SW 

11.745 

10.827 

OBS  WATTS 

81.000 

79.000 

DT  CALC 

10.581 

8.688 

SDT 

1.168 

1.082 

OBS  DT 

5.000 

8.000 

/7£ 


A30 


TEST  DATA  SHEET  FOR  DATA  SET  #4 


TEST  # 

K  ■ 

SK 

RE 

SR 

A 

SA 

I  LAYER 

II  LAYER 


2.729  E-0 8 
1.339  E-09 
7.921  E-03 
1.000  E-04 
5.314  E-04 
1.890  E-05 
2.365 
2.875 


3.209  E-08 
1.381  E-09 
8.278  E-03 
1.000  E-04 
5.203  E-04 
2.091  E-06 
2.365 
2.750 


3.291  E-08 
1.265  E-09 
8.675  E-03 
1.000  E-04 
5.406  E-04 
1.860  E-06 
2.000 
2.625 


(  50-60W) 
SUB  TEST) 


VAPOR  TEMP 

H 

SH 

ALPP 

SRALPP 

ALPPR 

SRALPR 

RL 

SRL 

CALC  WATTS  R 
SWR 

CALC  WATTS 
SW 

OBS  WATTS 
DT  CALC 
SDT 

OBS  DT 


182.000 

2.480 

0.042 

1.040 

0.030 

1.328 

0.064 

1.041 

0.403 

35.003 

26.410 

44.122 

30.879 

55.000 

-10.778 

8.668 

8.000 


195.000 
2.290 
0.  042 
1.040 
0.030 
1.328 
0.063 
1.166 
0,439 
68.111 
32.289 
85.004 
37.535 
56.000 
21.366 
10.562 
8.000 


188.000 

2.250 

0.042 

1.040 

0.030 

1.342 

0.076 

0.916 

0.442 

38.561 

31.797 

48.406 

37.576 

56.000 

11.904 

10.533 

8.001 


(CONTINUED) 

//3  . ___ 


A31 


TEST  DATA  SHEET  FOR  DATA  SET  #4 


(100-130  W) 
SUB  TEST) 


VAPOR  TEMP 

H 

SH 

ALPP 

SRALPP 

ALPPR 

SRALPR 

RL 

SRL 

CALC  WATTS  R 
SWR 

CALC  WATTS 
SW 

OBS  WATTS 
DT  CALC 
SDT 

OBS  DT 


200.000 

2.250 

0.042 

1.040 

0.030 

1.348 

0.085 

1.083 

0.440 

133.845 

29.706 

163.944 

40.669 

115.000 

45.000 

10.006 

18.000 


200.000 

2.230 

0.042 

1.040 

0.030 

1.345 

0.082 

1.117 

0.449 

286.371 

32.232 

103.088 

44.143 

120.000 

28.965 

10.608 

14.000 


200.000 

2.1D0 

0.042 

1.040 

0.030 

1.364 

0.098 

0.937 

0.464 

118.182 

33.356 

144.277 

45.878 

123.000 

40.421 

11.280 

14.000 


TEST  DATA  SHEET  FOR  DATA  SET  #4 


TEST  # 

K  • 

SK 

RE 

SR 

A 

SA 

I  LAYER 

II  LAYER 


62 

2.822  E-08 
8.494  E-10 
8.675  E-03 
1.000  E-04 
5.438  E-04 
1.001  E-05 
1.875 
2.625 


63 

3.231  E-08 
9.071  E-10 
8.675  E-03 
1.006  E-04 
5.295  E-04 
1.100  E-05 
2.250 
2.625 


( 50-60W) 
SUB  TEST) 


VAPOR  TEMP  176 
H  2. 
SH  0. 
ALPP  1 • 
SRALPP  0 . 
ALPPR  1 . 
SRALPR  0. 
RL  0. 
SRL  0. 
CALC  WATTS  R  23. 
SWR  23. 
CALC  WATTS  30. 
SW  30  • 
OBS  WATTS  53  • 
DT  CALC  *7  • 
SDT  8* 
OBS  DT  8  * 


.000 

170.000 

300 

2.240 

042 

0.042 

040 

1.040 

030 

0.030 

355 

1.331 

084 

0.065 

791 

1.124 

432 

0.448 

412 

158.707 

637 

29.374 

089 

73.392 

476 

34.452 

000 

55.000 

027 

18.544 

365 

9.631 

000 

8,000 

(CONTINUED) 


TEST  DATA  SHEET  FOR  DATA  SET  #5 


TEST  # 

64 

65 

66 

K  - 

6.737 

E-09 

6. 031 

E-09 

■6.421 

E-09 

SK 

2.587 

E-10 

4.131 

E-10 

4.344 

E-10 

RE 

3.888 

E-03 

3.888 

E-03 

3.888 

E-03 

SR 

1.000 

E-04 

1.000 

E-04 

1.000 

E-04 

A 

3.092 

E-04 

2.759 

E-04 

2.939 

E-04 

SA 

9.650 

E-06 

1.719 

E-05 

1.053 

E-05 

I  LAYER 

5.750 

5.250 

5,250 

II  LAYER 

5.875 

5.875 

5.875 

(  45-55W) 

SUB  TEST) 

VAPOR  TEMP 

175.000 

201.000 

185.000 

H 

4.310 

4.580 

4.490 

SH 

0.042 

0.042 

0.042 

ALPP 

1.040 

1.040 

1.040 

SRALPP 

0.030 

0.030 

0.030 

ALPPR 

1.308 

1.307 

1.305 

SRALPR 

0.055 

0.055 

0.055 

RL 

1.208 

1.208 

1.145 

SRL 

0.245 

0.230 

0.735 

CALC  WATTS  R 

61.077 

30.988 

45.768 

SWR 

9.506 

8.860 

10.015 

CALC  WATTS 

77.575 

39.683 

58.129 

SW 

11.395 

10.535 

11.815 

OBS  WATTS 

53.000 

45.000 

55.000 

DT  CALC 

9.561 

4.284 

6.803 

SDT 

1.587 

1.320 

1.650 

OBS  DT 

7.000 

5.000 

8.000 

Ml 


A3  5 


(CONTINUED) 


TEST  DATA  SHEET  CONTINUED  FOR  DATA  SET  #5 


(90-120W 

SUB  TEST) 

VAPOR  TEMP 

185.000 

H 

3.550 

SH 

0,042 

ALPP 

1.040 

SRALPP 

0.030 

ALPPR 

1.327 

SRALPR 

.  0.068 

RL 

1.208 

SRI/ 

0.308 

CALC  WATTS  R 

112.410 

SWR 

11.538 

CALC  WATTS 

140.247 

SW 

15.031 

OBS  WATTS 

104.000 

DT  CALC 

17.930 

SDT 

1.911 

OBS  DT 

15.000 

201.000 

180.000 

3.590 

3.750 

0.042 

0.042 

1.040 

1.040 

0.030 

0.030 

1.328 

1.327 

0.068 

0.068 

1.208 

1.200 

0.315 

0.302 

94.374 

87.366 

12.657 

11.942 

118.145 

109.325 

16.431 

15.490 

105.000 

110.000 

13.421 

13.764 

1.985 

1.997 

15.000 

19.000 

/  7 1 


TEST  DATA  SHEET  FOR  DATA  SET  #5 


TEST  # 

K  - 

SK 

RE 

SR 

A 

SA 

I  LAYER 

II  LAYER 


(45-55W 
SUB  TEST) 


VAPOR  TEMP 

H 

SH 

ALPP 

SRALPP 

ALPPR 

SRALPR 

RL 

SRL 

CALC  WATTS  R 
SWR 

CALC  WATTS 
SW 

OBS  WATTS 
DT  CALC 
SDT 

OBS  DT 


67  68 


6.400 

E-09 

■6.577  E-09 

3.475 

E-10 

3.111‘  E-10 

3.971 

E-03 

3.888  E-03 

1.000 

E-04 

1.000  E-04 

2.836 

E-04  ‘ 

2.893  E-04 

1.474 

E-05 

1.037  E-05 

5.165 

5.250 

5.750 

5.875 

190.000 

3.950 

0.042 

1.040 

0.030 

1.308 

0.055 

1.124 

0.263 

71.543 

10.360 

90.255 

12.109 

55.000 

10.302 

1.710 

8.000 


186.000 

4.380 

0.042 

1.040 

0.030 

1.308 

0.055 

1.208 

0.240 

52.877 

9.568 

67.271 

11.431 

55.000 

7.715 

1.487 

9.000 


(CONTINUED) 


m 


TEST  DATA  SHEET  CONTINUED  FOR  DATA  SET  t#5 


(90  -  120  W 
SUB  TEST) 


VAPOR  TEMP 
H 

SH 

ALPP 

SRALPP 

ALPPR 

SRALPR 

RL 

SRL 

CALC  WATTS  R 
SWR 

CALC  WATTS 
SW 

OBS  WATTS 
DT  CALC 
SDT 

or;s  dt 


no 


A3  8 


195.000 

186.000 

3.650 

3.490 

0.042 

0. 0*42 

1.040 

1.040 

0.030 

0.030 

1.327 

1.328 

0.068 

0.067 

1.208 

1.020 

0.301 

0,279 

91.632 

113.258 

11.601 

11.854 

114.474 

141.463 

15.102 

15.440 

95.000 

100.000 

13.395 

16.866 

1.822 

1.857 

15.000 

18.000 

APPENDIX  B 

COMPUTER  PROGRAM 


68.  05 


/ID  0 16367080, MEGS9  ROD ER T S , CHARLES 

/JOB  TIME=10,G0 

’  DIMENSION  AANI400)  , F ( 400 > , OBF ( 400 ) , D  IFF (4 00 > . X ( 400 ) 
COMMON  RHOF.UL . AL , H IL ♦ RE . SR . HFG . AK . SK » SAL 
DO  681  11=1,10 
READ, MV 
READ, SR V 
PV=RV/2. 

SRV=SHV/2, 

READ ,H1 L 
F1L=HIL/12. 

RE  AC  »  AL 


RE AD, SAL 
RC AD, AK 
READ , SK 
READ , Hr 
READ, SR 

READ  100, RHOF , HFG 
READ  100, SIGMA, UL 


READ  1 0 0  ,  H , 3H 
READ  1  00 ,TT ,ALPP 

100  FORM4T(E20.7.E20.7) 

200  FORVAT(  ;£20. 7, FP.O, 7,820.7*  1 
WR  I  TE( 6, 1  02) 

1 02  FORMAT (GX ♦ • AK* , 1 5X , ' SK • , 1 5X , • RE* , 1 5X , • SR • , 1 5X» ' H • , 15X, • SH • , 
WR I TF  C  6  t 101  )  AK ,SK  »  RR , SR i H  *SH 

101  FORMAT  (  2X ,  1PF.1  ?.£,  4X,  1PE12.5.4X,  1PE  1  2.  5 . 4X  ,  1PE  1 2.  5 , 4X  ♦  1  PE  I  ?. 
14X, IPF12.S,/) 

SRALPR=.03 


00  400  KL=  1,3 
WP.I  TE(  6 , 1  03  ) 

103  FORMAT!  6X»  *  ALPP  *  ,/) 
WR I TE (6,104)  ALPP 

104  FORVATC 2X, IPE12.5./) 
PR  INT.SRALPP 
WRITE<6 , 165) 

165  FORMAT!  6X,  '  TEMP’  ,/) 
WRITE (6, 166)  TT 
FORMAT! 2X» 1  PE  12. 5,/) 


166 


' v'VTT 

> a  ~  -  r-r ‘  s;  .-  -“v  <S; ^ 

"  <  :;  w  *V  -'V\rA  c?»i 


'  3=.5'?.iiq?'  '  ■  "  -  _■  *  /..-  :J;V'VV;^--:.V,;' 

GG=’J2^i’^'4.  -  -•  •  _  _  '  -V*  \  :-.  -  .;-;>?.r 

'  ^^ll^O^S'T^^At'K'.rG^UL  ■  •  ••'.•-  '  ;  .  ■’  V-;--'  ■  '. 

.G=.'sr'*f>j»'**V7. i.n.7  /  -  .  -;  }  '  \ 

'  rfs-ctff  -  _  '  . :  -  y :- 

Hl^«A,f:(6X?  • Aftf £■•'»!•  £>•  •■AWL*  ./>  ’  ' 

.  WTTF  i 0  i  i  56  )  Al.  i  AM:  ‘.  ,  ’  ' 

f*BH  VAT>2X  ,  1‘Pr.l  2i.S‘.  A-X  ,  1 FC 1 2 i 3 ./> 

prUNTiSAl-. 

.Vw  l-Tjify  .201  )  .  j  .  ' 

"0  'MAt(  2X«  «sr/  T  r‘l.f.y,  WATTS*  .  2X-.  •  PFf' MP AU J  LI T  Y.  F-TSO  •>.  2.X  ♦■•'C.P  I-T-tf.CAL 
1  •?  AO  I  US  *  IIS*  •  ♦  *?X  ♦  ♦  HP  1  OHT  •  •  3X. •4PFA*.aX.  ‘LIOUIO  FL.OW  PATH';/) 

.  N=?:0'0 

.xxr;x=KANF  ( r.?7a2  )■ 

on  io.  i  =  itN 

*'}!<*- (0Atjn(  0)  *SK)+AK 

FP-r.  ( r,  AI 1-3  ( .1  )  *3«  )  nt 

i?|i=(GAlJ~<0  )*SHHH 

'<AL  PP=  (  f.AUF  (  0  ) -Kf-HAL  ) +alp-‘> 

FAt..-(GAU'*»(C  >*SAL  >  +  AL 
A<,,=PK~i'AL/,VALf'c' 

:,f  2-r!^  ir*  { pH/1 )/r  tcpm 

GFAC=C»r  1  — ■"r .? 

0  s  A  P  =  A  M.  *  A  K  P  *  C,  F  A  C 

on  at- or  Af  *r 

?ir-  I T  F  <  C>  .  2-5  ? )  ',ir<AP.  pk.op  .ph.kal.palpp 
nv?vAT(.>x,ipr-:i  t  x  ,ipcir.5.6x.,  ipri  iPiiiz.n. 

.iox,  mr-  i  ?.r>«f.x,  iPr‘i?.!)) 

X  (  I  )  =  OBAR 

cr.Nnui:. 

PM  NT  (.0  0 

F,*>"  V  A  T(  •  SA*'Fl.L  31 7~  is*) 

i  r:  t  (.  \  o .  I*  ^ ^ 

rO’-’MAT  (  1  H*  1  i 

0-  100  I  <2-1  .N 
3Uv  =  sw.'+x<  i  2) 
fiX.l  4F-  >HM/.T 
SUV'.C=0  . 

DO  100  T.J-1.N 

SUMS  ,J=T.'JM  30  +  {  .X  {  T  •l)-'rXnAr?  )**2 

«uv‘*cn-oum‘:q/(.t-i  ) 

so?  ri=sor--T(  su.-ir.0N) 

Of.  I \T  "■»:>,  tlxnAH.sc?  I G 

ruPVATl //•  Mr  AN  VALtJ*'  ="•  .  T?  )  .3//.  *  STANDARD  DFV  I  AT  T  CN=  ♦  ,  F 20  .0// ) 
FF  =  SUMS  i/  7 


x?-n, 

X'i-0. 

X  =  o . 

on  ?o  t  ~  i  *  r j 
xi  =x  l  +>:  ( i ) 

X2  =  X  ;+X(  I  ) 

X  3—  X  3  T  X ( I  ) **3 
xi.  =X4  *y  (  I  )  *  *4 


m 


POC  does  not 


C M'3= X3/7- 3* (  XZ  ft  );*  <  x  1/7)+?*  c X 1  /Z) **3 

ALPHA=CM3/CS;QFiT  CrF*3.)) 

■PH  I  NT  6  66  e  ALPH  A3 

666  FORMAT.?//* COCFFI Clfc&T  OF  SKSWNFSS  =*  ,F20iQ//,> 

CM4=X^Zr4*<X  1/Z  >*  (  X3/Z  >  4A*  (X2/Z  )*{  X  l/Z  ) **2-3*1  X.i/Z>**4 
ALPHA4=CM4/(FF**2> 

PR  HIT  300, ALPHA* 

300  FORMAT?//* COEFFICIENT  OF  KURTHSIF  ,F20.8//> 

C  2=RM0F/,(  S  !  GMA*  12.)  ’* 

C3=?4. 

Q=C*AM-*(AK*AL/ALPP)*CC3/RE-C2*H> 

' nR INT ,0 
A4-=AK*AL 
A5=A4/ALPP • 

A3=03/R T— C2*M 
CALL  QUCTCRl ,0,»SK*1 ..RF> 

S 1  =  S  l  *C  3 
S2=C2*SH 

CALL  SUMO  IF C S3* SI ,52) 

CALL  PROD (S4,CK, PAL, AK, AL> 

CALL  0UCT(S5,S4 .SRALPP, A4 , ALPP) 

CALL  Pr.on{Sf.,S5,r>3,AS.A3) 

S6-=C*ANL*S6 
PR I NT. 36 
QB=C3/RE-C2'*H 
PQWA=C*ANL*(AK/ALFP)*OU 
PO WK=C  *ANL* ( AL/ALPP ) *OU 
PQWLP=-C*ANL*<  AL*AK/(ALPP**2>  >*Q0 
OC=C*ANL* ( AK*AL/ALPP) 

PQWH=QC*<-C2) 

PQWRE=QC*(-C3/(RF**2)  ) 

57={P0VJA**2)*(  SA,L** 2  )  +  {P0WK**2)  *<SK**2>  +  (PQWi;P**2  )  *  (SR  Al.PP**2  ) 
l+(PQWH**2)*(SH**?)+(PQVilF**2)  *(SR**2> 

S7=SQRT( S7) 

PRINT, 37 

CALL  RLCL5S (FX6ARtSQ3IG»  ALPPN  *  SRALPN  , H,  SH) 

Al PP=ALPPN 
SRALPP=SRALPN 
PR  I  NT  ,  SRAI.FP 

CALL  DCL  TAT(  PV  , ?-RV  ,C  XRAP  ♦  SOSI  G) 

4(.  0  CONTINUE 
65)1  CONTINUE 
CALL  FX IT 
END 


.  r  ' <•-  *  -f-  v/-'  ,,7-f. 


FUNCT  ION-  S  AUSCK) 

RANDCN  MIMHFR  GENE PAT»:n  PACKAGE.  RANrCK  )  GIVE*;  UNIFORM  RANDOM 
NUMBERS  *  GAUS'fK)  GIVES  .NOP  V  ALLY  PI  STRT-'nJTFO  (0*1)  NUMBERS. 

IF  iJA  l\F  HAS  OF  BN  TNTT  I AL  I  7cDi  GAUS  G WOULD  NOT  OF  T  MI  TI AL I  ZFD 


DATA  I  ;!"T  /0/ 

DATA  TW  OPT/  6.  2  S3  IB/ 

I F  {  JS'T) 10, 10,20 
ISFT-1 

A=SGFT(-2.TAl.OG(RANF(K>  )  ) 
B='?ANF(  0)*T'H0RI 
GAUS=  A*S IN(M) 

STPRF  =A^CP S  (P- ) 

RETURN 
ISFT=C 
GAUS= STOTT 
WFTURN 


FUNCTION  H  ANF ( N) 


C 

C 

c 

c 

c 

( 


c 


FORTRAN  ROUTINE  TC  GENFRATE  FLOATING  RANEjb.M  NUM3ER  3  OVFR  IHF 
.RANGE  OF  0. 

ROUT  INK  OY 

COU  I  VALENCE  (TX,RX) 

DATA  IX/32771/ 


o  to  i;.o. 

HARRYjp»  MIJFPHY « 


JP, 


22b 


<  f 

FEORUAPY 


1  <3*9. 


1  IF  (  I.  )  A  .  0  ♦  t 

2  TX=Mt)D(  10?I*IX+n»l  04TS7G) 
RANF^FLOAT  <iy)/l04°57£.0 
RFTUPN 

3  IX-VCC( N,  104  037c  ) 


A  RANF=RX 
RETURN 
FNO 


sun.^miTiNr  KSTsr(K,x,r:.xR4P,sosrG) 

OI  MFN<?T  ON  AAN'(400)  ,  F Vulo  )  ,GGF(  400  )  , D  IFF ( 400 > ,  X<  400) 
00  36  1=1. N 
56  FI  i  )  =  C. 

L  =  N  + 1 
M- N  + ! C 

00  114  I  =L . M 
4  X(  1  )  =  G. 

U=UX!1AP 
S  =  r,GS  IG 
PRINT  777, N 

77  FORMAT  (  /  7  •  S  AMPLF  SIZE=*.I5> 

7--t.0Ot+10 
7=0. 

5  DO  ?  1=1 .N 

1  r  (  X  (  I ) — Z  )  1 ,2 .2 

1  K=I— 1 
X  (  K  )  =  X  (  1  ) 

X{  I  )  =  Z 

Y=  1  . 

2  Z  =  X< I ) 

IF  (Y)  :» .6,4 

i  GO  TO  -3 

.  CONT  I  MJ!I 

0=0  . 

DO  S  1  1  =  1-,  N 

0=  1  . 

K=  I  +  1 

IF  <X(I))  60,51,67 

.2  IF  (  x(  I  )  — X  (  K  )  )  31  ,32,71 

17  00  S’*  .J=I,N 

l_  =  J  +  l 
L2= J+? 

^  X(L  )  =  X ( L 7  ) 

0=0+1. 

0  =  0+1  , 

!F<Y<I)-X(K)>  ^  1.2 7, 7  1 
.5  1  CONTI  VJO 
F<  I  )  r-  0 
1  COvT  l  n* JF 
L  =  o 
M  =  N— L 
0  F  7  =  0  . 

00  71  I = 1 , M 
r  =  \ 

OB  cf=F (  I ) 70 
Or  7  =o H z  +0  isr 
7  1  001-  (  I  )=CF  Z 

DO  41  I=1,M 
UU=(X(I  )-U)/S 
R  =  0  . ' 

IF  ( UC)  171,171,177 
131  UU=AFO(UU) 

R=1  . 

32  CONT  IN' if-! 

Z'.j=  1  ./(OQOT  (6.28217  )*FXP(UU*  +  2/7,  )  ) 

T=! ./(!.+• 22?67*UU> 


PX  =  1,V-  (  4  3<5T32ie*f-  .  1  20  t  ft7.*  (  T**2  )+  i<3372<3£d*{  T**3)-f  *  ZU 

IF  <R)  1  3.4  i  1 34  VI 35  •  '  ■ 

1.35  PX-  i  •, 

134  CONT  I  Mil; 

41  A  AN  (  I  )  =:,X 

DO  3 1  1=1, M 

.0  I F  F  {  T  )  =  AF!S<AAN(  D-fR^U  )  >. 

31  CONTINUE 
PRINT  333 

3HR  FORMA  7(  'OBSERVCD  DATA  FP.TOUcNCY  OOSERVFP  FUNCTION 

1  NORMAL  ^UNCTION  41 SOLUTT  O  I FFCPfc’NCE »  1 

on  i=i,*/ 

06  RR  I  NT  000,  v<  I  )  ,F(  I  )  ,00f-{  I  )  ,  AAN(  I)  ,0  I'-FI  I) 

090  FORMAT  (  1  H  ,  T 1 3 . 3  ,  SX  ,  F  1  0 . 3 , 5X  ,  <=  l  5.  o  ,  5  X  ,  F  l  5  .  P.  ,  5X  ,F  1  5  .■ 3 ) 

OIFVX=0. 

no  32  I  =  1 , M .  — 

•  IF  (DIF  MX  — D  IFF  (  I  )  >  36.33,  ‘V*  Reproduced  from 

36  D I F  v x  sf)  I  f F  (  I  )  best  available  copy« 

.3  3  CONTINUE 

PRINT  4  00  *  D IF  MX 

A00  FORMAT! //‘MAXI  MUM  ABSOLUTE  n  FFfRCNCF- * ,F20  »S) 

TCST  =  1  •  .‘>3/  (  SORT  {  B  )  ) 

PRINT  4  33,  T t: 3 T 

401  FORMAT  (  /  /  'KOLMOGOROV  CRITICAL  VALUE  AY  9,->  PFWCT  CL  =  '.F15.9/> 
RETURN 
CND 


m 


-  ■ 
,  u'DA‘v-^Ci^' 


1  rr 


I 


B7 


sun  pout-inf  oclt'at (rv,  s ov so’) 

D’l  Mfe-N  *  I  \"1N  f)  T  {  1  2  6  } 

CnMMCK  /CPl.K/REGESL  «  SRFSFL 
Cl=f:Qt,/l7.S7 
Q=C1*0 
SQ=C1*5C 
RKW”. 206766 
RO  I -.41/. 303 
N  —  1  C  C 
RK»=2in. 

RLE=14.S/12. 

RI.C--7  «/  12. 

DO  5  I-  l.N 

RRV=(GAU3(0)*5RV j+RV 

RVC=('.  I>93-P0  V  )  *(  (  GAUS(  0)*  .0  IS  )*•  25)  +PRV 
«V  I  VC-  ..5R3/RVC 

RT  VER  =  (  ,  303-PR V)  *  (  (G  VJF  (0  )  *  .0  3) +  .41  )  +  RPV 
PI  V^R  =  .:302/RI  V=R 
Rb*CL=(GAUS(0)*SPFS°>L)  fRE'-EGl. 
RE-sr-Ai.or,(nive)/(6.2,ms)-i='3Kw-*(WLF;-pna.) ) 

R"SFW  =  M.OG(  RIVER  )/(  6 ,2 0  1 1  ?-*PKW*Rf:CL  ) 
RIVC=.3P.3/RRV  - 

F?0-(GAUS(  0)*5O)+O  ^ 

KrtPF*ALOG(ROI  )/(6. ?931P*'-IKP*RLr;)  ' 

n R  w“  =  R'-  SH  P* RF  SE/  C  RFSER  +  RF^F. ) 
o R  WC  =  Al. PC  (  p  T  VC. )  /  (  6  .-2  2 3  1  8 7 P •<  **P LC  ) 

RRPG=AL0G  (  ROT  )/  (  6  .2  33  1 5-*PKP*RLC) 

«:U><R=RRDF  +  RRWir+nRWC+,!RcC 
DT  (  I  )=R0*SUMP 
CALL  wnAM(S'.DT*OnT,SDT) 

<>H  t  NT  lOO.DRT.GPT 

FORMAT (  //  '  MEAN  OT  =  1  ,F  lc.«  3//  1  »  STOV  D  T-  *  . F  1  5 .  0// ) 
RETURN 


l<$ 


suoroutinf  akval<n.x. al.sai_.ak.sk) 

DIMENSION'  X  (  40.0  ) 

PRIM  501- 

1  FORMAT (  ‘DATA  ISM 
DO  1 A  J=1  ,N 
PRINT  621. X(J) 
l  FORMAT! lH,ri5.3> 

PRINT  16 

FORMAT!  »SAMPLF  SIZE  TSM 
PRINT  610.N 
0  FORMAT!  1H  *  I  5 ) 

DO  23  1=1. N 
A=(GAUS(0)*5AL )4AL 

3  3  X<  T  )  =  6.  1  3  2C—  04*(  2.0/X! J ))  *  !  2 . 205F-0  3/ A )  *  (  1  */C2  .265  >  *  (  P.  ,/6?  .265 

1  (  1 ./12.2  > 

Z  =  N 
SUM=0. 

00  101  12=1 , N 
1  SUM=SUM+X( I  2) 

FXCAR=SUM/Z 
SUWSC=0 . 

DO  1<J  1  T3=1.N 

L  SUMSOr'.WMSO+IX  <  I  3  1-FX8AR  )  **S 
SUM  =  QN=-iiJMSO/(  Z-l  ) 

SOSIO-Opp T( SUMSON) 

PRINT  106..  rXP.AR,  SC3  IG 
:•>  FORMAT!//*  vr  AN  VAL'JF  =  *  «F2C.0//»  'STANDARD  RFV  T  AT  T  0N=  •  .  R20 . 8// ) 
AF  =  SUMSQ/ Z 
X1=0. 

X2=  0. 

X3=0  . 

X4“0  . 

DO  ?  1  ! =1 ,N 
X 1 =X 1 +X  !  I  ) 

X2  =  X2+  X ( I  )  **2 
X3=X2+X(T  >  *  *  3 
X4=X4  +X  (  I 

CM3-X?/Z-3*( X2/7) *<X1 /Z)+2*(X1 /Z)**.3 
Al.  PH  A. 1=  CM  3/  !SQ7T  (  AF  <  *3)  ) 

PRINT  2  03  .  ALPH A3 

3  ^ORMA  T  !  //  '  COPFF  1 C  I  FNT  CF  SK^WNF-SS  =',F?0.8//) 

CN'a=XO/Z-A*(X  1/Z  )  *<X3/Z)46M  X2/7  )*!  Xl/7  )  **2-7«! X 1/Z  >  **4 
ALPHA4=CM4/ C AF**?  ) 

PRINT  30 1 . ALPHA A 

1  FORMAT!//'  COEFFICIENT  CF  KUHTOSIF.  =',r?0.a//) 

AK=- X04" 

SK=SQStG 


RETURN 


SUBROUTINE  RVALU IN»X»RE,SR) 

DIMENSION  X ( 400 ) 

PRINT  000 

00  FORMAT!  'DATA  .!«»*) 

DO  -13  .1=1  »N 
3  PRINT  620,X(J> 

620  FOl,MAT(  1H,F15»8) 

PRINT  IS 

5  POPMAT(  •  SAMPLE  STZE  IS*) 

PRINT  (>  11«N 
ill  FORVATt  1H  ,  IS) 

C0NST='>83.*4.0234f:r'-03/6? .  ?3fjO 
DC  162  1=1  iN 
62  X(  I  ) =OONST/X< I ) 

Z=N  ' 

SUN  =0 • 

DO  102  I  2- 1  *N 
82  SUM=SUM+X( 12) 

EXOAR. -SUM/ 7 
SUM SO  =  0  . 

DO  IP  2  13=1  ,M 

02  SUM£'0  =  3U^S0+(X(  I  3  )-CXE  AR)  **? 

SUM  SO  Nr:  SUM  SO/  (  7-  1) 

SOS  I G=S  ORT  <  SUMSON ) 

POINT  1  C7.EX0AP  i  SOM  3 

07  F  0  R  M  A  T  (  /  /•  •  M  P  AN  VALUF  =  *  «T  20  ••)//*  'STANDARD  DC  VI  A  T  I  0N=<  *F?  0 . 3// ) 
AF=?UNSC/7 
X1=0. 

X2  =  0  . 

X3=0. 

X  A  =  0  « 

no  ??.  i  =i  .  n 

XI  =  X1  4X{  I  ) 

X?=X2  +X(  I  )**2 
X3=  X3  +  X ( I  ) **3 

??  XA=X44X( I )  **4 

C.M3  =  XJ/Z-3*(X?/Z)*(Xt/7)4?*<Xl/7.  )**3 
ALPHA  3=CM3 /(  SORT  (  AF**3)  ) 

PRINT  204, ALPHA.? 

?04  FOTMATI//'  COEFFIC  TENT  CF  SKE'aNEfSS  =  •  ,  F20 . 8/ / ) 

CM4-X4/7-4*  <  X  1  /7  >*(  X3/Z)  4f>*  {  X2/7)*  (  Xl/7)  **2-7*  (Xl/7)**4 

ALPHA4-CM4/I AF**2 ) 

PRINT  302  ,  Al  PHAA 

302  FORMA  T  (  //  *  COrrrF  I  C  TFNT  OF  KLPTOSIS  =*,r2C.O//) 

SR=SCSTG 
pf=exrar 
RFTURN 


Reproduced  from 
best  available  copy. 


SURRCUT  INE  MV  A LU ( N  »  X  v  AL  '*  3  AL  ) 
DIMENSION  X  {‘400  ) 

PPINT  600 

FORK AT ( 'SAMPLE  SIZE  IS*} 

PS l NT  6in,N 
FORMAT ( 1H» T6) 

PI= j. 14 t5P/F76 . 

r)l-.7£6 

DO  30  I=1,N 

X(  I  )  =  PI*(OI*-fc?-X{I)**2) 


SUM=0 . 

f^O  IfO  12  =  1.  N 
SUM=SUM  t-X  {  T  S) 

EXPAR=S  UM/Z- 
SUXSC=0. 

DO  too  13= 1 iN 

fUM«Q=SUMTO+(X{ 13 )-EXEAP 1 **2 
SUMSON=SUMvO/(7-l) 

:SOS  rO-SODT ( “OMSON ) 

PRIM  ViS  .EX'  AR.SOOIC 

FOSM  A  T(  //  *  MEAN  VALUE  ='  »e?C«3//.  '  STA  ND  A'RD  DFVI ATI 0N=»  .F20 . 3// ) 
F FsSUVSO/Z 
X1  =  0. 

X2=n, 

**=<*• 

no  20  i -l.-i 

XI  =X1  +X (  T  ) 

X2=X?*X(  I  )-'*+■} 

X3=X34X  (I  )**3 
X4  =  XA+X{  T  )  XAA 

CM3=X  1/7-  •»*(  X2/7  )*(X1/Z)+2*<X1/Z)**3 
ALPHA=r.v3/ (SORT  (FF*3)  ) 
os  I  NT  on6.4L.PHA3 

~\}RM  AT  (  //  ’  vIOFFF  I  C  IFNT  OF  SKFWNFSS  =*,F20«8//> 

0*4  =  <4/7-./i*(Xl  /?)*(X.3/Z)+6*{X2/Z>*(Xl/7  )  **2-3*(  XI  /Z  )  **4 
ALPHA4=EM4/(Fr**2) 

"PINT  300  ,  ALPHA  A 

r  OSVATf //«  C*'F  FFICI7NT  OF  KUFTCSTS  =*.F20.8//> 

SAL  =  SOS  16 
AL  =  rX  LA  ~ 

PE  TUP  N 

FMO  Itfd  . 


w  -  r-bto  ** 


600 

610 

i 

180 

1*30 

?5<5S 


SUBROUTINE  M  E  A  N  {  N  ,  X  .  F  X  8  A  R  ,  S  G  S I G  i 
DIMENSION  X(400) 

•  PP  I  NT  600 

FORMAT! ’SAMPLE  STZF  IS*) 

PRIM  610.N 
FORMAT!  IH,  I  5> 

■Z=N 
SUM=0. 

00  180  I?=1*N 
SOW -SUM  +X  (  12) 

EXHAR=^UM/Z 
SUM SO =0 • 
on  100  T  3= 1  ♦  N 
SUMSO  =  SLMSO+(  X{  13  j”EXQAR)**2 
SUMS0N=-SUMSQ/(Z-1) 

SQ  SI G  =  S0RT ( SUMSGN ) 

"PINT  B'3S»  FXDAP  *  SCSI  G 

FORMAT! //• mpaN  VALUE  - * « F2C • 8//« * STAMDARD  OF V  I  AT J PN= » » F ?0. 8/ / ) 

RETURN 

END 


(4$ 


Copy  ovcnMbW  to  BIX’  do03  nB} 
•>  f-’iy  ]-e  Mo  reproductions 


B12 


5UBR0UTINF  P.ECESSlOi  SO.  ALFPN.SRALPN.H.SH) 

DIMENSION  Y J 200) 

COMMON  RHQF,  UL«  AL  »H  IL  «RE  »  SR  »HFG  *  AK  •  SK  r  SAL. 

COMMON  /CeLK/RECFSL.SRCSSL 
AMO={ G*.056Q2/{HFG*60. ) ) 

SAMO=SQ*. 05692/1  HFG*60.  ) 

EL  =  14  .5/12. 

SEL  =  .  03 
FH=F/12. 

SFW=S>-/12. 
x=n  o 

0X=. 25/12. 

Cl=( LL*.5*AW0) /( (FH0F**2) *32.2*AK*AL >+{9.5/24 . }*FH 
A=(  AMO*UL)/<  AK*32.2*(RHnF**2)*EL*3.*AL/4..  ) 
fl=FH/-2. 


5 


3 

5 


w 


l  0 


N  =  5S 

DO  ?  1=1 .N 
X=X+px 

P0=-A/2.*(  X**2)+(EL*A+BJ*X+C1 
IF  ( ° C—  F 1 L  )  5,6.5 
OPO  =PD 
INT.PO 

nx=x 


cont  r  mje 
RFC.FSL=CL-OX 

AL f  F N  ='1»“C  E 5L*  *  2/  (  5  .-2 5  +£:  L  >  +  2  .  *El  /  3 .  +  .  5 

NN=  100 

DO  9  1=1, NN 

RK=  (  G  AIJS  (  0  )*5K)+AP 

P  A= ( GAUSI 0)*SAL ) +AL 

RAMO= (GAU5 (0 )+SAMO)+AMO 

RH= (G AUS( 0)*SFH) +FH 

WTH=<GAUS<  0)*.G52>+.75 

RC=  (UL'X.5*RAM0  )/(  (RFCF**2  )  *22  .?*RX*RA >  +  ( O. 5/24  .  )*RH 
A  A  =  (  P  AM 0 * Ut.  >/{  RK*2.2. 2* ( RHQr ** ? > *FL* R TH* P A > 

R0=RF/2. 

RPP=-  (A A/?  *  )*<  9X**2>  +  (FL*AA+R3)*nX  +  RC. 


Y(  I  )— FPO 
CGNT  INUC 

C  A  1  L  ^ANINN.Y.XOAP.SX) 

PD=  yw  A  l 

SPD=SX 


Reproduced  from 
best  available  copy. 


SLOPS =— A*OX+FL*A+D 
RB-SPO*  *2+SH**2 
SRF??L=cQRT(OP)/SLCPF 


00  10  I  =  1 , N  N 

RRFCSL~(  GAUSfO  )  *SfCESSL  )  +RECESL 
R  TH  RD  =  (  G  A  U  S  (  0 )  *  .  0 1 4  )  +  .  1  n 
REL=(CAUS(  0)*SEL  HCL 
RTTH=(GAUSCO>* .04 R  )+2 ,/3 , 
RALPPN=(RRFCSL**?)*FTHR0/RFL  +  PTTH*PEL+.-5 
Y (  I  )  =  RALPPN 

CALL  WFAN(NN,Y,X3APt«X> 

SPALPN=SX 
PRINT  ,3RALPN 
PRI NT  ,PX 
PR INT  ,SRFSSL 


RETURN 

END 


I 


B13 


Copy  cvc'i.ublt  to  DLX.  does  w>t 
peii.ji  tu'ly  kg.Lio  reproductioa 


SUBROUTINE  SUPDIFfSA.SE.SC) 

SA=£B**2+SC**2 

SA=SQRT(SA) 

RETURN 

END 


SU3PCUT INE  CUD7<SA»S3,SC.EEAR,C3AR) 

SA=( {CBAR**2)*  {SH**2 )+t 8BAR**2 )*?SC**2) )/{CBAR**4) 
SA=SGRT(SA> 

RETURN 


SU»3RCUT  INE  PROC  {SA . SB . SC .80 AR .CBAR ) 

SA-C  RBAR**2 ) *  < SC  442)  +  ( CB AE**2 ) * (SB**2 >  +  ( SD**2 ) *i SC**2) 
SA=  SORT C  SA  ) 

RETURN 

CN0 

1  best  available  jgpjNJgZP' 


5UORQWT INC  XSOAPCI SX2.SX.XDAR) 
SX2=4 .*X8AR*XnAR*SX*SX+2 .*SX**4 
SX2=SOPT( SX2) 

RETURN 

END 


Cc'.yy  cv  tr.  <?OS3  tl6f 

pc-ii..,!  tu  1  f  ].;  ‘  rc;;.:^cUofl 


APPENDIX  C 


HEAT  PIPE  EXPERIMENTAL  TEST  PROCEDURES 

The  following  is  a  step-by-step  procedure  for  the  testing 
of  a  water  heat  pipe. 

1.  The  heat  pipe  is  thoroughly  cleaned  using  trichloro¬ 
ethylene,  alcohol,  and  water. 

2.  A  wick  is  saturated  v/ith  rresh  distilled  water  and  in¬ 
serted  into  the  pipe. 

3.  The  wick  is  inserted  with  the  seam  facing  down  and 
is  forced  against  the  evaporator  plug. 

4.  The  pipe  is  then  filled  with  distilled  water  to  fur¬ 
ther  saturate  the  clearance  area  between  the  wick  and  the 
pipe  wall. 

5.  The  water  is  drained  and  30  c c  of  excess  water  is 
injected. 

6.  The  condenser  plug  is  inserted  and  sealed. 

7.  The  pipe  is  then  evacuated  to  23  inches  of  Hg. 

8.  The  evaporator  is  lowered  and  full  power  (100  watts) 
applied.  The  start  up  is  observed  to  avoid  premature  burnout. 

9.  When  the  pipe  pressure  reaches  0  psi,  the  condenser 
valve  is  opened  to  vent  any  non-condensible  residual  gases. 

10.  After  venting  and  resealing,  the  calorimeter  is 
turned  on  and  the  system  allowed  to  reach  an  equilibrium 
at  the  particular  power  level  to  be  tested. 

11.  The  evaporator  is  raised  in  increments  of  one  fourth 
inch  until  a  wick  burnout  is  attained.  A  burnout  is  defined 

M 


as  a  sudden  temperature  rise  of  more  than  10 °F  at  the  extreme 
end  of  the  evaporator  after  equilibrium  is  reached. 

12.  A  second  measurement  is  made  to  determine  the  burnout 
height  within  one  eighth  inch  of  the  true  height. 

13.  Data  recorded  at  this  time  are: 

a.  Heat  pipe  temperature  distribution 

b.  Calorimeter  flo w  and  temperature  rise 

c.  Pipe  vapor  pressure 

d.  Burnout  height 

14.  The  power  is  now  disengaged  and  the  evaporator  seal 
broken. 

15.  The  gravity  head  is  measured  at  the  various  burnout 
heights ,  as  shown  in  Figure  3.5. 

16.  The  permeability  is  measured,  as  described  in  Figure 
3.9.  Thirty  flow  readings  are  recorded.  Permeability  is  cal- 
culated  according  to  the  formula  K  =  uL'm/pAAp,.  where  Ap  is 
assumed  to  be  the  linear  head  loss  through  the  wick. 

17.  The  wick  is  removed  from  the  pipe  and  portions  of 
it  are  placed  in  the  wick  apparatus,  as  described  in  Figure 
3.12.  If  five  consecutive  identical  readings  are  observed, 
they  are  considered  to  be  the  mean  and  the  standard  deviation 
is  chosen  to  be  one  half  the  least  count  of  the  instrument. 
r^r  is  calculated  as:  r^r  =  2o/pGH. 

18.  The  data  are  compiled  for  reduction. 

The  procedure  for  wire  mesh  wick  manufacture  is  as  follows. 

1.  The  raw  wire  mesh  is  cut  to  size  and  the  retainer  rod 
is  welded  to  an  edge. 


C2 


2.  The  wire  mesh  is  inserted  into  the  mandrel,  wrapped 
the  trapping  apparatus  (Figure  3.15)  and  welded. 

4.  The  wick  is  cleaned  as  follows: 

a.  Ultrasonic  cleaner  with  alkanox. 

b.  Rinse  with  tirchloroethylene. 

c.  Rinse  in  methanol. 

d.  Rinse  in  distilled  water. 

5.  The  wick  is  oxidized  in  air  at  850°F  for  two  hours. 

6.  Storage  is  under  distilled  water. 


